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Pioneer V 
Paddlewheel Planetoid 
Is Vaulting 


‘Through Unexplored Space 
Toward The 
Orbital Path of Venus 


At this moment Pioneer V, one of the most advanced space 
probe vehicles ever launched, is on a course toward the path 
of Venus—26 million miles from earth. Blasted aloft March 11 
by aThor Able-4 rocket booster, this miniature space laboratory 
will reach its destination in about 130 days. 

The project, carried out by Space Technology Laboratories 
for the National Aeronautics and Space Administration under 
the direction of the Air Force Ballistic Missile Division, may 
confirm or disprove long-standing theories of the fundamen- 
tal nature of the solar system and space itself. 

Energy from the sun—captured by almost 5,000 cells 
mounted in the four paddles—is used to supply all of the elec- 
trical power to operate the sophisticated array of instrumenta- 
tion packed into the 94-pound spacecraft which measures only 
26” in diameter. 


SPACE TECHNOLOGY LABORATORIES, INC. 6p 


Los Angeles e Santa Maria e Edwards Rocket Base e Cheyenne 
Cape Canaveral e Manchester, England e Singapore e Hawaii 


May 1960 


By combining a phenomenal digital electronic brain (telebit) 
with a powerful radio transmitter inside the satellite, STL scien- 
tists and engineers expect to receive communications from 
Pioneer V at their command over interplanetary distances up 
to 50 million miles. 

STL’s technical staff brings to this space research the same 
talents which have provided over-all systems engineering and 
technical direction since 1954 to the Air Force missile pro- 
grams including Atlas, Thor, Titan, Minuteman, and related 
space programs. 


Important positions in connection with these activities are now availa- 
ble for scientists and engineers with outstanding capabilities. Inquiries 


and resumes are invited. 


P. O. Box 95004, Los Angeles 45, California 
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motors of unprecedented size... 


Why THIOKOL subcontracts to 
Industry, U.S.A. in the 
production of rocket powerplants 


of unprecedented reliability. _ ag 


Indevelopmentof advanced power- 
plants for missiles like Minute- 
man, Nike-Zeus, Subroc, and for 
research vehicles like Little Joe 
and X-17...THIOKOL draws on 
its own vast propulsion know-how 
plus the advanced technological 
background of scores of industrial 
organizations. 

The Allison Division of General 
Motors, Scaife Company, Goodyear, 
General Electric, Heintz Mfg. Co., 
Curtiss-Wright, RCA, Solar Air- 
craft, Borg-Warner Corp. . . . these 
are but a few of the many compa- 
nies, large and small, to whom 
THIOKOL has subcontracted in 
producing dependable propulsion 
systems. 

We have called upon the pres- 
sure vessel industry to whom met- 
als for strength are second nature 
to get rocket casings combining 


THIOKOL’s Utah Division is fully staffed and equipped to produce solid rocket 
of ICBM and even satellite proportions. 


7 


light weight and high tensile 
strength. 


We have called upon the elec- 
tronics industry whose art is 
instrumentation for the delicate 
devices required for precise test- 
ing and production controls. 

We have turned to the transpor- 
tation and construction industries 
for development of specialized 
equipment such as monorail sys- 
tems and movable cranes needed 
to process giant rocket motors 
with unfailing precision on an 
assembly line basis. 

Many industrial technologies 
are met in a rocket propulsion sys- 
tem. Recognizing this, THIOKOL 
calls on specialists to achieve high- 
est reliability, to meet the critical 
rocket power requirements of na- 
tional defense and space research. 


Vheokol. Chemical Corporation 


BRISTOL, PENNA. 
Plants in: TRENTON, N. J.; MOSS POINT, MISS.; DENVILLE, N. J.; ELKTON, 
MD.; HUNTSVILLE, ALA.; MARSHALL, TEXAS; BRIGHAM CITY, UTAH. 


®Registered trademark of the Thiokol Chemical Corporation for its liquid 
polymers, rocket propellants, plasticizers, and other chemical products. 
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At The Wooldridge Laboratories. 
integrated programs of research & development 


of electronic systems and components. 


The new Ramo-Wooldridge Laboratories in Canoga Park provide 
an environment for creative work in an academic setting. Here, 
scientists and engineers seek solutions to the technological prob- 
lems of today. The Ramo-Wooldridge research and development 
philosophy places major emphasis on the imaginative contribu- 
tions of the members of the technical staff." There are outstanding 
opportunities for scientists and engineers. Write Dr. Richard C. 
Potter, Head, Technical Staff Development, Department 11-E. 


$4 THE RAMO-WOOLDRIDGE LABORATORIES 
® 


8433 FALLBROOK AVENUE, CANOGA PARK, CALIFORNIA 


An electron device permits scientists to study the sehavior of charged dust particles held in suspension. 


not 
and 
mers 
inclu 
cylin 
ment: 
are Sl 
third 
Al 
been 
such ; 
absor] 
tion is 
to be. 
viseou 
Stokes 
effect 
condit 
velocit 
genera 
plied t 


May 1! 


2).2 
in tl 
rem: 
phe 
oi a 
The 
ratic 
dime 
disti 
raref 
va 
tm 


6" MMARIES of —" developments in the mechanics 
J of rarefied gases through 1956 are contained in (1 and 
2).2_ The present paper thus will concentrate mostly on work 
in the field since that date, except for a few brief introductory 
remarks. 

The general field of rarefied gasdynamics is concerned with 
phenomena related to the molecular or noncontinuum nature 
of a gas flow which may occur at sufficiently low densities. 
The relative importance of such effects is determined by the 
ratio of the molecular mean free path A to a characteristic 
dimension I of the flow field. It has been found convenient to 
distinguish between flows which are slightly rarefied (A/l ~ 
0.01 — 0.1), moderately rarefied (A/l ~ 0.1 — 10) and highly 
rarefied (A/l > 10). These are usually referred to as “‘slip 
flow,” “transition flow” and “free molecule flow,’ respec- 
tively. 

There is by now a rather extensive literature—mostly, but 
not entirely, experimental—dealing with aerodynamic force 
and heat characteristics of a variety of geometrical shapes im- 
mersed in a low density flow. Geometries so far investigated 
include: Flat plates, spheres, cylinders, cones, wedges, cone- 
cylinders, sphere- nosed cones, and some special develop- 
mental configurations. The data of this type available in 1956 
are summarized in (3); more recent data are included in the 
third section of the present report. 

A number of additional experimental investigations have 
been carried out for flow situations of basic scientific interest, 
such as Couette flow, shock structure, and the dispersion and 
absorption of ultrasonic sound waves. This type of investiga- 
tion is summarized in (2,4,5 and 6). The general conclusions 
to be drawn from these investigations are: 

1 The “slip flow’ regime is dominated by continuum 
viscous effects describable in terms of the normal Navier- 
Stokes formulation. The only important noncontinuum 
effect appears to be a slight modification in the boundary 
conditions for the gas flow past a solid surface, where slip 
velocity and temperature jump phenomena occur. These 
generally result in relatively small correction terms to be ap- 
plied to standard continuum viscous results. 

2 The “free molecule flow” regime seems to be quite satis- 


Received March 11, 1960. 

1 This article is a modification and extension of a paper pre- 
sented at the Sixth Midwestern Conference on Fluid Mechanics, 
Austin, Texas, Sept. 1959. 

2 Numbers i in parentheses indicate References at end of paper. 
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factorily understood in terms of the kinetic theory of gases. 
The prediction of heat transfer and aerodynamic force charac- 
teristics for convex bodies may be carried out in a quite 
straightforward way (except for the question of molecule- 
surface interaction, which will be discussed in some more de- 
tail further on). Experimental measurements appear to pro- 
vide quite good agreement with the theoretical predictions. 

3 Nosimple, generally applicable, analytic method has yet 
been developed for treating the “transition flow” regime. 
Experiments as well as linearized kinetic theory, however, 
suggest that a simple interpolation between “free molecule 
flow” and “slip flow” results is usually possible. 

With these general remarks we now turn to a more detailed 
discussion of recent developments in the field. a 


Limits of Continuum Flow 


There has been much discussion recently [see (7,8 and 9) ] of 
the correct limits on the gas density to insure that a flow is 
essentially continuum. The original demarcation between 
slip flow and continuum flow was suggested by Tsien (10) as 
corresponding to a value of \/6 ~ 0.01, where A is the molecu- 
lar mean free path and 6 is the boundary layer thickness. 
(For very low Reynolds number flows, the boundary layer 
disappears and only the characteristic dimension / of the flow 
field itself appears; the discussion in the present paper is 
directed toward the higher Reynolds number case.) For low 
speed flow past a flat plate 6/1 ~ 1/- Re, where Re is the 
Reynolds number and AX is essentially constant everywhere. 
Hence, using the relationship that \/l ~ M/Re, one finds that 
1/6 ~ M, ~/Re. The original Tsien definition of the de- 
marecation between slip and continuum flow was accordingly 
expressed in terms of the Mach number / and the Reynolds 
number Re by M /V Re~ 0.01. Most of the recent discussion 
referred to previously has been concerned with the extension 
of this basic notion to the high speed case. For very large M, 
the situation becomes more complicated, due to the depend- 
ence of both \,, the wall mean free path, and 6, the displace- 
ment thickness, or 6 the boundary layer momentum thickness, 
on M body geometry and the heat transfer conditions at the 
body surface. It seems reasonable to suspect that slip effects 
on induced pressure are governed by the parameter \,./6% 
[see, for example, (11) ], whereas slip effects on skin friction are 
determined by X,,/@ (2). Stagnation regions on blunt bodies 
will have different expressions for \,,/5, and for X,,/8 in terms of 
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Dr. Schaaf, a graduate in mathematics of the University of California, Berkeley, is now Pro- 
fessor and Chairman of the Division of Aeronautical Sciences there. 
contributions in the field of fluid mechanics, particularly on the subject of rarefied gas flow. He 
is the author of the section on this subject in the Princeton Series volume ‘‘Fundamentals of 
Dr. Schaaf has served as chairman of the Fluid Dynamics Division of the Ameri- 
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Mach and Reynolds number than will the downstream side of continuum flow are given in terms of an appropriate value for t 

slender objects. Also, adiabatic boundary layers—frequently ../dy or for X,,./8, and it turns out that all of these cases re- b 

the wind tunnel case—will, in general, have expressions for sult in a demarcation between the slip flow and continuum re- is 

\,,./d, and for \,./@ in terms of Mach and Reynolds number gimes at lower densities than the corresponding Tsien crite- c 

which differ considerably from highly cooled boundary rion. No experimental disagreement with the basic ideas 0 


layers—usually the high speed free flight case. Thus a large has yet been noted. 
variety of criteria for determining the onset of slip flow effects The emphasis on the presence of noncontinuum effects, such: 7 
may be established at high M. In every case, however, the as slip and temperature jump as the defining characteristic o: 


suggested modifications in a demarcation between slip and the slip flow regime has perhaps been somewhat unfortunate, 
in that it has tended to obscure the very important fact that p 


the dominant fluid mechanic effects in the slip flow regim: re 
are not these noncontinuum effects at all. The dominan: 


effects are associated with the presence of the very thick fi 
p= 5° laminar boundary layers which occur and which give rise 
aap pr skin friction of the same order as the pressure, and which W 
eee © Wind Tunnel Thermal Conditions | through the displacement effect, disturb the surface pressur» H 
& @® Free Flight Thermal Conditions from its inviscid value, oftentimes manyfold. These ar: in 
e|s completely continuum effects and are basically well under- Jes 
or 5 stood, although their application to determination of aerody- Vi 
S| x namic coefficients presents considerable computational dif- 
>} 1000t— _ a ficulty. A simple example will serve to illustrate the point. he 
Consider the skin friction on a wedge in low speed flow. Th ex 
so boundary layer for the no-slip case is that considered tre 
Ls Falkner and Skan (12), and leads to an expression for th cul 
skin friction coefficient in the form du 
MACH NUMBER res 
: s Modifying the boundary layer analysis to include the effect 0! ais 
Fig. 1 Viscous effect on the axial forces of a wedge at zero slip (2, 13), one obtains 0 
angle of attack ] 
T T T V Re V Re 
——— £=30° inc =} : lhe slip correction is thus seen to be of the order of M/ / Re, blu 
in © Wind Tunnel Thermal Conditions as expected. More important, however, is the fact that for I 
a 200}-  @ Free Flight Thermal Conditions ara 4 this term to amount to a correction of the order of a few per vat 
Ie ip. . cent, it is necessary that Re be very small indeed, i.€., of the (20 
order of 10? or so. Under these conditions, the skin friction was 
3| 5 re with or without a slip correction amounting to only a few tha: 
212100 os J per cent or so—will be orders of magnitude larger than at nor- nar 
S mal density. Most of the experimental results so far ob- 
Z 7 tained in the slip flow range have reflected these continuum mes 
nae - viscous effects, rather than noncontinuum slip flow effects. in I 
sa The more recent of these are summarized in the next section. inte 
Sse re ae a If slip corrections to boundary layer effects are significant, A 
| | | D | then other corrections, for example departure from strict utili 
aan ona 15 20 boundary layer theory itself, must be considered. For the out 
flat plate case this is apparently very difficult [see, for ex- with 
Fig. 2 Viscous effect on the normal forces of a wedge at zero ample, (14)]; for the blunt body case in highly cooled flow the cinit 
angle of attack correction to boundary layer theory due to the vorticity mad 
generated behind the detached shock wave (15) appears to T 
dominate the slip effect, so that even lower densities are neces- plate 
' T T sary for this case before any sizable departure from con- heat 
tinuum flow will become of importance. However, a calcu- the 
lation of the “shock layer” heat transfer with and without a 
ne eens modification in the boundary condition to include the tem- Tr 
° @ Free Flight Thermal Conditions Tier 7] perature jump effect has been carried out by Hoshizaki (16). of th 
E He finds a substantial effect due to slip which serves to il- Thee 
al 8 lustrate the point that noncontinuum effects are often of im- on 
3] 2 portance when considering modifications of boundary layer resul 
= F100 = theory. In connection with the foregoing discussion, it is 
S [oS appropriate to emphasize the importance of both geometry Ki 
S, and heat transfer effects on the viscous modifications to aero _ 
a dynamic coefficients in a low density continuum flow. Thes« Se 
a > 4 coefficients depend not only on Mach and Reynolds numbers, deali 
but also on the heat transfer conditions at the model or proto- tion 
| | | Gf | type surface. Lift, drag and moment characteristics of simple Boltz 
mace 10 LUMBER 15 20 wedge-shaped airfoils calculated on the basis of a no-slip mary 
laminar boundary layer and Newtonian flow are presented in and [ 
Fig. 3 Viscous effect on the moment coefficent of a wedge at Figs. 1, 2 and 3 for two different heat transfer conditions, one for sc 
zero angle of attack corresponding more or less to wind tunnel test conditions, and the k 
444 | ARS JouRNAL May 
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the other to the free flight case. The very great difference 
between the two, particularly for the slender wedge airfoils, 
is quite striking. It is also apparent that viscous effects be- 
come of importance at higher Reynolds numbers on slender 
objects than on blunt objects. 


Recent Investigations 


A number of experimental investigations have been com- 
pleted during the past 3 years in the slip flow and transition 
regimes. These will be summarized briefly here. 

Lift and drag forces on a slender, sharp cone-cylinder con- 
figuration have been measured (17) at M = 4 and Re = 
1 X 10%. It was found that the measured drag coefficient 
was approximately 10 times the inviscid drag coefficient. 
However, the simplest laminar boundary layer correction, 
including both skin friction and displacement pressure effects, 
leads to good agreement with the experimental data (see 
Vig. 4). 

Pressure distributions on sharp leading edge flat plates have 
heen measured by a number of different investigators. The 
«xperiments which have penetrated farthest into the slip and 
transition regimes are reported in (18). These and other re- 
sults published in available publications will not be repro- 
duced here. It was observed that the data followed con- 
tinuum theory until rather high values of the interaction 
parameter, X ~ M?/+/Re, but departed at a point cor- 
responding to X,,/64 ~ 0.35 (18). This was presumably due 
to the effect of slip near the leading edge (11). 

Pressure distributions on sphere-nosed cones (19) have been 
determined at M ~ 4 and M ~ 6 at low Reynolds numbers. 
It was observed that—in contrast to the sharp leading edge 
configuration of the cone-cylinder and the flat plate—neither 
induced pressure nor slip effects were present, due to the 
bluntness of the geometry investigated. 

Heat transfer to cross-stream cylinders has been investi- 
gated by several groups (20,21,22). The results obtained in 
(20) extend farthest into the slip flow regime. Here again, it 
was noted that the bluntness of the configuration insured 
that the results were within the scope of applicability of ordi- 
nary boundary layer theory. 

Drag coefficients for cross-stream cylinders have been 
measured and reported in (23). The results are presented 
in Fig. 5, which serves to illustrate the applicability of the 
interpolation process for transition regime results. 

A study of the flow in the vicinity of a very short flat plate 
utilizing a “free molecule pressure probe” has been carried 
out in (24). The results were in reasonably good agreement 
with a simple diffuse reflection model for the flow in the vi- 
cinity of the plate, and constituted the nearest approach yet 
made to the direct observation of a slip velocity. 

The flow field in the vicinity of the leading edge of a long flat 
plate has been studied (25), utilizing a “free molecule flow 
heat transfer probe.”’ Various features of the flow, including 


the interaction between the shock and the growing viscous 
layer, were revealed in a qualitative manner. 

In addition to these experimental investigations, a number 
of theoretical calculations have been made for rarefied flows. 


These are mostly based on the Navier-Stokes equations and 
The 


on slip and temperature jump boundary conditions. 
results of most interest are listed in (26 to 32). 


Kinetic Theory Developments 


Several new methods have been developed recently for 
dealing with basic flow situations, particularly in the transi- 
tion flow regime, in terms of the fundamental Maxwell- 
Boltzmann equation of kinetic theory. An excellent sum- 
mary of these methods is given in (33). Gross et al. (34,35,36) 
and Uhlenbeck and Chang (37,38), have developed methods 
for solving the linearized Maxwell-Boltzmann equation, i.e., 
the kinetic theory of infinitesimal disturbances in a gas, but 
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with no restrictions as to the gas density. The essential dif- 
ference between the two methods is in the use of “full range” 
expansions, i.e., expansions of the molecular velocity distribu- 
tion function in terms of a set of eigen functions of the linear- 
ized collision operator taken over the entire range of molecular 
velocities, or in “half range” expansions, i.e., expansions of the 
distribution function which explicitly distinguish between 
molecular velocities directed toward or away from surface 
elements. Most of the simple basic problems, such as Cou- 
ette flow, the impulsive flat plate, heat transfer between 
parallel plates, shear flow near a wall and the dispersion and 
absorption of sound waves, have been investigated along 
these lines. Perhaps the most important conclusions from 
these studies have been: 

1 There are no unexpected phenomena in transition flow. 
The results usually constitute reasonable interpolations be- 
tween the free molecule flow and continuum or slip flow limits. 

2 Slip flow and temperature jump boundary conditions 
along with the full Navier-Stokes equations seem to be valid 
for predicting macroscopic flow conditions in the slip flow 


regime. 


3 The mathematical dependence of the solutions on 


seems to be nonanalytic in both \ and 1/), so that simple per- 


turbation expansion procedures for modifying either con- 
tinuum or free molecule flow results appear questionable. 

A somewhat more powerful method has been suggested 
by Lees (39). This involves solution of the Maxwell trans- 
port equations with suitable parameterized forms for the 
molecular distribution function. The method is in the 
spirit of the well-known momentum-integral boundary layer 
method, and also to some extent follows the earlier work of 
Mott-Smith. An entirely different approach to low density 
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kinetic theory has been formulated by Bhatnager, Gross and 
Krook (40). They suggest replacing the momentum collision 
integral term in the Maxwell-Boltzmann equation by a re- 
laxation expression which requires the distribution function 
to decay toward local Maxwellian equilibrium at a rate pro- 
portional to its instantaneous departure from it. This modi- 
fied Boltzmann equation (sometimes called the ‘Krooked 
Boltzmann equation’) is more tractable for analytical work. 
Solutions for typical problems have been in good agreement 
with results of the more exact methods (40,41). 


Near Free Molecule Flow 


The flow regime at densities which are just beginning to 
depart from strict free molecule flow has been the subject of 
a number of recent studies. The case of a sphere in very 
high speed flow, for example, has been investigated by Char- 
wat and Baker (42) among others (43,44,45). Molecules re- 
emitted from the upstream surface occasionally deflect free 
stream molecules away from the surface which they would 
otherwise have struck, as well as occasionally deflecting free 
stream molecules into the surface which they would other- 
wise have passed by. It appears that the former effect is the 
dominant one, so that there is a net “shielding effect”” which 
reduces drag and heat transfer coefficients if the density in- 
creases. The effect begins at surprisingly low densities for 
the high speed flow case. At a Mach number of 20, for ex- 
ample, the free stream mean free path must be a hundred 
times the sphere diameter to insure strictly free molecule flow. 
This result is thus in line with the preceding discussion, which 
indicates that high speed noncontinuum effects for blunt 
bodies begin at lower densities than had been originally sup- 


posed. 
gard 
Free Molecule Flow 


The aerodynamics of free molecule flow, i.e., flow at such 
low pressures that the molecules re-emitted from the surface 
have no effect on the incident flow, had been regarded as 
pretty well understood [see (2) for example] until very re- 
cently. However, the importance of and the consequent 
need to learn more about the exchange of momentum and 
energy between gas molecules and surfaces has only lately 
been realized, mostly in connection with the problem of pre- 
dicting drag forces on artificial satellites. The classical free 
molecule flow analyses obtain expressions for heat transfer 
and aerodynamic force characteristics in terms of momentum 
and energy exchange ‘accommodation coefficients” which 
range from the limiting condition of completely diffuse re- 
flection for which the exchange is complete, to the case of 
specular reflection for which there is no exchange. For some 
applications, particularly at very high velocities, it develops 
that the actual physical interaction is probably such as to give 
gross aerodynamic characteristics which are not intermediate 
to these two supposedly limiting types of interaction. The 
theoretical drag coefficient of a sphere, for example, in high 
velocity free molecule flow is exactly 2.00 for both the ideal 
specular and diffuse types of surface interaction. The theo- 
retical drag coefficient for a more reasonable model for the 
surface interaction now appears (46) instead to be in the 
range 2.2 to 2.6, depending primarily on the degree of dis- 
sociation experienced by incident nitrogen molecules upon 
collision with the surface. The need for more information, 
particularly in the satellite range of velocities, had led to the 
development of new experimental equipment in the form of 
high energy ion and molecular beams. No experimental data 
of aerodynamic use are yet available, however, in the case 
of satellite energies. At Berkeley, molecular beam (47,48) 
and rotating cylinder (49) investigations of surface-molecule 
interactions have been carried out. These studies have 
suggested a new approach to the “‘accommodation coefficient” 
method of free molecule flow analysis. The new method (50) 
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involves the use of three empirical constants to describe the 

surface interaction, but should be more in accord with the 

underlying physical reality than the older method. bod 


Symposia 
The rapid growth of the field of rarefied gasdynamics as a 
research area has given rise to a series of international sym- 
posia. The first of these was held in Nice in the summer o! 
1958; the proceedings of this meeting will be published 
shortly. The second symposium will be held Aug. 3-6, 1960. 
at the University of California, Berkeley. It will be preceded 
by a short course in Rarefied Gas Dynamics, conducted b) 
Dr. Schaaf, to be held July 18-29. 
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Drop Size 


7 graphic problems encountered in the development are discussed. 


FFICIENT power generation by fuel oxidation demands 
high rates of heat release within efficient combustion 
reactors, often of limited volume. In combustion reactors 
where fuel, or propellants, are originally present in the liquid 
phase, information on the rate and mechanism of the atomiza- 
tion and vaporization of these materials is of interest, partic- 
ularly in reaction propulsion when liquid fuels are used (1-4). 


Presented at the ARS 13th Annual Meeting, New York, N. Y., 
Nov. 17-21, 1958. 

1 Research Scientist, Radiation Laboratory. 

2 Professor of Mechanical Engineering. 

3 Numbers in parentheses indicate References at end of paper. 
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Fluorescent Technique for 
- Determining the Cross-Sectional 


“Liquid Sprays O. A. UYEHARA? 


The development of a unique laboratory instrument which measures cross-sectional drop sizes 
in dynamic liquid sprays by a fluorescent technique is described. 
ing a short duration, high intensity radiation flash into a thin sheet of radiation which intersects © 
the spray nearly perpendicularly to the spray axis. 
primary radiators due to a fluorescent dye additive, are photographed and subsequently measu red. 
The error in the measurement of the drop size was estimated from actual camera calibration data — 
through the use of a specially made reticle to be less than 10 per cent for 10-micron diameter drops, 
and was found to decrease sharply for larger diameter drops. 
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This technique consists of focus- 


The spray drops within this sheet, acting as — 


The fluorescent, optical and photo-— 
Drop size distribution data aad © 


spray photographs showing drop sizes are presented. 


Consequently many studies of the atomization and vapori-— 
zation processes have been conducted. 
numerous studies the state of our knowledge is such that it is — 
still difficult to predict rates of vaporization of sprays for a 
particular nozzle, fuel and set of operating conditions because 
of incomplete information on drop size distributions from 
nozzles and rates of vaporization in sprays. 


In spite of these 


With the feeling that a basic understanding of drop size 


distributions and rates of vaporization from liquid sprays 
was essential to the understanding of heterogeneous combus- 
tion, a research program intended to furnish information on 
these points has been in progress for some time (5-9) in the 
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‘Mechanical Engineering Department at the University of 
Wisconsin. These studies have resulted in the experimental 
confirmation of theoretical relations that predict the life his- 
tory of single drops of various liquids, subjected to an atmos- 
phere at various temperatures, pressures and velocities. The 
accuracy of the predictions, for the sizes studied, was nearly 
as great as the accuracy with which the physical properties 
are known. The application of these techniques to densely 
populated sprays (10) indicated that the assumption of a 
single drop surrounded by an infinite volume of air, unaffected 
by the vaporization of the other drops, does not agree with the 
experimental facts. In order to develop a correction factor 
for the effect of one drop on another, experimental: drop size 
distributions in vaporizing sprays at cross sections normal to 
the spray axis and located at various distances downstream 
from the nozzle are needed. 

Techniques previously used to measure drop size distribu- 
tions in sprays may be classified into six major groups: 
Collecting drops on slides or in cells and microscopic size de- 
termination, freezing of drops and mechanical sieving, direct 
photographic techniques, scattering or light absorption tech- 
niques (including Tyndall and Mie scattering, and other opti- 
cal methods), electronic or radiographic techniques, and selec- 
tive impaction. Reference (4) has a bibliography of 120 
papers written on these techniques, and discusses their ad- 
vantages and disadvantages in detail. All of these techniques 
may be criticized because the functional relation with drop 
size is either complicated, or based upon conjectural and tenu- 
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Fig. 2 Photograph of apparatus showing radiation sheet in 


ous assumptions, the accuracy of the method cannot be deter- 
mined, an undesirable physical disturbance of the spray is 
produced by the measuring technique, the physical phe- 
nomena of the measuring characteristic is inadequately under- 
stood, or the practical difficulties have seemed insurmount- 
able. In addition, rigorous direct methods of calibration un- 
der operational conditions are in most cases unavailable. 

Of the methods indicated, direct photography shows the 
most promise. However the methods used to date obtained 
cross sections of portions of the spray that were parallel rather 
than perpendicular to the spray axis. They all suffered from 
inherent uncertainties as to uncertain depths of field and im- 
age blur because of drop motion perpendicular to the cameras’ 
lens axis. 

The technique that was developed and is being reported 
herein is believed capable of determining, with good accuracy, 
cross-sectional drop sizes in vaporizing sprays. It is not 
completely free of all the mentioned criticisms, but it does 
offer the following advantages over previously used tech- 
niques: It is capable of direct and fairly accurate calibration 
for drop sizes under actual spray conditions; it produces mini- 
mum interference with the spray; it does not influence drop 
size distribution; it has the sampling plane perpendicular to 
the spray axis; it permits only a well-defined volume within 
the spray to be analyzed at one particular setting; and it pro- 
duces little or negligible drop motion-blur. 


Technique 


The apparatus, which utilizes the fluorescent characteristic 
of certain liquids or liquid dye mixtures, is shown schemati- 
cally in Fig. 1. A fluid pressurized by nitrogen is ejected 
from a nozzle, forming a conical spray which, further down- 
stream, is drawn, together with entrained air, into a duct and 
fan assembly and subsequently discharged from the test area. 
A high pressure capillary mercury-arc lamp is flashed by capac- 
itor discharge. The resulting radiation is focused into a thin 
sheet by two plano-parabolic quartz lenses having a rectan- 
gular midstop aperture. Ata predetermined location between 
the nozzle and the duct, this sheet intersects the spray nearly 
perpendicularly to the spray nozzle axis and irradiates a nar- 
row cross section of the spray. The spray drops within this 
sheet act as primary radiators due to a minute quantity of a 
fluorescent dye additive. 

This dye partially absorbs the radiation and reradiates it. 
The images of the fluorescent drops are enlarged and photo- 
graphed by a high resolution camera whose optical axis is 
perpendicular to the sheet of radiation and centered upon the 
spray area to be photographed. 

The apparatus is rigidly fixed in space with the exception of 
the nozzle, which can be positioned accurately upon machined 
ways by three micrometers in three mutually perpendicular 
directions. This construction insures nearly perfect optical 
alignment, since the spray can be completely photographed 
by simply moving the nozzle, while the optical systems remain 
fixed in space. 

Fig. 2 is a photograph of the apparatus showing the radia- 
tion sheet intersecting the spray. 

If satisfactory drop images are to be obtained, only the half 
of the spray on the camera side should be photographed. 
Reference (11) presents the details of this problem. If it is 
not desired to assume symmetry in the spray, the nozzle can 
be rotated and the other half of the spray studied. 

The brightness (power per unit solid angle per unit pro- 
jected area) that a given fluorescing drop exhibits is directly 
proportional to the product of the brightness in the sheet of 
radiation and the solid angle of radiation convergent upon the 
drop. The excitation brightness must, of course, lie within 
the absorption band of the fluorescent dye. In addition, the 
sheet of radiation illuminating the spray should be quite thin 
and have a thickness no greater than the depth of focus of 
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the camera lens. The radiation sheet should also be carefully 
positioned so that it will fall within the depth of focus. 

These objectives were met in the present apparatus by using 
a plano-parabolic lens with a rectangular midstop, plus a re- 
flecting spherical mirror, by using a source of high effective 
brightness, and by maximizing the pulse-light circuit param- 
eters so that electrical energy is supplied to the radiation 
source at a maximum rate. 


Fluorescence 


The radiation falling on the drops is absorbed and re-emitted 
because of the added fluorescent dye. The use of this 
dye causes the drop to act as a primary light source and thus 
enables better drop images to be obtained than if reflected 
light had been used. In addition, calibration of the optical 
system in place and under actual conditions is possible. 

Inasmuch as the drops are to receive radiation from a small 

olid angle and to emit over a solid angle of 47, it is necessary 
to conserve radiation in all ways possible. It is obvious that 
the fluorescent dye must efficiently convert the radiation 
rom the illuminating source to fluorescent radiation. In 
addition, the dye must absorb at the wave lengths present in 
the illuminating source and, because of the unavailability of 
suitable lenses, emit in the visible region. Also, the fluores- 
cent induction and decay times should be very short (1077 
sec or less), so that the drops will not be still fluorescing when 
they have passed out of the ‘‘depth-of-focus” of the camera 
lens. Finally, if a dye additive is used, it must have a high 
solubility in the fluids of interest. 

Fig. 3, taken primarily from data given in (12), indicates 
how well an aqueous uranin solution meets these requirements. 
It will be noted that the relative quantum yield (Q,) is very 
high throughout the whole absorption band; the absorption 
and fluorescent spectra (relative intensity 7, not to common 
scale) are strong, and the fluorescence is in the visible spec- 
trum. Further, the induction and decay times are approxi- 
mately 5 X 10~* sec, and this dye is highly soluble in glycer- 
ine, water, alcohol and alcohol-ether mixtures (which repre- 
sent a tremendous range of liquid properties). Other fluo- 
rescent dyes, such as rubrene (tetraphenylnaphthacene), are 
soluble in many fuels and exhibit characteristics that are as 
favorable to this technique as those of uranin (11). 

It should be noted that in spray fluids composed of fluo- 
rescent dyes (solutes) in solution with nonfluorescent spray 
liquids (solvents), the fluorescent brightness of the spray 
vapor is usually negligible, since the fluorescent dyes usually 
have much less volatile than the solvents; i.e., the boiling 
point of rubrene is much higher than that of benzene. There- 
fore, the spray vapor usually is not photographable. Further, 
due to the great difference in boiling points, the dye concen- 
tration in the liquid drop may increase as the drop vaporizes. 
As an example, if a drop is reduced in diameter by a factor of 
10 through vaporization, and if high internal circulation (7) 
exists within the drop, the dye concentration would increase 
roughly by a factor of 1000. In that case, if the initial con- 
centration were close to its saturation value, the liquid prop- 
erties of the drop would change as the drop vaporizes, with a 
resulting tendency for the drop to gel. In the data taken to 
date, the initial dye concentration by weight has been less 
than 1 per cent and usually about 10~? mols per liter of spray 
liquid, and the diameter change was much less than 10. Of 
course, any system improvement which further increases or 
conserves radiation would permit the use of smaller dye con- 
centrations. 

Fig. 4 shows the relative effect of dye concentration on 
brightness. C. and By are reference concentrations and 
brightness, respectively. Curves similar in shape to those 
shown in Fig. 4 but having different coordinate values are ob- 
tained for different liquids and dyes. Thus it is important to 
select the proper concentration for a given dye and liquid. 

Table 1 presents data on the effect of a 5-gm per liter uranin 
dye concentration on selected properties of different fluids. 
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Fig. 3 Dilute aqueous uranin solution spectral characteristics 
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Fig. 4 Relative fluorescent brightness vs. dye concentration 


The quality of a photograph of the drops depends upon the 
radiation from and the motion of the drop, the camera lens 
and the film. 

Speaking of the camera lens, a low f-number is desired to 
gather the maximum amount of radiation. At the same time | 
a long focal length lens is desired so that the lens does not inter- 
fere with or be wetted by the spray. Most readily available 
lenses having small f-numbers also have a short focal length. 
Also, equal manufacturing precision will result in less optical 
errors in short focal length lenses. The lens must also have 
high resolving power if small drops are to be faithfully re- 
corded. Fig. 5, taken using a high priced camera lens oper- 
ated under reasonable photographic conditions, indicates that 
the problem is not simple. 

A sensitive film is desired, since the radiation from the drop 
is very low. For the same reason, not much enlargement of 
the drop image is desired. However, if the drop image is not 
enlarged, grain size and defects in the film will become impor- 
tant. The problem in photographing the drops then is to 
maximize several conflicting requirements and characteristics 
(11). 

Since the technique demands both a true size and shape 
record of spray drops, and since a number of lenses did not 
product such images (see Fig. 5), an extensive lens testing 
program was conducted to optimize the camera optical sys- 
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tem. In order to both test and calibrate the optical system, 
an annular reticle having increasing concentric radii from 5 to 
250 u, accurate within a 1 per cent error (relative to a ruled 
stage micrometer of accuracy within 100 A) was used upon 
high resolution spectroscopic plate emulsion and a high reso- 
lution microscope objective. Such a reticle is a very critical 
test for a well-corrected optical system (11). The image of 
this reticle formed by an aerial camera lens is shown in Fig. 6. 
For comparison, the image of the same reticle produced by 
the finally adopted camera optical system is shown in Fig. 7. 


Table 1 Effect of dye additives on spray liquid properties 


Viscosity 
Measured Referenced 
viscosity® viscosity, 29 
Fluid, 20 C (c.p.) (c.p.) 

ethyl alcohol 1.2 1.200 
ethyl alcohol (dyed) 1.2 
water 1.0 1.005 
water (dyed) 1.0 
glycerine 1.5 X 108 1.490 x 108 
glycerine (dyed) 1.5 X 103 


Surface tension 


Measured 
surface Referenced 
tension,” values, 29 
Fluid, 20 C dynes/cm dynes/cm 
ethyl alcohol 23 22.3 
ethyl alcohol (dyed) 23 
water 73 72.75 
water (dyed) 7 
glycerine 64 


glycerine (dyed) 64 


Specific gravity 
Measured Referenced 
specific gravity values, 29 
100° (aleoholometer) 100° 
99° (aleoholmeter) 


Fluid, 20 C 


ethyl alcohol 
ethyl] alcohol (dyed) 


water 0° (alcoholmeter) 0° 
water (dyed) 0° (alcoholmeter) 
glycerine 29.9° (Baumé) 29.92° 


glycerine (dyed) 29.9° (Baumé) 


2 Determined by Saybolt viscometer using a Universal 
tube. 
> Determined by open-to-air capillary tube. 


Fig. 5 Spray photograph showing coma distortion of fluorescing 
drops by camera lens 


450 


It is obvious that any system which records an optical im- 
age must be precisely calibrated under actual data sampling 
conditions over the entire operating range, if the recorded data 
are to be reliable. This apparatus was calibrated by deter- 
mining the relation between the true drop size and the ap- 
parent drop size recorded on a film produced by this reticle. 
In this test the reticle was successfully positioned throughout 
the optical field volume (that spray volume determined by 
the camera’s field of view superimposed upon the thickness of 
the radiation sheet). The same type of film was used as is 
used in recording spray drop sizes. The negatives were then 
developed, processed and visually measured from their magni- 
fied images using the same procedures as used for spray photo- 
graphs. Since the size of the reticle was known and the 
camera and read-out magnifications were determined inde- 
pendently, the error could then be determined. A graph of 
apparent size vs. actual size of the diameters of the concentric 
circles in the reticle was plotted from this data and is shown in 
Fig. 8. The region between the two solid lines in the graph 
indicates the spread or indeterminancy of the experimental 
data for any given actual diameter. It may be noted that, 
since the curve as a whole is linear, diffraction effects are unde- 
tectable even at 10 uw diameters. From this graph it is evi- 
dent that, percentagewise, the size indeterminancies are 
largest for the smallest drops, and that this indeterminancy 
decreases rapidly for larger drops. The indeterminancy at 
10 uw is about 10 per cent. However, it should be cautioned 
that extrapolation of this linear curve to smaller sizes is not 
at all possible, since both aberration and diffraction effects 
would cause the curve to swing upward. Part of the inde- 
terminancy may be attributed to human error in averaging 
the edge gradient in the film image during read-out. The 
same statement, in general, may be made for drop counting 
and measuring. 

Unfortunately, fluorescing spray drops within the radiation 
sheet apparently fail to fluoresce with uniform emittance 
(emitted radiant flux per unit cross-sectional area) both across 
the cross section of any given drop and between drops. The 
nonuniformity across the cross section is attributed to both a 
fluorescent trapping effect and the variation in the optical 
thickness of a drop throughout its cross section, with conse- 
quent production of an emissivity gradient across the drop. 
The nonuniformity between drops is attributed, among other 


Fig. 6 Image of annular reticle formed by an expensive photo- 
graphic lens 
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Fig. 7 Image of annular reticle formed by camera of apparatus 


things, to the light source producing a nonuniform brightness 
across its cross section due to electron temperature gradients in 
the source’s plasma, and to the excitation optical system 
producing an imperfect image of the source, thus producing 
a nonuniform radiation field. 

This random variation in intensity of radiation from the 
drops produces images of varying intensity on the film with 
consequent “edge effects.” These edge effects (which are 
further increased by “graininess” in the film) cause some un- 
certainty in determining drop sizes from the original negative, 
and are magnified when a reproduction of the negative is 
made. Since no way of calibrating these edge effects is 
available, the original negatives are used for drop counting 
purposes. 

The foregoing discussion indicates that all data must be 
taken under the same conditions at which the apparatus was 
calibrated. If not, then not only are the spray drop size and 
size distribution data relative, but they can be varied for 
given spray conditions by simply changing the camera aper- 
ture, camera lens, film exposure, film developing and process- 
ing (including printing), or by a variable discrimination of 
the optical density value that represents the spray drop image 
boundary (edge gradient effect). 
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Data 


Fig. 9 is an unretouched photograph of an ethyl-alcohol 
spray injected at 10 psig into static ambient air, taken at a 
location on the spray axis 3 in. downstream from the nozzle. 
The other operating conditions are listed in Table 2. The 
photograph is the result of five successive exposures (flashes) 
on the same film. Multiple exposure is desirable from a 
statistical point of view. It may be noted that some drops 
are not completely within the radiation sheet and conse- 
quently appear out-of-focus, whereas the grainy appearance is 
caused by the negative emulsion material. What at first may 
be interpreted to be spray breakup, represented by the some- 
what superimposed drops in Fig. 9, is probably due to the 
multiple exposures on this photograph, since the location is 
thought sufficiently downstream to avoid such phenomena. 
Fig. 10 is a photograph having five exposures of a 100-psig 
alcohol spray taken under the same spray conditions. 

To date, spray drop photographs of similar quality have 
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Fig. 8 Calibration curve for camera of apparatus showing 
apparent reticle diameter vs. actual reticle diameter 


Fig. 9 Spray photograph of ethanol at 10-psig injection pressure 
with five flashes superimposed 


Fig. 10 Spray photograph of ethanol at 100-psig injection 
joy pressure with five flashes superimposed 
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been obtained by this technique for glycerine, water, ethyl 
alcohol, and isobuty] alcohol-diethy] ether liquids with uranin 
as the fluorescent additive. In addition, benzene, kerosene, 
diethyl ether, and petroleum lubricating oil with rubrene as 
the fluorescent additive have been successfully photographed. 
Such fluids give a wide range of viscosity, surface tension, 
vapor pressure and density. Further, all of the fluids in 
which one dye is soluble are mutually soluble and therefore 
permit a continuous alteration of the above properties, if 
desired. 

In presenting drop size distribution data it is convenient to 
use a reduced temporal drop size number represented by 
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Fig. 11 Drop size distribution for ethanol spray at 10-psig 
injection pressure 
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Fig. 12 Drop size distribution for ethanol spray - Saati’ 


ie., the number of drops (dN) in diameter range (dD) per 
total number of sampled drops (N) averaged over a time in- 
terval (t) sufficiently long to average out spray fluctuations. 
When this drop size number is plotted against drop size (di- 
ameter), a reduced temporal drop size distribution is obtained. 
In practice, these infinitesimals are finite. The drop size 
classification so obtained are represented in Fig. 11, which 
shows the percentage of drops in a particular size classifica- 
tion. The data for this bar-graph type drop size distribution 
curve were obtained for ethyl alcohol atomized at conditions 
indicated in Fig. 11 and Table 2. The size distribution for 
4000 drops, plotted here, represents over 1000 independent 
exposures—a number of exposures felt large enough to aver- 
age out any temporal spray fluctuations. 

Fig. 12 is a reduced temporal drop size distribution bar 
graph for ethyl alcohol injected at 100 psig. Fig. 13 is a 
similar graph for an equal concentration of glycerine and 
water atomized at 100 psig under the same conditions. It is 
evident that this last distribution is very similar to that of 
ethy! alcohol injected at 10 psig, as shown in Fig. 14 by curves 
AandC. These continuous curves were obtained by drawing 
a relatively smooth curve through the midpoint values of the 
bar graphs of Figs. 11, 12 and 13. It is interesting to note 
that each curve has three peaks with those of curve B being 
reversed in order of magnitude to those of curves A and C. 
Since most naturally occurring distribution curves are smooth 
curves if a sufficient number are considered, this may indicate 
that 4000 drops are an insufficient number to give a completely 
accurate distribution curve. 

In the data presented herein, all of the counted drops that 
appeared to be smaller than 10 yu in diameter were measured to 
be greater than about 5 yw, but it is believed that even in this 
5-10 uw range a high percentage was missed. The many un- 
knowns in this problem might be circumvented by not includ- 
ing any population values in the 0-5 uw range, since these values 
are, at best, in very great error and therefore could be mis- 
leading to spray atomization and vaporization studies. 

It is self-evident that unless a certain minimum number of 
drops are measured the drop size distribution obtained is a 
function of the number of drops measured. In order to esti- 


Table 2 Operating conditions for spray data 


capacitor voltage 8-10 kv 7 

capacitance 0.15 pf 

spark gap setting such as to cause self-breakdown at 
11 kv when in series with cold 
tube (around 0.10 in.) 


source BH-6, nonpre-ionized 
BH slit width 800 
radiation sheet thick- = 
ness 200 
radiation sheet, loca- 
tion 3 in. downstream from nozzle, 15 
deg off normal from spray axis 
radiation sheet area 5 X 5mm 
spray fluid tempera- 
ture (upstream from 
nozzle) 70 F 
nozzle Monarch ‘‘F-80” 30-deg solid cone, 
0.60 gal/hr at 200 psig for fuel 
oil 


camera magnification 20.00 + 0.05 
photographic film Royal-X Pan (Kodak) 


development 12 min at 68 F, DK-60a, continu- 
ous agitation 
intensification chromium intensifier followed by 


Monckhoven’s intensifier 


Visual size classification is made from measuring film 
image projected at 28.2 <X upon a white low-gloss plane 
screen (12). 
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mate if this certain minimum number were included in the 
foregoing data the following analysis was made. The first 
100 drops sampled were classified into their k predetermined 
size classification with the result that m,; drops would occur in 
the ith size classification for this first group of 100 drops, 
represented by subscript 1 (¢ ranges from 1 tok). Therefore 
the summation of m, from 1 to k will yield 100. Each n; is 
the percentage of drops in the ith size classification for this 
first group of 100 drops. This, divided by 100, is the reduced 
temporal drop size number for the ith size classification based 
upon 100 drops sampled. Let this be denoted by M,,. 
Mp, for any size classification is zero. Denote (M,; — Mo,)/ 
Mi, as percentage change in W,; for the first 100 drops sam- 
pled, which is 100 per cent. Now 100 additional drops are 
classified into the same k classifications. Consider again the 
ith classification, which has ne; drops in it. The sum of m; 
and ne; is the number of drops in the ith classification for 200 
drops sampled. This sum is divided by 200, and this per- 
centage is called M2;. Then (M2; — M,;)/Mbe,, the percent- 
ige change in W2; in the ith size classification for 100 addi- 
tional drops sampled, is calculated. 

In general then, the following procedure can be used. 
Consider a given ith drop size classification obtained for a given 
total number (VN —j) of sampled drops and the percentage 
of total drops contained therein. Count j additional drops, 
which in the case noted was 100, and calculate the new per- 
centage in this given ith drop size classification based upon all 
drops counted so far, i.e., N drops. It is evident that if any 
comparison in the rate of convergence of Wy, with N is to be 
made for various total numbers (N) of drops counted, the 
value of 7 must be fixed. Calculate the absolute value of the 
difference between these two percentages and divide by the 
new percentage. This ratio is termed percentage change in 
Wy; for this given ith drop size classification for N drops 
sampled, and serves as the ordinate for the graph shown in 
Fig. 15. As is evident, the percentage change in W x, is dif- 
ferent and converges differently for each drop size classifica- 
tion, as illustrated for the ith groups of 10, 50 and 250y for 
this aleohol spray injected at 10 psig. The low value of per- 
centage change in Wy, occurring at 1200 drops is misleading, 
since other size (ith) groups had a relatively high value at this 
point. 

It can be seen, therefore, that in determining the cutoff point 
for the number of drops to be included in the analysis, the 
envelope of least upper bound for all size classifications 
should be considered; whereas, due to fluctuations in the rate 
of convergence between size groups, it is incorrect to state that 
the ith drop size classification that approaches zero most 
slowly determines the required minimum number. It may 
be stated that the values of W y, for all drop size classifications 
were approximately within these curves given in Fig. 15 at 
cutoff; i.e., the curves shown there serve as least upper bound 
envelopes at cutoff. 
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An investigation was conducted on the flow over a flat plate with a sharp leading edge in a hyper- 
_ sonic shock tunnel. A flow Mach number range of 9.77 to 19.6 with a reflected stagnation tempera- 
ture of 2500 R, and a flow Mach number of 7.8 at 7200 R were used. For these flow conditions the 
= pressure in the test section varied from 20 microns to 1 mm of mercury, and the Knudsen 
_ number based on the leading edge radius varied from about 0.3 to 2.9. Slip flow effects at the sharp 
4a leading edge were observed showing that the strong shock wave and the boundary layer were delayed 
_ in forming on the plate. The slip region increased with the free stream flow Mach number and with 
_ the mean free path of the air ahead of the model. In the strong interaction region near the leading 
edge, the viscous region extended from the plate surface to the shock wave, and there was no appar- 
ent inviscid region between the shock wave and the boundary layer. The induced pressure on the 
= varied inversely as the square root of the distance from the leading edge as originally predicted 
The strong interaction parameter M,*/V Rey varied from 2 
_ to approximately 140, and at the highest flow Mach number tested, the maximum induced pressure 
_ near the leading edge was nearly 29 times the free stream static pressure. Cooling of the boundary 
layer by the plate surface greatly affected the magnitude of the induced pressure at high flow Mach 
- mumbers. Even for a reflected stagnation temperature of 7200 R, with appreciable real gas effects, 
_ the strong interaction between the shock wave and the bou ndary layer was present. 


control problems will be important, so that a glide-type 
vehicle looks most promising. Also, to decrease the dispersion 
of long range ballistic missiles, the nose cones can be equipped 
with fins for terminal guidance. Thus, there is a great need 
for fundamental knowledge regarding the aerodynamic 
characteristics of lifting surfaces at high flow Mach numbers 
and correspondingly high stagnation temperatures. 

For a flat plate with a sharp leading edge in a hypersonic 
flow, a shock wave and a viscous boundary layer will be 
formed at the leading edge. Recent experimental results 
(1-3)* for a flat plate at hypersonic Mach numbers have 
indicated a delay in the formation of the strong shock wave, 
where the boundary layer was merged with the shock wave for 


* Numbers in parentheses indicate References at end of paper. 
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some distance back of the leading edge before separating. 
For a given Reynolds number with insulated wall conditions, 
the boundary layer thickness is approximately proportional 
to the square of the Mach number. At hypersonic Mach 
numbers the large boundary layer growth causes a strong 
interaction between the leading edge shock wave and the 
boundary layer. Since the pressure across a shock wave is 
proportional to the square of the flow deflection angle and 
the Mach number, the induced pressures due to the inter- 
action of the shock wave and the boundary layer on a flat 
plate can be much greater than the free stream pressure at 
high Mach numbers. This increase in the surface pressure 
for a flat plate will increase the skin friction as well as the heat 
transfer. For lower supersonic flow conditions these induced 
pressure effects are small, and at low subsonic speeds the 
classical laminar boundary layer theory applies, which pre- 
dicts no pressure increase. 

The case of hypersonic flow over a flat plate with a thick 
leading edge, where the thickness is much greater than the 
mean free path of the gas in the free stream, is entirely dif- 
ferent from that of the sharp leading edge. With a blunt 
leading edge, a detached shock wave is produced ahead of the 
plate and the flow Mach number after the shock is subsonic 
with high pressures. For this condition the viscous effects 
become small and the inviscid flow over the blunt leading edge 
predominates as discussed in. (4 through 6.) Since the flow 
is subsonic ahead of the leading edge, there is no possibility 
of the existence of slip flow for a reasonably dense gas, as in 
the case of a sharp leading edge at hypersonic Mach numbers. 

A number of theoretical papers (7-13) for a flat plate with 
infinitely thin leading edge at hypersonic Mach numbers have 
been prepared in recent years. At the high flow Mach num- 
bers, the induced pressure near the leading edge caused by 
the interaction of shock wave and boundary layer will be 
large. This is usually referred to as the strong interaction 
regime. Two general approaches have been utilized in de- 
veloping the strong interaction theories. In (8 and 9) the 
boundary layer is assumed to extend from the plate surface 
to the shock wave. It was originally shown in (9) that for a 
strong interaction region the induced pressure varies as 
X~-"/:, where X is the distance from the leading edge, and the 
important parameter was M 3/ Rex, where M, is the free 
stream Mach number and Rex is the Reynolds number based 
on the distance X. The other solutions (10-12) were derived 
on the assumption that the boundary layer and the shock 
wave were distinctly separated, even in the vicinity of the 
leading edge. A theoretical result for the hypersonic viscous 
flow over a noninsulated flat plate has been derived in (18). 
Recent experimental results (1,2) for a flat plate with a sharp 
leading edge have indicated that the shock wave and the 
boundary layer are merged for some distance back of the 
leading edge before separating. The shock wave-boundary 
layer interaction induced pressure results obtained at very 
high flow Mach numbers (to be presented in this paper), 
indicate the existence of the slip flow, and the merged 
region of the shock wave, and the boundary layer for large 
values of the interaction parameter M,°/ V Rex. 

Some experiments have been conducted in air to determine 
the magnitude and nature of the shock wave-boundary layer 
interaction phenomena in hypersonic wind tunnels (14,15). 
In (14) a careful investigation of the flow between the shock 
wave and the flat plate surface has been conducted at Mach 
5.8. The interaction parameter was only of the order of 3, 
so that the induced pressure due to the interaction at the 
leading edge was quite small. By heating the air to 1360 R, 
a hypersonic flow over a flat plate was investigated at Mach 
9.6 by measuring the surface pressure (15). The highest 
interaction parameter for these investigations was about 10, 
and the maximum measured induced pressure ratio near the 
leading edge was about 9. A similar flat plate with a sharp 
leading edge was studied in the University of California Low 
Density Tunnel over Mach range 2 to 5.8 for rarefied condi- 
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tions. The flow field between the shock wave and the plate 
surface was surveyed by means of a free molecular probe (1). 
These results indicated the possible existence of the slip flow 
and the merged region between the shock wave and the 
boundary layer in the vicinity of the leading edge. The other 
method of investigating the flow over a flat plate was to 
measure the induced pressure on the surfaces of a number of 
flat plates with orifices located at different distances from the 
leading edge. Induced pressure results (3,16) have also 
indicated the existence of the slip flow at the leading edge, 
with the largest interaction parameter being about 10 and a 
corresponding induced pressure ratio of 8 at flow Mach 5.8. 
An experimental investigation of a flat plate with a sharp 
leading edge in a helium tunnel has been conducted, and the 
induced pressure results are presented in (17). For helium, 
the ratio of specific heats y is 1.67, which is greater than the 
value for air, y = 1.4. 

The present investigation was conducted in order to extend 
the available knowledge in the literature regarding the slip 
flow at the leading edge of a sharp flat plate, the shock wave- 
boundary layer interaction phenomenon at higher flow Mach 
numbers and the interaction parameters. A study was made 
to determine the existence of the leading edge slip, and the 
nature of the shock wave-boundary layer interaction at high 
stagnation temperatures where the real gas effects for air are 
appreciable. A brief description of the hypersonic shock tun- 
nel facility and the instrumentation used to obtain surface 
pressures, and Schlieren photographs of the shock wave and 
the viscous region will be presented. The leading edge slip 
flow phenomenon for a flat plate with a sharp leading edge 
was investigated by measuring the surface pressure near the 
leading edge, and obtaining Schlieren photographs of the 
delay in formation of the strong shock wave and boundary 
layer at the leading edge. Shock wave-boundary layer inter- 
action results over flow Mach range 9.77 to 19.6 and an inter- 
action parameter range of 2 to 140 will be discussed for con- 
ditions of nearly perfect gas flow. The real gas effects at high 
stagnation temperatures upon the induced pressures on a flat 
plate will be presented. : 


Facility, Instrumentation, Model and Test 


Conditions 


Hypersonic Shock Tunnel 


The investigation was conducted in the hypersonic shock 
tunnel with different throat diameters, for a conical nozzle 
with a 24-in. exit diameter, to obtain the various hypersonic 
flow Mach numbers in the test section. A detailed description 
of this facility is presented in (18). Fig. 1 shows the nozzle 
and the dump tank with windows located at the test section 
of the nozzle. Under perfect gas conditions for air (y = 1.4), 
the nominal nozzle Mach numbers of 10, 16 and 20 had throat 
diameters of 1.0, 0.33 and 0.19 in., respectively. This conical 
nozzle was attached to the end of the 4-in. ID, 103-ft long 
constant area driven tube. The shock wave was reflected at 
the entrance to the nozzle, for the present investigation, to 
increase the pressure and temperature after the reflected 
shock and to minimize the attenuation effects due to the 
boundary layer in the driven tube. 

At the entrance to the nozzle a scored aluminum diaphragm 
is used to permit the shock tube and dump tank to be evacu- 
ated to different pressures. Usually the dump tank and the 
nozzle are evacuated to very low pressures, often as low as 4 u 
of mercury, to facilitate the flow establishment in the test 
section and to minimize the strength of the starting shock 
wave. This diaphragm opens upon the arrival of the incident 
shock wave, with the bursting pressure controlled by using 
different material thicknesses and cut depths. 

Observation windows are located at the exit of the nozzle 
to permit optical investigations of the flow fields around 
models. The windows are optically selected thick plate glass. 
A sting, which holds the model in the test section, was 
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mounted separately to the floor and vibration isolated from 
the dump tank and floor by the use of Teflon bushings. This 
hollow sting was also used to carry the electrical leads from 
the pressure gages mounted in the model to the oscilloscopes. 


Instrumentation 


All the static and impact pressures were measured by 
piezoelectric pressure gages (18,19). The standard Kistler 
SLM quartz pressure pickups were used to measure the re- 
flected stagnation pressure at the entrance to the nozzle, and 
the impact pressure in the test section. These gages, large in 
size and with a small output of about 5 mv per psi, cannot be 
used to measure the very low pressures on the surface of a flat 
plate. Small barium-titanate pressure pickups have been 
developed to meet the requirements of small size and high 
output. These gages were mounted in the flat plate model 
and were dynamically calibrated by placing the model at the 
end of the constant area driven section, and using the pres- 
sure jump across a weak incident shock wave. Both the 
Kistler and the barium-titanate pressure pickups are sensitive 
to their position and type of mounting in the models. 

A modified Tektronix 535 oscilloscope (18) was used to 
determine the shock wave position as a function of time in 
the 103-ft long driven tube. Heat and ionization gages are 
used as time of arrival pickups, and the output from these 
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gages are fed into the Tektronix scope. For the normal 
operation of the shock tunnel, the shock velocity measure- 
ments are supplemented by a Berkeley counter which meas- 
ures the shock wave traversal time over the last 6} ft of the 
driven tube with 1-microsec accuracy. 

A single pass Schlieren system with parabolic mirrors is 
used to obtain the Schlieren photographs of the flow field for 
the flat plate model. The light source is provided by a suit- 
ably timed spark of 10,000-v capacitor discharge, with a 
duration of the order of 0.4 microsec. 


Model 


The steel flat plate model shown in Fig. 2 was of 5-in. width 
and 11-in. length, of which the last 5 in. was a lucite extension 
mounted flush with the top surface. This extension was used 
to minimize the trailing edge effects upon the induced pres- 
sure measurements. The bottom of the leading edge was 
beveled 15 deg. By observing the leading edge under a high 
power microscope, the leading edge thickness was estimated to 
be 0.001 in. over the central part of the plate. Thin side 
plates were attached to the flat plate, as shown in Fig. 2, to 
prevent the strong shock wave, produced by the wedge hous- 
ing the pressure gages, from disturbing the flow on the top 
surface. Since the boundary layer at hypersonic Mach 
numbers is very thick, it is necessary to use great caution 
to minimize any disturbances that may affect the boundary 
layer. 

For the measurements of the pressure induced by the inter- 
action of the shock wave and the viscous layer, there were 
five orifices in the top surface located between X = 0.193 and 
3.90 in. The shielded electrical leads from the pressure gages 
were all housed within the wedge and in the hollow sting, 
both of which were maintained at the same low pressure as 
the dump tank. Each pressure gage was carefully shock 
mounted with rubber and also insulated electrically from the 
model. A separate cylindrical housing with a hemispherical 
nose (see Fig. 2) was used to hold the Kistler gage for measur- 
ing the impact pressure at the leading edge location. All the 
gages were dynamically calibrated in the model at the end of 
the driven tube before and after each series of induced pres- 
sure measurements in the test section. The output from each 
pressure gage was amplified and then fed into a Tektronix 
535 oscilloscope equipped with a Polaroid attachment for 
recording the induced pressure. 


Test Conditions 


The flow conditions in the test section (20) at the exit of 
the conical nozzle are essentially determined by the area 
ratio for the nozzle, and the conditions behind the reflected 
shock at the downstream end of the driven tube. The equi- 
librium stagnation temperature behind the reflected shock 
wave is controlled by the strength of the incident shock 
wave, which is produced by the combustion driver, and the 
initial temperature and pressure in the driven tube. Thus 
for a hypersonic shock tunnel the stagnation temperature at 
the entrance to the nozzle can be varied over an extreme 
range, without difficulty, for any given area ratio for the 
nozzle. 

Pressures behind the reflected shock wave at the entrance 
to the nozzle are measured directly with the standard Kistler 
quartz pickups, whereas the reflected temperatures are cal- 
culated values based on the shock Mach number at the end 
of the driven tube and the equilibrium thermodynamic data 
for air (21-25). All reflected equilibrium stagnation tem- 
peratures referred to in this paper are determined in this 
manner. 

To determine the flow Mach number in the test section 
with the conical nozzle, an axial survey was conducted with 
static and impact pressure probes. For the condition of a 
reflected stagnation temperature of 2500 R at the entrance 
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to the nozzle, the slight variation in the ratio of specific heats 
was considered in the determination of the flow Mach number 
for a measured static pressure in the test section. Even at 
this relatively low reflected temperature, the slight real gas 
effects caused the flow Mach number to lower slightly, and 
increased the static temperature and the static pressure com- 
pared with the perfect gas expansion with constant y = 1.4. 
By this method the flow Mach numbers in nozzles of nominal 
Mach number 10, 16 and 20 were determined to be respec- 
tively 9.77, 15.6 and 19.6 for a reflected stagnation tempera- 
ture of approximately 2500 R. 

At an elevated reflected stagnation temperature of 7200 R 
for the nominal Mach 10 nozzle with a 1-in. throat diameter 
ind 24-in. exit diameter, the static and impact pressures were 
measured along the axis of the nozzle. Again, by assuming 
‘thermodynamic equilibrium during the expansion process in 
he nozzle, the flow conditions in the nozzle were calculated 
‘or isentropic expansion of a source flow. The measured 
static and impact pressures agreed very closely with calculated 
values for an equilibrium flow (20). By this method the 
flow Mach number in the test section was determined to be 
7.8, with a static temperature of 810 R for a reflected stagna- 
tion temperature of 7200 R with a corresponding reflected 
pressure of 500 psi. Hence, it is apparent from this result 
that the real gas effects decrease the flow Mach number in 
the test section for a given area ratio, and the static tempera- 
ture is then higher due to the energy release from the recom- 
bination of the plasma during the expansion process. A more 
detailed investigation is being conducted to determine the 
state of the hot plasma as it expands in the nozzle, for various 
reflected stagnation temperatures and pressures. a, 


Leading Edge Slip Flow for a Flat Plate es 


Mach Number, Reynolds Number, Mean Free Path and 
Knudsen Number in the Test Section 


The flow Mach numbers in the test section of the nozzle 
were determined by measuring the static and impact pres- 
sures on the axis, and correlating these pressure results for 
air (21-25) with the assumption of equilibrium flow. In 
order to minimize the real gas effects, one series of investiga- 
tions for the flat plate was conducted at a reflected stagnation 
temperature 7’; of approximately 2500 R, with nozzle area 
ratios for nominal Mach numbers of 10, 16 and 20. 

Calculation of the Reynolds number for the different test 
conditions was based on the free stream condition existing in 
the test section. For low static temperature ranges, the vis- 
cosity coefficient was obtained from results presented in (16). 
The Sutherland equation was used to determine the viscosity 
for a static temperature in the range of 810 R. These values 
of viscosity coefficients were used in calculating the mean free 
path of the gas from the equation based on the kinetic 
theory for gases (26) 


0.499 por 


where 


= density 


average molecular velocity 
= mean free path 


ll 


All of these quantities are for the ambient condition in the 
test section. To determine the effective sharpness of the 
leading edge, the Knudsen number K, = /t, based on the 
ratio of mean free path to thickness of the leading edge, was 
calculated for each test condition. 

By taking a Schlieren photograph of the flow over the flat 
plate at different delay times, it could be seen that the 
boundary layer on the plate required approximately 200 
microsee to become fully established after arrival of the 
starting shock wave at the plate leading edge. Approximately 
the same time interval was required for establishment of the 
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detached shock wave;at the impact probe, which is mounted 
on the bottom of the flat plate. For these investigations the 
flow duration after flow establishment is approximately 2.5 


millisec. The shock wave and the boundary layer require 
only a small part of the total testing time to be fully estab- 
lished. Response time of the static pressure gages in the 
model is approximately 100 microsec, and, thus, the surface 
pressure measurements should be very close to a steady flow 
condition. 


Surface Pressure Distribution Near the Leading Edge 


In Figs. 3 and 6 the ratios of surface pressure to free stream 
pressure for different Mach numbers are plotted as functions 
of X~—!/2, which is the variation originally derived in (9) for 
the strong shock wave-boundary layer interaction region on 
a flat plate with a sharp leading edge. For both of these 
figures the pressure ratio varies nearly linearly with the fore- 
mentioned parameters for X greater than the second pres- 
sure orifice located 0.712 in. from the leading edge. There is 
a sharp drop off in the increase rate of pressure ratio from the 
X~!? variation as the leading edge is approached. The first 


pressure orifice is located at X = 0.193 from the leading edge. 
This deviation from linear variation becomes greater with 
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Fig. 3 Ratio of surface pressure to free stream pressure vs. the 
inverse of the square root of distance from the leading edge, 
T; = 2500R 
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Fig. 6 Ratio of surface pressure to free stream pressure vs. the 
_ inverse of the square root of distance from the leading edge, 


increasing flow Mach number for nearly perfect gas condi 
tions, as indicated in Fig. 3. Also the asymptotic value for 
the pressure ratio near the leading edge increases approxi- 
mately as the square of the flow Mach number. At the higher 
stagnation temperature of nearly 7200 R, where el 
real gas effects are present, the induced pressure ratio near 
the leading edge also approaches an asymptotic value for a | 
flow Mach number of 7.8. (See Fig. 6.) 

In Figs. 5 and 8 the induced pressure ratios are plotted in 
terms of the hypersonic parameter M,°/ vy Rex for reflected 
stagnation temperatures of approximately 2500 and 7200 R. 
Plotted in this manner, the tapering off of the induced pres- 
sure as the leading edge is approached is quite evident, and 
the magnitude of departure from the straight line for a given 
Mach number is much greater at the highest flow Mach 
number investigated. The asymptotic value for the induced 
pressure at the leading edge increased with the flow Mach 
number. At the higher stagnation temperature of approxi- 
mately 7200 R the curvature of the induced pressure ratio, 
in terms of the interaction parameter near the leading edge, 
is well defined because of the lower free stream Reynolds 
number. 

This drop off in the increase rate of the induced pres- 
sure and the approach to an asymptotic value at the leading 
edge could be caused by the slip flow existing at the front part 
of the plate (3,16) at low density conditions for flow Mach 
numbers up to 5.8, theoretically postulated in (7). 


Schlieren Photographs of the Flow 


The Schlieren photographs presented in Figs. 9-12 of the 
flow over a flat plate were obtained with a spark light source 
of 0.4 microsee duration. Figs. 9-11 are for the flow with a 
reflected stagnation temperature of approximately 2500 R 
for flow Mach 9.77, 15.6 and 19.6. For these test conditions 
the air can be considered a nearly perfect gas with y = 1.4, 
and the real gas effects are small. In these figures the strong 
shock wave and the boundary layer are delayed in forming at 
the leading edge. The maximum shock wave angles are lo- 
cated downstream rather than at the leading edge X = 0, 
as predicted by the classical boundary layer theory with no 
slip at the surface. Also, the shock wave and boundary layer 
are merged together before becoming separated. The merged 
region increased with the flow Mach number as indicated by 
the Schlieren photographs. In hypersonic flow, the outer 
edge of the boundary layer has a large density gradient (14,27), 
so it appears in the photograph as a light region. The loca- 
tion of the reverse curvature in shock wave and boundary 
layer also moved downstream from the leading edge with 
increasing flow Mach number. 

' Because of the low density and Reynolds number at a 
reflected stagnation temperature of approximately 7200 R 
and flow Mach 7.8, the boundary layer is not distinct in the 
Schlieren photograph presented in Fig. 12. Again, the shock 
wave and baci layer are delayed in tiie and the 
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Fig. 8 Ratio of surface pressure to free stream pressure vs. the 
interaction parameter, 7’; ~ 7200 R 
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Fig. 10 Schlieren photograph of flow over a flat plate at Mach 
15.6, 7’; = 2500 R 


maximum shock angle location is downstream of the leading 
edge. This location at this temperature condition was 
greater than that of the case of a reflected temperature of 
2500 R and flow Mach 9.77. The primary cause for this 
difference in the location of the maximum shock wave angle 
(7) must be due to the larger mean free path in the test sec- 
tion. See Table 1. 

All the Schlieren photographs of the flow over the present 
flat plate with a sharp leading edge have exhibited the delay 
in formation of strong shock wave and boundary layer at the 
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Fig. 9 Schlieren photograph of flow over a flat plate at Mach 


Fig. 11 Schlieren photograph of flow over a flat plate at Mach 
19.6, 7; = 2500 R 


Fig. 12 Schlieren photograph of flow over a flat plate at Mach 
7.80, 7; = 7200 R 


leading edge, for both the nearly perfect gas condition y = 
1.4 and the higher temperature case with appreciable real gas 
effects. The shock wave and boundary layer were merged 
together before becoming separated. The location of the 
maximum shock wave angle and the merged shock wave and 
boundary layer region increased with the flow Mach number 
and the mean free path in the test section of the nozzle. The 
shock wave angle at the leading edge was always appreciably 
lower than the maximum shock angle slightly downstream. 
This delay in the strong shock wave and boundary layer 
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Table 1 Some representative test conditions for 
different flow Mach numbers in the test section 
M, R mmHg Re/m in. Kn 
9.77 2375 0.765 43,400 0.335 X 10-* 0.335 
15.6 2385 0.104 33,400 0.695 X 10-* 0.695 
19.6 2475 0.020 15,670 1.85 X 10-3 1.85 
7.80 7130 1.128 3,880 2.99 X 10-3 2.99 


* T; is the refiected equilibrium stagnation temperature at 
the entrance to the nozzle. 


formation could be caused by the existence of slip flow on the 
surface. A free molecular probe survey of the leading edge 
region, at flow Mach 5.8 and with rarefied flow conditions, 
has indicated a similar delay in the formation of a strong 
shock wave, and the existence of a merged shock wave and 
boundary layer (1). 

By measuring the shock wave angle from the Schlieren 
photograph as a function of the distance from the leading 
edge, the pressure ratio across the shock can be calculated by 
using the oblique shock wave relationship for perfect gas, as 
well as for real gas with thermodynamic equilibrium (28). 
From the leading edge the pressure ratio across the shock 
wave increased with X up to a peak value corresponding to 
the maximum shock wave angle. After this maximum value 
the pressure ratio continuously decreased as the shock wave 
angle became smaller. For flow Mach 19.6 the pressure ratio 
near the leading edge was approximately one-half the maxi- 
mum pressure ratio value, slightly farther downstream from 
the leading edge. The result for the stagnation temperature 
of approximately 7200 R and flow Mach 7.8 has indicated a 
much lower pressure ratio at the leading edge compared to the 
peak pressure ratio. 

This increase in the pressure ratio behind the shock wave to 
a Maximum pressure ratio could be caused by the slip flow 
existing in the leading edge region (2,3,16). As the slip 
velocity at the surface is brought to rest from the free stream 
value in the slip distance ¢ (7), the pressure ratio must 
increase from a finite value to a maximum value at the end of 
the slip distance. Beyond this value of X the velocity at the 
surface will be equal to zero, which is the case for the conven- 
tional boundary layer in a continuum flow. ‘ 


Slip Flow Theory 


In the vicinity of the leading edge of the flat plate, the 
surface pressures and the Schlieren photographs at various 
flow Mach numbers and stagnation temperatures have indi- 
cated the possible existence of the slip flow region, as discussed 
previously. A first-order slip flow theory has been developed 
(7) to estimate the slip distance, the delay in the formation of 
the strong shock wave and boundary layer, and the possible 
variation of the slip flow region. The slip distance based 
upon the kinetic theory of gases is 

¢~NMA 

where WN is the number of collisions required to bring the mole- 
cules to rest on the surface. It has been found experimentally 
that the value of N is approximately five collisions; M, is the 
free stream Mach number, and A is the mean free path in the 
test section. For a Reynolds number based on the free stream 
conditions and the slip distance ¢, the no slip regime is defined 
as 


Slip flow exists when 
as indicated in Fig. 13. 
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In the slip distance ¢ where the velocity at the surface is 
brought to rest, the pressure ratio must increase from a 
finite value at the leading edge to a maximum value at X = 
¢ where the velocity of the gas at the wall is equal to zero. 
Also, in the slip region the strength of the shock wave must 
increase from a finite value at the leading edge to a maximum 


shock angle at the end of the slip region. The shock wave 
and the boundary layer are both merged and delayed in form- 
ing in the slip region. The present experimental results for a 
flat plate with a sharp leading edge have indicated the exist- 
ence of the slip flow region at high flow Mach numbers and 
reasonably large mean free path. Rarefied flow results for a 
flat plate up to flow Mach 5.8 have also indicated the existence 
of slip flow (1,16). 


Shock Wave-Boundary Layer Interaction 


Induced Pressure Distribution for Reflected Stagnation 
Temperature of Approximately 2500 R 


The actual measured static pressure ratios, surface pressure 
to the free stream pressure, are plotted as a function of 
1/VX, X and M,3/+/ Rex, in Figs. 3-5, respectively, for a 
reflected stagnation temperature of 2500 R. For these 
investigations the flat plate surface temperature was close to 
room temperature throughout the test duration. By assum- 
ing a laminar boundary layer and a recovery factor of 0.85, 
the ratio of the wall temperature to the recovery temperature 
is approximately 0.25 for Mach range 9.77 to 19.6 at the 
previously noted reflected stagnation temperature. Thus, the 
cooling effects of the plate surface upon the induced pressures 
caused by the shock wave boundary layer must be considered. 
A theoretical result for the effects of a noninsulated flat 
plate, for both cooling and heating of the boundary layer, is 
presented in (13). The real gas effects upon the induced 
pressure are small for this particular stagnation temperature 
compared to the cooling effects. 

It is evident from Figs. 3-5 that the maximum pressure 
ratio near the leading edge increases approximately as the 
square of the free stream Mach number. For the present 
flat plate, the first static pressure orifice is located at 0.193 in. 
from the leading edge. At this location the induced pressure 
ratios are close the asymptotic value and have deviated 
appreciably from the X~—/? variation which seems to apply in 
the strong interaction region. This similar variation of the 
induced pressure at flow Mach range 3.5 to 6 has been ob- 
served in the low density tunnel (3,16). For all the Mach 
numbers that were investigated, the induced pressure ratio 
from the second static pressure orifice, located at 0.712 in. 
back from the leading edge, does vary inversely as the square 
root of the distance (see Figs. 3 and 5) as originally predicted 
in (9). The asymptotic values of the induced pressure ratio, 


as X becomes large, approach unity for all the Mach numbers 
investigated. These results also indicate that the lucite 
extension at the end of the flat plate and the fins at the sides 
of the plate, as in Fig. 2, were effective in reducing the ex- 


ternal disturbances on the boundary layer growth. 
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Fig. 13 Hypersonic flow near the leading edge of a flat plate 
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The induced pressure ratios have been plotted as a function 
of the distance X from the leading edge for the different Mach 
numbers in Fig. 4. For flow Mach 9.77 the peak pressure 
ratio near the leading edge is nearly 6, and at 4 in. back from 
the leading edge the pressure ratio is approximately 2.5. 
Thus, at even this relatively low hypersonic Mach number, 
the effect of the shock wave-boundary layer interaction is 
appreciable. For flow Mach 19.6 the maximum pressure 
ratio near the leading edge is nearly 29, and for X = 4.0 in. 
the induced pressure is 10 times the free stream value ahead 
of the flat plate. The skin friction and heat transfer (13) for 
the plate will be much higher for this flow condition compared 
with values determined by the compressible boundary layer 
without pressure gradient. The results for flow Mach 15.6 
lie between the induced pressure results for flow Mach 9.77 
and 19.6. 

In Fig. 5 the induced pressure ratio for the different Mach 
numbers are plotted as functions of the interaction parameter 
M,3/ V Rex, where M, is the free stream Mach number and 
Rex is the Reynolds number based on the free stream condi- 
tion and the distance X from the leading edge. For each flow 
Mach number the induced pressure ratio approaches an 
asymptotic value of different magnitude, indicating the 
existence of the slip phenomenon near the leading edge for 
the present plate with a sharp leading edge. After the leading 
edge slip region, the induced pressure ratios for all the Mach 
numbers vary linearly with the strong interaction parameter. 
The slope of the straight portion of the induced pressure ratio 
in Fig. 5 decreases as the Mach number increases. The 
theoretical slopes predicted for the insulated case (9,10) are 
plotted in Fig. 5. At flow Mach 9.77 the analytical slopes for 
insulated conditions are reasonably close to the observed 
experimental results for noninsulated conditions. The results 
obtained in the low density tunnel (3,16) for induced pressures 
on a flat plate agree reasonably well with the present Mach 
9.77 results. 

In (13) the effects of noninsulated conditions on the induced 
pressure, skin friction and heat transfer have been determined 
analytically for a sharp flat plate with perfect gas assumptions. 
For a stagnation temperature of 2500 R, the gas can be con- 
sidered nearly perfect without an appreciable error. The 
theoretical slope for the induced pressure ratio as a function 
of the interaction parameter is plotted in Fig. 5 for the actual 
temperature ratio for the plate surface. It is interesting to 
note that the result for the highest flow Mach number, 19.6, 
has the closest agreement with the analytical prediction for a 
noninsulated surface. At high flow Mach number, the sur- 
face cooling of the boundary layer decreases the induced pres- 
sures produced by the shock wave-boundary layer interaction. 
But at the lower flow Mach 9.77, the cooling effects of the 
plate surface does not seem to be appreciable. For flow Mach 
15.6 the slope of the induced pressure ratio curve in the strong 
interaction region is approximately in the middle of the slopes 
for Mach 9.77 and 19.6, as indicated in Fig. 5. Hence for 
this flow Mach number of 15.6 the surface cooling effects are 
noticeable. Also, the induced pressures are decreased in com- 
parison with the insulated plate values, but the slope is 
appreciably higher than the theoretical value for the non- 
insulated flat plate. 

Real Gas Effects on Shock Wave-Boundary Layer 7 oe 
Interaction 


For a reflected stagnation temperature of 7200 R and pres- 
sure of 500 psi, the flow Mach number in the test section was 
7.8 with a static temperature of 810 R. Behind a normal 
shock wave in the test section, for this condition, the density 
ratio is increased approximately 50 per cent over the perfect 
gas case y = 1.4, due to the real gas effects (20). By assum- 
ing the boundary layer to be laminar and with a recovery 
factor of 0.85, the ratio of the wall to the recovery temperature 
is 0.0872. The boundary layer is relatively cool because of 
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the high stagnation temperature for the free stream. In 
Table 1 the other physical parameters in the test section for 
this temperature are presented. Because of the high stagna- 
tion temperature, the free stream Reynolds number per inch is 
quite small compared with the results for the lower tempera- 
ture. 

Even at this relatively high temperature and with, com- 
paratively, very cooled surface, the induced pressure caused 
by the interaction of the leading edge shock wave with the 
boundary layer is present, as indicated in Figs. 6-8. The 
pressure ratio near the leading edge again approached an 
asymptotic value, and in the strong interaction region behind 
the leading edge slip phenomenon, the induced pressure ratio 
varied as a function of X~!/?, as indicated in Figs. 6 and 8. 
The observed pressure ratio varied from approximately 4.4 to 
2 over the front 4 in. of the flat plate. Since the Reynolds 
number per inch is relatively low compared with its value at 
lower stagnation temperatures, the pressure over the surface 
did not vary as much as for flow Mach 9.77 (see Fig. 4). In 
Fig. 8 the induced pressure result has been plotted as a func- 
tion of the interaction parameter M,3/~/ Rex. Near the 
leading edge the surface pressure was obtained for the value 
of the interaction parameter of 17.3. The leading edge slip 
phenomenon existed down to values for the interaction 
parameter of 10. Below this value the induced pressure ratio 
varied linearly with the, strong interaction parameter. The 
theoretical value of the slope of the induced pressure ratio 
has been calculated by using the results of (13) and assuming, 
as an approximation, the ratio of specific heats equal to 1.4. 
In the free stream this assumption is reasonable because the 
static temperature is only 810 R. For this test condition the 
theoretical value of the induced pressure is less than the ex- 
perimental results as shown in Fig. 8. In comparison, the 
theoretical insulated plate slope (9) has been plotted in that 
figure. The low ratio of wall to recovery temperature has a 
great effect on the induced pressure, even with an appreciable 
amount of real gas effects present in the boundary layer. 


Pressure Ratio Across the Shock Wave 


The pressure across the shock wave was determined from 
the Schlieren photographs, Figs. 9-12, of the flow over the 
flat plate at various flow Mach numbers and stagnation tem- 
peratures. In these photographs the shock wave can be 
observed as a sharp discontinuity in the density gradient, and 
the outer edge of the boundary layer appears as a light region 
where the density gradient is large for a hypersonic laminar 
boundary layer. The shock angle was measured from the 
Schlieren photographs. The free stream Mach number and 
the ratio of specific heats y = 1.4 were used to calculate the 
pressure ratio across the shock wave from an oblique shock 
relationship. Since the free stream temperature in the test 
section is 810 R for a reflected stagnation temperature of 
7200 R, the ratio of specific heats of 1.4 is valid as long as the 
shock wave angle is not large (28). 

All the Schlieren results have indicated that the pressure 
ratio at the leading edge based upon the shock wave angle is 
appreciably lower than the pressure ratio calculated for the 
maximum shock wave angle located slightly downstream of 
the leading edge. After the maximum shock wave angle, the 
shock wave becomes weaker as the distance from the leading 
edge increases. The shock wave and the boundary layer seem 
to be merged before separating into a distinct shock wave and 
boundary layer. As the flow Mach number was increased, 
the merged region extended over a larger portion of the plate, 
as in Figs. 9-12. Thus, the pressure ratio increased to a 
maximum value and then decreased in a manner similar to 
the surface pressure ratio. Similar results were observed for 
the pressure ratio determined from the shock wave angle 
obtained from Fig. 12, where the stagnation temperature is 
7200 R. Again, the calculated pressure ratio at the leading 
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edge was less than slightly downstream where the maximum 
shock wave angle was located. 

In the region where the shock wave and the boundary layer 
were merged, the pressure calculated from the shock wave 
angle and the surface pressure was approximately the same. 
This indicates that the pressure in the strong shock wave- 
boundary layer interaction region is nearly constant from the 
surface to the shock wave. This assumption was utilized in 
the analysis of the strong interaction between the shock wave 
and the boundary layer in (9). When the shock wave became 
separated, the pressure ratio across the shock wave was 
greater than that measured at the surface for a given distance 
from the leading edge, as previously observed at lower Mach 


numbers (14). 


The flow over a flat plate with a sharp leading edge has been 
investigated in a hypersonic shock tunnel with air at flow 
Mach 9.77, 15.6 and 19.6 for nearly perfect gas conditions. 
At flow Mach 7.8 the reflected stagnation temperature of 
7200 R was utilized to determine the real gas effects on the 
shock wave-boundary layer interaction. Schlieren photo- 
graphs and surface pressure distributions for the flat plate 
were obtained for each test condition, and the interaction 


parameter M,°/+~/ Rex has been extended to 140. 

Slip flow phenomenon was observed at the leading edge of 
the flat plate. The strong shock wave and the boundary 
layer were delayed in forming, and the surface pressure near 
the front deviated from the strong interaction variation. 
The important parameters governing the slip distance are the 
free stream Mach number, mean free path and the number of 
collisions required to slow the velocity at the wall to zero. 

In the region of the strong interaction between the shock 
wave and the boundary layer, the induced pressure on the 
flat plate for nearly perfect gas condition y = 1.4 varies 
inversely as the square root of the distance from the leading 
edge, as originally predicted by Li and Nagamatsu. At 
higher flow Mach numbers the surface cooling effects upon 
the induced pressure become appreciable and the results 
agree reasonably well with the noninsulated theory for the 
strong interaction phenomenon. For flow Mach 9.77 the 
cooling effects were not large and the results did not agree 
with the calculated value. The maximum induced pressure 
ratio near the leading edge increased approximately as the 
square of the Mach number for the present test conditions. 

Even for a reflected stagnation temperature of 7200 R, 
with appreciable real gas effects, the strong shock wave- 
boundary layer interaction was observed. The variation of 
the induced pressure on the plate was predicted reasonably 
well with the noninsulated theory, using the ratio of the wall 
temperature to the recovery temperature. 

The pressure ratios determined from the shock wave angle 
and the free stream Mach number have indicated the in- 
crease in pressure ratio from the leading edge to a maximum 
value slightly downstream. Beyond the location of the maxi- 
mum value, the pressure ratio decreased similarly to the 
surface pressure variation. When the shock wave and the 
boundary layer were merged, the pressure from the shock 
wave to the surface for a given location from the leading 
edge was approximately constant. After the shock wave 
became separated from the boundary layer, the pressure 
behind the shock wave was greater than the surface pressure 
at the same location from the leading edge. 

Higher induced pressures at high flow Mach numbers, 
caused by the interaction of the shock wave and the boundary 
layer, will increase the skin friction, heat transfer and force 
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Conclusions 


on a lifting surface, as compared with the values calculated by 
the compressible laminar boundary layer theory with no 
pressure gradient. 
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Numerical Between Exact 
and Approximate Theories of Hyper- 
sonic Inviscid Flow Past Slender Blunt 


numerical calculations. 


HE INTEREST in hypersonic flow around blunt bodies 

has been stimulated in the last few years by the inter- 
continental ballistic missile and satellite programs. Meteor 
physicists have been interested for many years in the trails 
left by meteors (1). Recently, interest in the wakes left in 
the atmosphere by a body moving at hypersonic speed has 
also been stimulated by the desire to monitor objects entering 
from space. Among effects that can be observed during at- 
mospheric entry are the radiation emitted by the hot gas left 
in the trail and the reflection of microwaves by the ionized 
material. In order to design adequate detectors of these 
effects, an estimate is required of the flow field surrounding a 
blunt object. 

Although a survey of the present status of inviscid hyper- 
sonic flow theory has been made by Hayes and Probstein (2), 
it should be mentioned that numerical calculations of the 
entire flow field about slender blunt nosed bodies at hypersonic 
speeds in equilibrium air have never been reported in the 
literature. Ferri (3) and Casaccio (4) have given results for 
the pressure distribution on a hemisphere-cylinder at M., = 
20 in a perfect gas to a position on the cylinder about 15 
nose-radii away from the nose. Numerical calculations are 
deemed necessary for evaluating the merits of approximate 
theories suggested in the literature for determining the flow 
field away from the blunt nose. 

For the purpose of carrying out this evaluation, we will con- 
sider the features of the equilibrium flow about a hemisphere- 
cylinder about 75-radii long flying through the atmosphere at 
60,000-ft altitude at 17,500 fps. The reason for choosing this 
particular slender body is that it is plausible to assume that 
the features of the gas flow far away from the nose of the body 
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Nosed Bodies 


fied on rational or theoretical grounds. 
of energy is radically different than that given by blast wave theory. 
those considered here are of interest, the only reliable approach at the present time is to carry out 
It was found that for certain flight velocities the pressure on the body 
does not decay to free stream pressure monotonically but overexpands. 


SAUL FELDMAN? 


Avco-Everett Research Laboratory 
Everett, Mass. q 


"This paper presents numerical results of exact calculations of the inviscid equilibrium flow about 
a long hemisphere-cylinder in motion at hypersonic velocity. 
wave as well as free layer theories of hypersonic flow. 
that the second-order blast wave theory can be used for the purpose of finding the shock shape and 
the body pressure distribution. However, this procedure is definitely empirical and cannot be justi- 


A comparison is made with blast 
As a result of the comparison, it is concluded 


We show that the presently calculated radial distribution 
If body shapes other than 


chosen here will not be too different from the inviscid flow 
surrounding the viscous wake of a short blunt body (since the 
viscous wake core could probably be assumed to be replaced 
by a solid body without introducing very large errors). At 
sufficiently high altitudes the flow will not be in thermody- 
namic equilibrium. However, this effect will not be considered 


here. 


Theories for the overall features of the flow have been pro- 
posed by Lin (5), Lees and Kubota (6), Lighthill (7), and 
Hayes and Probstein [(2), pp. 81-92]. Earlier Russian work 
has been summarized by Sedov (8). These researchers at- 
tempt to determine the general features of the flow field by 
using constant energy solutions or Newtonian theory. Lin (5) 
solved the problem of a strong cylindrical blast wave for a 
perfect gas of constant specific heat and used his results, 
together with Hayes’ equivalence principle, for determining 
the shape of the strong portion of the shock in front of a body 
moving steadily at hypersonic velocity. The basic assump- 
tion was that the gas energy, in a plane normal to the body 
axis, per unit length in the axial direction, is a constant. 
This was equated to the product of body drag force and flight 
velocity. The blast wave theory is valid for a body of zero 
thickness and is of asymptotic nature, its results being invalid 
locally near the nose. 

The shock shape for an axial symmetric body obtained from 
the “blast wave” theory is a parabola. Sakurai (9,10) also 
solved the problem of the axisymmetric blast wave and, in 
addition, analyzed the two-dimensional blast wave case. He 
used values for the specific heat ratio y of 1.2, 1.4 and 1.67. 
He also solved higher order approximations that account for 
the free stream pressure being finite instead of zero. A sum- 
mary of the results for different specific heats y is given in 
Table 1, where r is the radial ordinate of the shock, R denotes 
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Table 1 Summary of first-order blast wave results for 
axisymmetric flow 
7 Paxis/Po Pehock/Paxis r/R 
1.20 0.0925M.*/(x/R) 3.32 0.920 (2/R)'/2 
1.40 0.133 M 2.68 1.104 (x/R)'/ 
1.667 0.177 M..*/(x/R) 2.19 1.273 (2/R)'/ 


the nose radius, M., the free stream Mach number and z the 
axial distance measured from the stagnation point. The re- 
sults presented are: 

1 The ratio of the pressure along the axis of symmetry of 
the flow paxis to the free stream pressure po. 

2 The pressure ratio Pehock/Paxis, WheTE Pshock denotes the 
pressure behind the shock at the axial position given by the 
z/R coordinate. 

3 The shock ordinate r nondimensionalized with respect 
to the nose radius R. 

For the purpose of evaluating the merits of the blast wave 
theory, we will now make a comparison between some of its 
results with some results from a characteristic calculation. 
We will compare the ratio of pressure behind the shock to the 
pressure at the axis of symmetry (as obtained from blast wave 
theory) and a characteristic calculation. Ferriand Pallone (11) 
calculated the flow around a double wedge with a leading wedge 
half-angle of 45 deg and trailing wedge of 12 deg in a free 
stream Mach number of 20. They used a specific heat ratio 
y of 1.4 and carried out the calculation to about 35 body 
thicknesses downstream of the leading edge. A comparison 
of the pressure ratio mentioned previously, as a function of 
position, with the pressure ratio of 2.56 that one can obtain 
from the two-dimensional blast wave results of Sakurai (9,10) 
(with y = 1.4), shows that the agreement improves when 
away from the nose region (Fig. 1). From blast wave theory 
one would expect worse agreement, since the shock is not 
strong there. Thus, blast wave theory is unsatisfactory in 
this case. Lees and Kubota (6) used the work of Sakurai (9,10) 
in order to improve Lin’s results by accounting for the fact 
that the free stream pressure is finite and not zero. Their re- 
sults for y = 1.4 are 


laxi 0.133 M.? 


2\'/2 2/R 
= 1.104 (7) (: + 0.81 


Lighthill (7) and Hayes and Probstein [(2), p. 85] determine 
the pressure distribution on the blunt portion of a body by 
assuming Newtonian flow with a centrifugal correction. For 
a hemisphere, this assumption leads to vanishing pressure 
at a point on the surface whose normal makes an angle of 60 
deg with respect to the flight direction. Within the New- 
tonian framework the shock is, of course, parallel to the body 
up to the same point. In order to determine the shock shape 
beyond 60 deg, these authors make the assumption that the 
pressure downstream is zero. The physical meaning of the 
assumption is that the mass of the gas that goes through the 
shock moves parallel to it after crossing it, and that the cen- 
trifugal force on the shock layer permits the pressure to drop 
from the value that it had behind the shock to a value of zero 
behind the shock layer. It should be noted that the free layer 
theory is strictly valid in the limit obtained when yy approaches 
unity. 
using this so-called ‘free shock layer” theory is cubic. This 
should be compared with the parabolic shock shape obtained 
from first-order blast theory. 

The results of numerical calculations about a hemisphere- 
cylinder flying through the atmosphere at an altitude of 60,000 
ft and at a velocity of 17,500 fps will be discussed next. 
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The shape of the shock layer that one obtains when _ 


Results of Numerical Calculation for the Flow 
Field and Comparison With Approximate Theories 


In order to calculate the supersonic flow field around the 
body it was necessary to first determine the subsonic flow 
field. This was done using a method first suggested by 
Dorodnitsyn and then applied to the blunt body problem by 
Belotserkovskii [(2), pp. 214-226]. This numerical calcula- 
tion was carried out for a hemisphere to first approximation 
only, for an ideal gas with a ratio of specific heats y of 1.2; 
this value of y was chosen as being reasonable for a simplified 
calculation representing the behavior of air at the high flight 
velocity of interest here. Although this calculation was not 
very accurate, that is probably of no consequence here. The 
output of the calculation downstream of the sonic “‘line” and 
limiting characteristic was used as the input to a characteristic 
calculation (12) using the properties of dissociated air in 
thermodynamic equilibrium. The important features of the 
flow field are given in Figs. 2 through 5. 


Streamlines (Fig. 2) 


One interesting observation that can be made from the 
streamline pattern is that the gas that crosses the shock be- 
tween the axis of symmetry and a point on the shock inter- 
sected by the perpendicular to the body at the sonic point, 
occupies about half the shock layer thickness when the gas 
expands around the body to free stream condition. The rest 
of the shock layer is occupied by the gas that goes through 
the remaining portion of the shock. It is interesting to note 
that the gas close to the body at the sonic point expands so 
drastically as it flows around the body, that, for example, a 
streamline at a distance of 0.012 nose-radii away from the 
surface ends up 2 nose-radii away from the body when far 
away from the nose region where the pressure is near free 
stream. 

The calculated flow field has left, as may be expected, only 
vestiges of the results of the “free layer’’ theory of Lighthill (7) 
and Hayes and Probstein [(2), p. 85]. 


Density (Fig. 3) 


It should be noted that near the hemisphere-cylinder tan- 
gency point (Fig. 3 top), the density across the shock layer 
increases by a factor of about 33 or 4 when going from the 
body to the shock in a direction perpendicular to the axis. 
However, at distances far from the nose region (Fig. 3 bottom) 
this factor becomes 20 or 30. 
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Fig. 3 Density, temperature and pressure profiles. Top: In the neighborhood of the nose. Bottom: Far from the nose 
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The temperature is nearly constant across the shock layer 
in the neighborhood of the nose region (Fig. 3 top). Far away 
from the nose (Fig. 3 bottom) the temperature is nearly con- 
stant near the shock for about one half the shock layer thick- 
ness and then rises to a high value as the body is approached. 
The temperature along the cylindrical portion decays at a 
rather slow rate. It should be noted that the high temperature 
region occurs where the density is low. Fig. 3 (top) also 
shows the entropy profile across the shock layer at a point not 
far from the hemispherical nose. It should be noted that the 
entropy gradient has a very slow variation across the shock 
layer; the entropy changes across the shock layer only by 
about 30 per cent. Reference is often made to the large en- 
tropy gradients that exist in connection with curved shocks at 
hypersonic speeds. We here conclude from Fig. 3 (top) that a 
far more physically significant thermodynamic variable than 
the entropy gradient is the density gradient, since large 
density gradients can exist when a gas goes through a curved 
shock at hypersonic velocity. a 


Temperature (Fig. 3) 


Pressure (Fig. 3) 


The pressure increases across the shock layer when going 
from the body to the shock by a factor whose value becomes a 
maximum of about 6 when near the hemisphere-cylinder 
tangency region; this factor changes to about 2 when 55 nose- 
radii away, and from there on it decreases further. 

The velocity profiles (Fig. 4) change curvature from con- 
cave to convex at about 3 nose-radii away from stagnation 
point. Beyond 10 nose-radii the velocity at the body is about 
80 per cent of free stream, and in about half a shock layer it in- 
creases to just about the free stream value. This means that 
a body flying through the atmosphere at a velocity of 17,500 
fps would leave a jet of gas whose velocity could change from 
about 20 per cent of the flight velocity when near the axis to 
zero velocity when the gas is 5 or 10 body-radii away from the 
jet axis. Of course, the velocity on the axis, where viscous 


Velocity (Fig. 4) 


effects are important, would be larger than 20 per cent of free 
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stream velocity; thus, due to mixing, the velocity of the in- 
viscid wake mentioned previously could be increased. 


Body Pressure Distribution 


Fig. 5 shows the pressure distribution along the body as 
determined by the characteristics calculation. Free stream 
pressure is reached at about 80 nose-radii from the stagnation 
point. From there downstream, the flow overexpands. The 
manner in which the pressure recovers to the free stream value 
was not computed. For the purpose of comparison, the pres- 
sure distribution obtained by using first- and second-order 
blast wave theory is also given. In the region away from the 
nose the second-order blast wave theory with y = 1.4 seems 
to disagree at most by 25 per cent. From Fig. 5 the limitation 
of the constant y assumption in the blast wave theory, when 
applied to the present problem, is quite clear, since in the 
actual flow y varies approximately between about 1.2 and 1.4. 


T 7 7 4 
Alt.= 60,000 ft | 

Vey * 17,500 ft/sec 


! | 
CHARACTERISTICS | 
~ 
-SECOND ORDER BLAST 
Pst. 
ion 
FIRST ORDER BLAST~ 
WAVE 


° 10 20 30 40 50 #60 70 #480 #90 100 
x 
R 


Fig. 5 A comparison of body pressure distribution with blast 
wave results 
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Shock Shape 
A comparison was made of the results of the characteristic 
calculation and a least square fit for the shock shape down- 
stream of the sonic region, using a quadratic and a cubic equa- 
tion; the cubic equation gave a shape which was indistinguish- 
able from the characteristics results. Fig. 6 shows the shock 
shapes obtained from first- and second-order blast wave 
theory with y = 1.4, and from the free layer theory of Light- 
hill (7) and Hayes and Probstein [(2), p. 85]. When using 
the free layer theory, an infinite density ratio was assumed 
across the bow wave. Thus, the shock is parallel to the body 
up to 60 deg where the theory yields a vanishing value for the 
pressure. The free layer theory yields a shock shape which 
disagrees radically with the characteristic calculation. The 
principal reason for the large disagreement is the assumption 
of zero pressure behind the shock layer beyond 60 deg. How- 
ever, the free layer shape could be improved further by using 
a finite instead of zero pressure behind the shock layer. 
Although the general shock shape obtained from the second- 
order blast wave is nearly parallel to the numerically calcu- 
lated shape, there is no valid reason for using a value of 1.4 for 
y except when about 20 radii away from the stagnation point, 
where the gas temperature is sufficiently low. In order to 
coincide with the calculated shock shape, the blast wave 
shock would have to be displaced outward by about 0.4 R. 
One reason for the inadequacy of blast wave theory is that 
it assumes that the energy per unit length, in a plane normal 
to the body axis, is a constant. Fig. 7 shows this energy dis- 
tribution as obtained from the characteristics calculation non- 
dimensionalized with respect to the blast wave first-order 
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of the thermal and radial kinetic energy (based on velocities 
relative to a coordinate system fixed on the body) per unit 
length within the shock envelope; and curve (b) corrects for 
the fact that the free stream energy is not negligible when 
away from the nose. The energy level is not only incorrect 
but also is not constant. (A discussion of the fallacy involved 
when assuming in blast wave theory that the energy per unit 
length is equal to the drag will be given in a future paper.) 
However, since we radius of the shock wave and the pressure 
depend on the } and 4 powers, respectively, of the energy, 
the errors in these quantities due to evaluating the energy 
incorrectly are partly concealed. 

One would conclude that blast wave theory might be used 
more as an empirical result rather than a justifiable theory. 


Concluding Remarks 
oncluding Remarks 


In view of the fact that the density is low in the neighbor- 
hood of the body as compared to regions closer to the shock, 
it seems reasonable to assume that a certain amount of change 
in the body shape will not affect the general features of the 
flow field about the hemisphere-cylinder. In order to check 
this assumption, the shock shape around a hemisphere-cone- 
cylinder (a fairly short cone) is presented in Fig. 8 together 
with the shock wave shape of Fig. 1. It can be seen that the 
shock shape is not affected much by the change in nose shape 
and cylinder radius. 

As a result of the comparison made between the results that 
different theories yield for the flow about a long hemisphere- 
cylinder, we conclude that the second-order blast wave theory 
could be used for the purpose of finding the shock shape and 
the body pressure distribution. With this information and 
conservation of mass consideration [(2), pp. 210-214], one 
could find approximate values of the remaining quantities of 
interest in the flow field. This procedure is definitely em- 
pirical and cannot be justified theoretically. 
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_ at the inlet of the plenum chamber. 
tinea rized viscous compressible flow equations in the entire subsonic diffusion regimes with assump- 
= of one-dimensional quasi-steady flow and linear velocity distribution. 
_self-excited shock oscillation is not possible in absence of viscous dissipation. 
= to be due to further entropy generated by viscous dissipation in addition to shock oscillation. 
Neutral stability boundaries are obtained in terms of the ratio of perturbation entropy to the corre- 


sponding shock velocity and are in agreement with Stoolman’s experimental data. 
analysis confirms Trimpi’s theory of instability and also justifies the Sterbentz-Evvard theory in 


If body shapes other than those considered here are of in- 
terest, the only reliable approach at the present time is to 
carry out numerical calculations. These are not too difficult 
using the presently available high speed computers. 
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Literature Review | 


HE PROBLEM of instability of supersonic inlet diffusers oe 


(commonly known as ramjet buzz) was first observed by | 
Oswatitsch (1)4 in 1944, and has received extensive attention 
during the recent years. Experimental research on this sub- 
ject has been carried out by many investigators for the pur- 
pose of explaining the causes and nature of the instability. 
Ferri and Nucci (2) found that the vortex sheet generated at 
the three-shock intersection impinging on the cowling lip 
caused the choking of the entering flow to initiate buzz. 
Trimpi (3) observed that the flow separation on the central 
body is another factor contributing to the instability. Trimpi 
(4) proposed that the stability of flow is dependent upon the 
instantaneous ratio of the increase of entropy to the decrease 
o! mass flow, and that it is not necessary in general to block 
or even to choke the inlet. Dailey (5) still emphasized that 
the instability is caused by a strong interaction of the sub- 
critical shock with the boundary layer on the central body, | 
and that this interaction blocks or chokes the inlet. Con- 
troversy in opinions between Trimpi (6) and Dailey (7) indi- 
cites that further investigation is needed. 

Theoretical analyses have also been carried out extensively 
concerning the buzz frequency, amplitude and stability limit. 
Sterbentz, Evvard and Davids (8,9) applied the Helmholtz 
resonator approximation to determine the buzz frequency and 
amplitude. These authors concluded that the instability 
would occur if the curve of diffuser pressure recovery vs. cap- 
tired mass flow has a positive slope of sufficient magnitude. 
Trimpi (4) improved that analysis by using quasi one-dimen- 
sional unsteady flow, and proposed his forementioned theory 
for instability. In fact, the instantaneous entropy generated 
should be contributed by both viscous dissipation and shock 
oscillation. Stoolman (10) first raised the interesting ques- 
tion of whether a self-excited shock oscillation can exist in the 
absence of viscous dissipation. Stoolman [and later Mirels 
(11)] formulated a stability criterion in terms of the acoustic 
impedance at the inlet of the plenum chamber. The acoustic 
impedance of the external diffusion regime for a normal shock 
inlet diffuser in the quasi-steady flow was obtained from 
numerical integration by Stoolman. Analytical solutions 
for Stoolman’s original equations were presented by Chang 
and Hsu (12) and Edwards (13) with an additional assump- 
tion of linear axial velocity distribution. It was also indi- 
cated (12), contrary to Stoolman’s conclusion (10), that a 
self-excited shock oscillation is not possible without introduc- 
ing viscous dissipation. 


Present Problem and General Approach 


The main purpose of this paper is to show that the in- 
stability of a supersonic inlet diffuser is essentially dependent 
on the viscous dissipation of the flow, and to justify and unify 
the various existing theories of instability for supersonic inlet 
diffusers. 

The subcritical operation of a typical supersonic inlet dif- 
fuser is shown in Fig. 1, where, for a small disturbance gen- 
erated downstream of the shock system, the instantaneous 
shock positions and captured flow boundaries are shown with 
solid and dotted lines, respectively. Assume the free stream 
captured flow area to be always inside the circle of triple- 
shock intersection. The Mach number M, immediately be- 
fore the normal shock may be considered approximately uni- 
form by averaging the Mach number distribution in the coni- 
cal flow field because the variation is rather small. Therefore, 
the diffuser model under consideration can be simplified to a 
normal shock inlet diffuser as shown in Fig. 2. The viscous 
foree acting on the surface of the central body and the free 
stream boundaries in the external diffusion regime are still 
retained in the present analysis. The length of the plenum 


* Numbers in parentheses indicate References at end of paper. 
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Normal shock inlet diffuser at subcritical operation 
(simplified model ) 


chamber is assumed much longer than the portion of inlet and 
the exit nozzle. Therefore, the flow in the inlet and exit 
nozzle can be assumed to be operating quasi-steadily (11). 
The neutral stability criterion obtained by Stoolman (10) 
and Mirels (11) for small Mach number in the plenum cham- 


ber is 


This criterion indicates that the particular mechanism lead- 
ing to instability is actually the amplification of acoustic 
waves in the plenum chamber after successive reflections from 
the inlet and exit sections. In the first bracket of Equation 
[1] is the amplitude ratio of the reflected wave to the imping- 
ing wave at the exit nozzle; it is always less than 1. The 
second bracket contains the amplitude ratio of the reflected 
wave to the impinging wave at the inlet of the plenum cham- 
ber. This ratio is possibly greater than 1 if the viscous term 
is considered in the equation of motion. Obviously instabil- 
ity may occur if |¢s|> [(y — 1)/2]M.. Therefore, the sta- 
bility problem simply reduces to finding the acoustic imped- 
ance in the inlet of the plenum chamber. In order to obtain 
the inlet acoustic impedance, it is necessary to solve the vis- 
cous compressible flow equations in the external and internal 
diffusion regimes. 


[1] 
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Analysis of External Diffusion Regime 


Fundamental Flow Equations 


In the quasi one-dimensional unsteady and viscous com- 
pressible flow,® the equations of motion in the external dif- 


fusion regime are taken, respectively, as follows: At 
Continuity 
re) 
— (pA A)=0 
(pA) + (puA) 
fu’ ose’ p’ u 
= | — s—{— —-1)>-2 
+02 (Z) +5 2 (2) 5 
Motion 


D 2 


act 
| 
| 
~ 
w 


Entropy 


‘ o\ Ss fey 1) 
ot : or} 2D 


The conventional equation of continuity (3) in one-dimen- 
sional unsteady flow is used here.6 The change of flow area 
with respect to time is taken into consideration as a flexible 
wall tube to account for the mass spillage over the cowling. 
The transient effect of viscous dissipation, such as flow sepa- 
ration, shock boundary layer interaction, etc., will be taken 
approximately into account by c;(z, t). Equation [4] expresses 
that the Eulerian derivative of entropy is contributed essen- 
tially by those viscous dissipations (14). 

Assume the flow to be slightly perturbed from the steady 
state and represent the small perturbations as 


= p(x) + p’(z, t) 
= u(x) + u’(z, t) 
= p(x) + p(z, t) 
= Sx) + S'(z, t) 
= G(x) + cy'(z, t) 
= A(x) + t) 
= D(x) + D(z, t) 

6] 


Equations [2, 3 and 4] together with Equations [5 and 6] 
can be reduced, respectively, to the following forms by re- 
taining only the first-order terms in perturbations 


5 The justification of quasi one-dimensional assumption has 

— by Trimpi (6) in his reply to Dailey’s comment 

6 The mass flux term plv — u(dR/dx)|2xR in Stoolman’s 
equation of continuity (10) is equivalent to p(0A /d¢t) in Equation 
[2]. The approximation of the mass flux term made later by 
Stoolman is equivalent to the neglect of v with respect to u(0R/dz). 
It seems too far from the realistic model. Consequently, the 
solution of v in Stoolman’s original equation retains a singularity 
(10,12) atA = 0. 


Eliminating the terms involving ¢, between Equations |S 
and 9] and replacing p’/p through Equation [10], Equation 
[8] becomes 


dx ¢y dx (y — oz \er — 1)a? Ot \ cy 


where 4? = yp/p. Equations [7 and 11] are the two first- 
order differential equations with two dependent variable: 
u’'/a and p’/p if S’/c, and A’/A are given. Because of lin 
earity of the previous equations and the conditions of Equa 
tion [5], the harmonic form of oscillatory time dependence car 


be chosen 
Se 
p’/p = = [ 


p'/p = ¢(x)e = (ze =[¢(x) — yo(x) [13 
and 


& 
A'/A = {14 


where o, v, ¢, € and 6 are, in general, complex functions of 2 


-and w. Making use of Equations [12, 13 and 14], Equation: 


[7 and 11] can be reduced to the following system of ordinary 
differential equations with dependent variables in o and v 


+ iw( +6)=0 [15] 
"=" te 


( di iwa? 


dz — 


a 
game 


where 6 and € are assumed to be given functions, as is shown 
later. 


Linear Velocity Distribution 


In view of Stoolman’s experiments [Fig. 10 of (10)], the 
axial velocity gradient is known to be usually very large in 
the external diffusion regime just preceding the flow insta- 


bility. It is probably reasonable to assume that # decreases 
linearly with z. Therefore, we choose 
di = a = a* [17] 
de z* 
| 
_[p p’ 2u’ Cy’ : 
2D p D p u Cy 


Moreover, we introduce a new independent variable z 


2 2 
(2) = i? <1 18] 
Zz a 


in terms of which we have 


ol or] cy 2D 7) 
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also we shall define a reduced angular frequency 


wa* 
>0 
a* dii/dx 


[20] 


Substituting Equations [17-20] into Equations [15 and 
16] we obtain, respectively 


do dv dé 
— — — 9 
2z iB(o + 6) = [21] 
— be) do, dv 


where b (y — 1)/(y + 1). Eliminating v in Equations 
21 and 22], we have a second-order differential equation in 


2(1 — 2) (2 — #601 -+ B — iB)\(1 — 
dz? dz 4 
(1— b)(1 — bz) d% 
—-(2 


Equation [23] is a nonhomogeneous complex hypergeo- 
metric equation, with singularities at z = 0, 1 and © which 
correspond to M = 0 and 1 in the actual flow field. 


Quasi-Steady Flow and Constant 6 


It may be possible to solve Equation [23] when e(z) and 
6(z) are given. However, in this paper, only the simple case 
of quasi-steady flow is treated. In other words, @ is con- 
sidered so small that at each instant those nonsteady per- 


turbation ratios p’/p and u’/a may be represented approxi- 
Equation 


mately as their amplitudes o and », respec tive ly. 


then reduces to 


ao do 
(l-2)— = 
2) dz 
de 1—bzed d6 d*6 
+2 —| [24 


The area fluctuation function 6(z) is rather difficult and 
complicated to impose. However, in view of Stoolman’s 
experimental results [Fig. 7 of (10) ] showing that the effect of 
external diffusion on the flow instability is about the same as 
internal diffusion, we may assume d6/dz = 0 for a first ap- 


proximation. Equation [24] then reduces to 
, do de 1 — be d*e [25] 
dz? “dz dz 2b dz? 


The significance, and the consequence, of this assumption 
should be examined. The more realistic fluctu: ating model 
has the boundary streamlines swinging about a pivoted point 
near the cowling lip as shown in Fig. 2. Thus, d6/dz is a 
positive quantity, whereas de/dz must be a negative quantity 
in order to produce flow instability, as will be shown later. 
Assume that the second derivatives of 6 and e€ are much smaller 
than those first derivatives and can be neglected in Equa- 
tion [24]. Therefore the flow area fluctuation in the external 
diffusion regime has, in fact, a damping effect on the flow in- 
stability. For this reason the assumption of constant ampli- 
tude will not deteriorate our qualitative investigation of the 
‘iscous induced instability. It should be noted, however, 
that the present result will lead to a smaller value of € for 
neutral stabilities. 


May 1960 


oscillating shock that € = e for z = z2, we _ 
€ 
22 — 
for < 2s. 
iB(1 — b)? iB(1 — b)*%1 — bz) | de 
1 — b)z 
1 
— b)6 + (1 — — +1 - [23] 
dz dz? 


Linear Entropy Perturbation €(z) 


For a qualitative analysis, the perturbation entropy dis- 
tribution is assumed to be linear with z, primarily for mathe- 


matical simplicity 
d ( € ae 
dz 


where £/# and m are constants. Integrating Equation [26] 
and inserting the boundary condition immediately behind the 


[26 ] 


The first term on the right-hand side of Equation [27] is the 
perturbation entropy due to shock oscillation, which remains 
constant with fluid particles moving downstream. The 
second term on the right-hand side is the perturbation entropy 
due to viscous een which is assumed to increase lin- 


sarly with the decrease of z _ 7 


Substituting Equation [26] or [27] into Equation [25], we 
have 


Solutions and Boundary Conditions 


2 


Solution for o can be easily obtained in terms of M* by 
using Equation [18] and Prandtl’s shock relation M,*/.* = 1 


ue: 


where 


u 
Us 


Solution for »y may be obtained by substituting or [29] 
into Equation [22] in which 6 is taken to be zero 
BM,*? 


| 1 y+1- 
+ 


The arbitrary constants B and C can be determined from 
the boundary conditions immediately behind the oscillating 
normal shock. Consider a normal shock oscillating slightly. 
The upstream flow conditions remain constant except for the 
magnitude of the relative velocity, whereas the flow condi- 
tions behind the shock are disturbed and vary with the shock 
velocity ratio Assuming w/t. <1 and using Equation 
[6], Stoolman (10) obtained the following relations from Bur- 
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[30] 


| 
| 
a 
| 
| 
a- 
5 
8 
— 
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gers’ general expressions (15) by retaining only the first-order 
terms in perturbations (12) 


i (: + a) [31a] 
w/ te + 1)M;? 


w/in (2yM,? — y + 1 M,?) 
[31¢] 


(S'/co), 


y+12yM2 


w/the 


where w is the shock oscillating velocity, positive toward 
urstream.’ It is important to note that entropy will increase 
if the normal shock is mov ing upstream due to any disturb- 
ance, such as isentropic compression waves originating from 
the plenum chamber. 


7 It is equivalent in that the free stream velocity is accelerated 
to u + w with respect to a stationary shock. Therefore, it is in 
accordance with our coordinate system that flow velocity is posi- 
tive toward downstream. In this way, ¢«/(€/i#2) is expressed as 
a positive quantity. The present sign convention differs from 
that in (10 and 12). 
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Fig. 3 Velocity perturbations in the subsonic diffusion regime 
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Fig. 4 Density perturbations in the subsonic diffusion regime 
(quasi-steady) 


arbitrary constants B and C can be determined. 
the complete solutions for vy and o can be expressed as 


v= M,* 2 2 
_& +o fi- (y= EL 
= y+1 y-1 af 
1 
= 2 


with the relation given in Equation [13], we obtain 


1—? {ro — 1) 


Equations [32 and 33] are plotted in Figs. 3 and 4, respe: 
tively, against velocity ratio X with M, and m as the p: 
rameters. Also g/v is plotted in Fig. 5. Since the ink 
acoustic impedance can be expressed in terms of vy and ¢ 


(5). 


it is clearly indicated in Fig. 5 that flow instability will no 


. occur for a nonviscous fluid, since ¢ < 0 for m = 0, and that 


instability may occur if m is sufficiently large. 


Internal Diffusion Regime 


It was noted in Equation [25] that the terms involving th 
flow area fluctuation disappear when the assumptions o! 
quasi-steady flow and constant 6 are made. In other words 
the flow equation in the external diffusion regime becomes ex- 
actly the same as in the internal diffusion regime. Therefore, 


the solutions in Equations [32, 33 and 34] can be directly 
extended into the internal diffusion regime, provided the as- 
sumptions of linear velocity variation and parabolic pertur- 
bation entropy distribution, with respect to the distance, are 
also used. 
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Fig. 5 Modified acoustic impedance in the subsonic diffusion 
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With the boundary conditions given in Equations [31], the 
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Final Results 
Equation of Neutral Stability Limit 


Ah 
The stability criterion in Equation [1] can be further ex- 
pressed as follows, by using Equation [35] and 4/,* instead 
of M, 


The perturbation entropy generated due to viscous dissipa- 
tion at the inlet of plenum chamber is taken directly from the 
last term in Equation [27] and may be expressed as 


m 


Substituting Equations [32 and 34] into Equation [86] and 


using Equation [37], we obtain the equation of neutral sta- 
bility limit in terms of A, and e€,/(&/ie) 


(anes 


[37] 
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Fig. 6 Stability boundaries of viscous dissipation entropy with 
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For a given free stream Mach number r2/(§/i2), o2/(E/i2) 
and €/(¢/i#.) are known quantities in Equation [31]. Thus 
the critical plenum chamber velocity ratio \, can be solved 
in Equation [38] for any given values of €,/(&/#2) or vice 
versa. 


Stability Boundaries 


From Equation [38], we plot the critical values of €,/(£/i2) 
against A, for neutral stability in Fig. 6. Also, stability 
boundaries are replotted in Fig. 7 as lines of constant €,/(&/i2) 
with M7.* and M;,* as the ordinate and abscissa, respectively. 
The value of €,/(£/i#2) on each curve may be considered as 
the sensitivity of a diffuser for neutral stable operation. In 
other words, the described value indicates a certain shock 
velocity which has to spill a certain amount of mass flow, due 
to a small entropy generated by viscous dissipation, so that a 
neutral stability can still be maintained. The neutral condi- 
tions previously mentioned can only be achieved if the diffuser 


Y- 


1 


2 


(y= &u 
Vy + 1) 


operating conditions (M,*, .7.*) lie on those curves of Fig. 7. 
Stable or unstable conditions will be obtained if the diffuser 
operating conditions are above or below those curves respec- 
tively. It should be noted that all the stability boundaries 
in Figs. 6 and 7 are valid only for small 17, according to 
Equation [1]. For a further understanding of the stability 
boundaries, we plot [e,/(&/#2)]A. = 0, €s/(€/i#2) and also 
their sums vs. M,* in Fig. 8 for the particular case of zero mass 
flow. The main purpose is to show that the entropy generated 


*? 


due to shock oscillation alone is insufficient to cause 
instability. It also indicates that the maximum value of 


léu/(E/ti2) |Xce = O required for instability occurs near 


M,* = 1.36 or (M, = 1.5). 


Comparison With Stoolman’s Experimental Data 


Discussion 


Stoolman (10) tested four normal shock diffuser models 
with various internal contraction ratios (43/A4 = 0.28 to 1.0) 
at M, = 1.6 (reported at 1.8), 2.0 and 2.45. The Reynolds 

-numbers based on diffuser length are all about in the order 
Comparisons of Stoolman’s data with the present 
analysis are made in the following discussion. 
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Fig. 8 Stability boundary of total entropy generation for zero 
mass flow 
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No significant effect of diffuser geometric contraction 
ratios on stability limit 


For a given free stream Mach number and captured mass 
flow, Stoolman’s result does not reveal any significant effect 
of diffuser geometric contraction ratios on the stability limit 
[Fig. 7 of (10)]. It immediately confirms a few of our im- 
portant assumptions: Stoolman’s result indicates that the 
effect of external diffusion on stability limit is approximately 
the same as internal diffusion. In other words, the assump- 
tions of quasi-steady flow and constant amplitude of flow area 
fluctuation are quite close to the actual flow conditions; 
also, Stoolman’s result shows that the velocity distribution 
in the external and internal diffusion regimes is not very 
significant to the stability limit. Therefore a simple linear 
velocity distribution may be chosen for qualitative investiga- 
tion. 


Agreement with Stoolman’s data on neutral stability 
boundary 


Stoolman’s experimental data of the neutral stability limits 
for M, = 2 and 2.45 are also plotted in our Fig. 7. It is 
interesting to note that the previously mentioned data fall 
exactly on the curve of €,/(&/#e) = 0.099. If we assume that 
the diffuser sensitivity ¢€,/(E/i#2) for Stoolman’s model is 
close to this value, then Stoolman’s observation of stable op- 
eration at M, = 1.6 or 1.8 with zero mass flow (Fig. 7) indeed 
confirms the present analytical result. a 


Comparison With Stability Theories of Trimpi and 
Sterbentz-Evvard 


The present result may also be interpreted as follows: For 
a given free stream Mach number and captured mass flow, 
the flow stability is dependent on the instantaneous ratio of 
the total perturbation entropy generated by viscous dissipa- 
tion and shock oscillation, to the corresponding shock veloc- 
ity. Since the shock velocity is directly related to the in- 
stantaneous captured mass flow, the present analysis indeed 
confirms Trimpi’s proposed theory of instability (4). 

In the Sterbentz-Evvard analysis, a ratio dp,-/dm was in- 
troduced which, as stated in their paper, represents the pres- 
sure change in the ramjet due to the diffuser shock position 
which changes during the “unsteady” operation of the dif- 
fuser. It has to be justified whether this ratio should be in- 
terpreted as the ratios between the “instantaneous” or the 
“time averaged”’ values of them in an unsteady flow with fixed 
exit throat area, or between their ‘“‘steady”’ values with small 
changes of exit throat area. In order to answer this we ex- 
plore the meaning of the inlet acoustic impedance, which is 
shown in Equation [35] to be independent of time. The 
acoustic impedance actually represents the magnitude of the 
amplitude ratio of the perturbation quantities ¢ and v either 
in an “unsteady” state, or in a “steady” state with a small 
change of exit throat area. For this reason Mirels (11) can 
further express the stability criterion in terms of the “steady”’ 
diffuser characteristic curve (m/poc) (dfoc/dm) directly from 
the relation with the inlet acoustic impedance. However, it 
is not true when the stability criterion is measured by the 
slopes of “time average’’ diffuser characteristic curve during 
an “unsteady” operation. Therefore, it is clear that the 
original sense of dp,-/dm in Sterbentz-Evvard’s analysis is 
exactly the same as Trimpi’s idea, but the interpretation of 
this ratio as the slope of diffuser characteristic is only appli- 
cable in the “steady”’ range of subcritical operation. The ar- 
guments noted also explain why some experimental results 
(5,16) have shown diffuser instability at zero or negative slopes 
of the “time average” diffuser characteristics, and why some 
others may have intermittent buzz in the subcritical range. 
This is because the instantaneous total pressure recovery with 
respect to its corresponding mass flow cannot be expressed on a 
time average scale of the diffuser characteristic curve. — 
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Nomenclature 

A(a, t) = instantaneous flow area captured by the inlet 
diffuser 

a = local sound velocity 

-B,C ~~) = arbitrary constants 

b - = (y — 1)/(y + 1), constant for gas property 
(6 = for air) 

b) = 1/p(u?/2), coefficient of friction 

Cp, Co = specific heat at constant pressure and volume, 
respectively - 

D(z, t) = 4A/P, hydraulic diameter of flow area | 

M = u/a, Mach number 4 

M* = u/a*, dimensionless velocity 

m = slope of €/(t/t:2) vs. z curve, m > 0 i 

m(t) = instantaneous captured mass flow 

P(2, t) = perimeter of A(z, ¢) 

p(x, t) = instantaneous pressure 

R(2, t) = radius of mean flow area A(z, t) of 

r = radial coordinate normal to diffuser axis 

S(a, t) = entropy 

t = time 

u(x, 0), v(x, t) = instantaneous velocity components in 2, 7 di- 
rections, respectively 

w = £et!t, shock oscillating velocity (positive to- 
ward upstream ) 

x = coordinate along the diffuser axis (positive 
toward downstream) _ 

Zz = (x/x*)? = (a/a*)? = M*?, new independent 
variable 

B = —war*/a* > 0, reduced angular frequency 

ratio of gas specific heats; y = 1.40 for 
air 


- = time average value, function of z only 


Conclusions 

1 The determination of the self-excited shock oscillation 
of a supersonic inlet diffuser is not possible without considera- 
tion of viscous dissipation. The entropy generated due to 
the shock oscillation alone is not sufficiently large to cause 
flow instability. 

2 Instability is found to be due to further entropy gener- 
ated by viscous dissipation. The present stability parameter 
(or diffuser sensitivity) expressed as €,/(£/i2) is the amplitude 
ratio of the perturbation entropy generated by viscous dissi- 
pation, to the corresponding shock velocity. 

3 Neutral stability boundaries are obtained in terms of the 
diffuser sensitivity ¢,/(t/i#) for the entire range of free 
stream Mach numbers, and are in agreement with Stoolman’s 
experimental data. 

4 The free stream Mach number, which is rather insensi- 
tive to instability, is in the neighborhood of 1.5. 

5 The present analysis confirms Trimpi’s theory of insta- 
bility and also unifies the Sterbentz-Evvard theory with that 
of Trimpi. However, the interpretation of Sterbentz- 
Evvard’s theory as the slope of diffuser characteristic is only 
applicable in the “‘steady”’ range of subcritical operation, but 
is not applicable in the range of the time averaged values of 
an unsteady state. 


v(x), o(x), g(x) = time independent parts of u'/d, p’/p and p’/p, 
respectively. See Equations [12, 13] 
g(x) — yo(x), time independent part of 
See Equation [13] _ 
5(x) = = time independent parts of A’/A. See Equa- 
tion [14] 
(4/yp)(p'/u'), acoustic impedance 
instantaneous density 


t) = shear stress 

= M*/M,* = t/t, velocity ratio 

w = angular frequency 

= amplitude of shock velocity 


Subscripts 


* = sonic condition 


| 
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/ = nonsteady small perturbation, function of z 
and ¢ 

ce = sections refer to different locations in the dif- 
fuser. See Figs. 1 or 2 

¢ 7 = plenum or combustion chamber 

o = local stagnation condition 

s = due to shock oscillation 

. = due to viscous dissipation 
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volume of experimental data has been obtained on the crater- 
ing effects of hypervelocity particles, and the results of these 
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mates to yield conservative figures for the probability of pene- 
trating the skin of a satellite. 
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Satellite Pressure Losses 
Caused hy Meteoroid Impacts 


In order to predict the frequency of hull penetration of a satellite capsule by met 
are made of the frequency of encounters with meteoroids, and the cratering effect of each impact. 
The cratering effects are based on the correlation of recent laboratory experimentation using hyper- 
velocity particles, and very conservative estimates of impact frequency are employed. The resulting 
calculations of percentage of hull area covered by holes, and loss of internal pressure ys. time, are ex- 
pected to be conservative for design purposes. 
sule having a given surface area and time of exposure to meteoroids results in predictions of: Prob- 
ability of penetration of several hull designs, time required for air pressure to drop to 4 atm, and 
weight of reserve air needed to replace leakage and maintain internal pressure at 1 atm. ae 
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eoroids, estimates 


Application of this method to a manned satellite cap- 


Manned space vehicles pose the problem of maintaining air 
pressure within the hull. This problem is examined from the 
following standpoints: 

Thickness of hull sufficient to withstand meteoroid bom- 
bardment for the duration of the mission of the space 
vehicle. 

Weight of air needed to replace leakage when this does 
occur. 

Time required for air pressure to drop below a tolerable 
level. 

This last item is of practical interest to the passenger, who 
must decide whether he has time to complete his mission, or 
whether he must return to Earth immediately. 
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Calculation of Hull Penetration 


Laboratory data on hypervelocity toate (2) — the 
empirical formulas* 


_ The assumptions and extrapolations (velocity and size) 
necessary in order to arrive at Equations [1 and 2] are con- 
sidered in (2). Note that the expression for depth of crater 
is expected to be correct within a factor of 2 or 3 for untested 
materials. A further assumption is that the skin of the satel- 
lite behaves as a semi-infinite medium, i.e., craters on a thin 
skin are the same as on a thick plate, and the skin is pene- 
trated when h exceeds the thickness of the skin. 

The meteor bumper configuration proposed by Whipple 
and others, consisting of a bumper shield placed some dis- 
tance from the main pressure hull, is known to be quite ad- 
vantageous on a weight basis. The thin bumper shield serves 
to shatter the hypervelocity particle, so that the material 
sprays out from the hole in the bumper and covers a large 
area on the main hull. The loss in concentration of energy 
permits the use of thin hulls, so that a saving of total weight 
is effected. Unfortunately, there is virtually no experi- 
mental performance data for double wall construction at 
present, and the following estimates must be based on solid 
hull construction. 


* Equations [1 and 2] are valid with any consistent set of units. 
The dimensions given in the Nomenclature apply to subsequent 
calculations. 
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Probability of Encounters With Meteoroids 


The estimates of frequency of encountering meteoroids at 
each kinetic energy level have a spread of about three orders 
of magnitude, depending on the source of the estimate [see 
(3, 4); note that the figures in (1) fall near the average of the 
highest and lowest estimates]. Fig. 1 presents these esti- 
mates. Although the recent rocket and satellite meteoroid 
impact data (which fall in the meteoroid size range below 
that which may be expected to penetrate satellite hulls) 
show high impact frequencies, authoritative expectations 
place the larger meteoroids near the lowest curve (Watson 
curve). Noting, however, that meteor showers may increase 
impact frequency considerably (5), the most conservative 
estimates of Fig. 1 are used, on the kinetic energy basis 


N; = 1/U [3] 


where N; is the number of impacts having kinetic energy 
greater than U (in.-lb) per 10 m? of exposed area per hr. No 
correction for satellite attitude appears to be warranted at 
this time. 

Combining Equations [2 and 3] 


N, is the number of holes per hr in a 10-m? hull of thick- 


where 
ness h. 


Area Covered by Holes 


For purposes of calculating the loss of air within the space 
vehicle, each meteoroid crater of depth greater than hull 
thickness h, is assumed to be a sharp edged orifice of area 
mh*(h > h). Fig. 2 illustrates this assumption, which is 
estimated to be conservative; i.e., it should lead to high 

Re estimates of air leakage. The total area covered by holes 


larger than is 
a= 
Nh 


[5] 


15,500E \E, hy hy 


per hour. Expressing the hole production rate 
a 
— = 
A 1 


_ with ¢ in hours, the penetration constant is 


E,\o.73 1 
C, = 4.86 X 10-° (2) (7 in.?/in.? per hr 
1 N 


Thickness hy must be taken as infinity if all meteoroids are 
to be considered. Cutting off the calculation at probability 
level N;, as in Equation [4] 


2 


E 0.09 
hy = [9] 


Values of hy are shown in Table 1 for several values of N, and 
E, and Fig. 3 presents penetration constant C;, as in Equation 
[8]. N> is shown for a 28-hr mission. 
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Table 1 Values of hy 


N, E, psi X 10~* ~ 
28 hr 1 5 10 15 30 45 
0.001 0.614 0.415 0.351 0.318 0.268 0.244 
0.01 0.285 0.193 0.163 0.148 0.125 0.113 
0.1 0.132 0.0900 0.0760 0.0684 0.0581 0.0528 


Pressure Decay 


“or an infinite pressure ratio across a nozzle of area a 


aw ( 2 | 207 
— = —aP 
dt +) 


for y = 1.4 = 
dw ae a P 
— = —0.578 AP ¢ 9) 11 


Assuming the expansion to be a throttling (Joule-Thomson) 
process, the internal temperature remains constant 


(12 
dt Wo dt 
Py 0 


Combining Equations [11 and 12] 


dP 


oY APo ) 


Using the convenient spherical shape 


4 
— 
r e 
A = 4rr? 
Equation [13] now becomes a 
dP 3 
— = —0.578 (2 ygRT ~ at (14) 
r A 


Noting that VW ygRT = C, is equal to the speed of sound in 
air at temperature 7’, and introducing Ct equal to a/A from 
Equation [8], Equation [14] becomes 


dP "oC; 
771.734 tdt [15] 
Integrating, we arrive at 
P 
0-867CoCit /r [16] 
0 


Fig. 4 consists of plots of P/Pp vs. time, for various values 
of CoC/r, as in Equation [16]. Note the interesting feature 
of Equation [16]; the larger the space vehicle, the slower the 
pressure drops. This behavior is occasioned by the geometric 
fact that surface area to volume ratio decreases as size in- 
creases. 

Of practical interest is the time for internal pressure to drop 
from 1 atm (at zero time) to } atm 


aa 
0.867 =log3 = 1.10 
ty = 0.969C,-1/2 


where C, = 0.74 CoC,/r in units of (hr)~?, (C2 = 10°C; for 
70 F andr = 3'). Equation [17] is plotted in Fig. 5. 
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Estimates for a Manned Satellite Capsule 

It is of primary practical interest to the passenger of a 
manned satellite capsule to learn the time required for the 
pressure to drop below a tolerable level. Fig. 6 is obtained 
by taking 3 atm as the tolerable pressure level, using the values 
of penetration constant C, obtained from Fig. 3 for various 
values of 4, FE and N» (see Eq. [9]), using C2 = 10°C; which 
implies that T = 70 F and r = 8 ft; and, finally, calculating 
time from Equation [17]. Note, in Fig. 6, that even the 
very thin hulls (e.g., 0.010 in.) permit ample time for return 
to Earth. 

In order to compare heavy hulls (with no penetrations) with 
light hulls permitting pressure drops, we compute the weight 
of air needed to replace leakage and maintain pressure at 1 
atm. Equation [11] at constant pressure becomes 


dw Ty 
dt RT A RT 
Integrating, we have 
AW = 0.289 gy/RT- [18] 


Expressing weight loss as a fraction of Wp « and : adopting the 
spherical shape such that AP)/W,) = 3RT/r, Equation [18} 


becomes 
AW C,t? 
= 0. — = 0.867 [19] 
Wo r 
For Co = (1128) (3600) ft per hr, r = 3ft,¢ = 28hr 
AW/W, = 9.19 X 10°C, [20] 


Using Equation [20] and Fig. 3 for C;, one may rapidly esti- 
mate air weights AW (noting that Wy = 8.50 lb for r = 3 ft, 
P) = 1 atm, T = 70 F) for various values of shell thickness 
E and probability N;,. 

When the sum of weight of replacement air plus hull weight 
(10-m? area) is plotted vs. hull thickness, it is found that there 
is a hull thickness which represents minimum total weight. 
The dashed lines (no. 1) in Fig. 7 represent those hull thick- 
nesses, with elasticity moduli of beryllium, steel, titanium 
and magnesium having been taken as 45 X 108 psi, 30 X 10° 
psi, 15 X 108 psi and 5 X 108 psi, respectively. Dotted lines 
(no. 2) represent thicknesses for P dropping to } Pp in 28 hr, 
and solid lines (no. 3) represent hull thicknesses for no penetra- 
tions in 28 hr. 

Fig. 8 shows the corresponding weights, with curves num- 
bered: 

1 Shell weight for minimum total weight (shell plus air), 
which is attained by replacing leakage to maintain P = P». 

2 Total of shell plus air w eight for condition 1. 

3 Allowing P to drop to 3 Po in 28 hr with no air replace- 
ment. 

4 No inet in 28 hr. 


Conclusions 


Calculations of pressure loss vs. time, as well as satellite 
capsule hull weights required for several types of protection 
against pressure loss, have been based on the most conserva- 
tive estimates available for impact frequencies. Besides 
being conservative on impact frequency, the calculations do 
not reflect the expected advantages of the meteor bumper 
configuration, nor do they account for schemes such as self- 
sealing walls which should provide tremendous advantages 
in saving weight. On the other hand, much more empirical 
data on impact probabilities are required in order to place 
confidence levels on the estimates. It is anticipated that 
impact frequency data will be forthcoming in the next years 
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of satellite research, as well as performance data on meteor 
bumpers and self-sealing double hull construction. 

Even with the most conservative calculations, an examina- 
tion of Fig. 6 shows that “explosive decompression” is not 
properly descriptive of the slow leakage process. Only if the 
pressure hull has too high a design stress will there be danger 
of catastrophic loss of pressure on meteoroid impact for any 
reasonable probability level. 

The thickness and weight estimates of Figs. 7 and 8 show 
that weight saving is achieved by allowing meteoroids to 
puncture the hull, and by replacing the leakage, at least for 
the example of A = 10 m?, r = 3 ft, T = 70 F with a 28-hr 
mission. Whether it will be considered practical to permit 
hull penetration is questionable. Also, that other capsule 
geometries and missions may lead to other conclusions on 
weight is not considered. General conclusions on capsule 
design will be possible after extensive design studies are made, 
using such methods as have been described. 


Nomenclature 


= total area, holes produced by meteoroid impacts _ 


A = surface area exposed to meteoric bombardment 

» = speed of sound in air 

, = penetration constant, (Aa/A)/hr 

C, = leakage constant, 0.74 hr-? (70F, C. = 10°C, 
forr = 

E = modulus of elasticity of hull material, psi 

Ey = reference modulus, 108 psi 

4g = acceleration of gravity at Earth surface, 32.2 fps? 

h = depth of crater, see Equation [2] 


deviation. 


HE TEMPERATURE distribution on the surface of an 
Earth satellite or space vehicle is receiving considerable 
attention. Schmidt and Hanawalt (1)* discussed the case of 
a nonrotating cylindrical shell exposed to solar and Earth 
radiation. They pointed out the importance of the ratio of 


Received May 25, 1959. 
1 Professor of Applied Mathematics. 
? Professor of Mechanical Engineering. 


3 Numbers in parentheses indicate References at end of paper. 


May 1960 


hy = hull thickness which would ‘es ave ro penetrated by 
largest meteoroid (defined by impact probability Nz) 


N; = number of impacts with meteoroids larger than a given 
size 

N», = number of impacts with meteoroids larger than that just 
able to penetrate thickness h 

P = internal pressure, any time 

P, = initial pressure, time zero 


r = radius, spherical capsule 

R = gasconstant 

t = time, hr 

YT = air temperature inside capsule 

U = kinetic energy of impacting particle, in. db 


V = volume of crater, in.® 


W = weight of air inside capsule 
W, = initial weight of air inside capsule 
y = ratio of specific heats, taken as 1.4 
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Solar heating in a space vehicle idealized as a thin walled circular cylinder rotating with uniform 
velocity about its geometric axis is studied for a situation in which heat transfer by convection and 
heat exchange within the cylinder is negligible. The nonlinear problem is approximated through 
a perturbation analysis and detailed estimates are made of the parameters of interest for various 
ranges of speed of rotation. An estimate of the error between the exact and perturbational ap- 
proach is made in the case of an illustrative example which might be expected to entail an extreme 


absorptivity to emissivity on the temperature of the shell. 
Seavery (2) and Hibbs (3) studied the temperature problem 
of orbiting satellites; Haas and Drummeter (4), Drummeter 
and Schach (5), Hanel (6) and Sandorff and Prigge (7) in- 
vestigated the problem of temperature control of satellites 
and space vehicles. 

The problem considered here is that of a thin walled circular 
cylinder rotating with uniform velocity about its geometric 
axis and exposed to solar radiation. When steady state is 
achieved, the heat received by the rotating body is equal to 
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the amount of heat reradiated into space. On an application 
of this analysis to a space vehicle, the heat transfer by convec- 
tion can be neglected as soon as the vehicle leaves Earth’s 
atmosphere. Furthermore, it is assumed that no heat ex- 
change takes place inside the cylinder. 

Simple physical considerations indicate that the part of the 
cylinder receiving the radiation will be hotter than the 
“shadow side.” This temperature difference will be largest 
for the nonrotating cylinder and will decrease for higher rota- 
tional speeds. There should also be a displacement in the 
location of the line of the highest and lowest temperature on 
the surface of the cylinder due to rotation. For the stationary 
cylinder, the highest and lowest temperatures should be lo- 
cated at points that are closest and farthest from the heat 
source. In the case of a rotating cylinder, these positions 
should be shifted in the direction of the rotation. An analysis 
of the problem confirms these predictions and gives quanti- 
tatively the temperature distribution along the circumference 
of the cylinder. 

The parameters defining the problems are: 


= radius of the cylinder, ft 
= wall thickness (assumed to be small compared with the 
radius), ft 


k = heat conductivity of the wall material, Btu/ft hr R 
c = specific heat of the wall material, Btu/Ibm R 
p = density of the wall material, lbm/ft? 
a = thermal diffusivity of the wall material, ft?/hr 
s 
Fig. 1 Thin walled cylinder receiving solar radiation. Solar 


coordinate system v, 7 stationary in space 


— 


Fig. 2 Heat flux (per unit cylinder length) traveling along the 
developed stationary cylinder surface 


Fig. 3 Heat balance for a volume element of the cylinder wall 
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; 


a = average (with respect to wave length and angle of in- 
cidence) absorptivity of the wall material 

eé = average total hemispherical emissivity of the wall ma- 
terial for the spectrum of wave length radiated by the 
cylinder 

K, = energy received from the radiating source by a plane 
perpendicular to the “line of vision,’’ Btu/ft? hr 

v = circumferential velocity, ft/hr 


It is assumed that the cylinder axis is perpendicular to the 
line connecting the center of the cylinder cross section with 
the heat source, and that the cylinder length is large compared 
with the radius and small compared with the distance between 
the cylinder and the heat source. It is further assumed that no 
heat flow takes place in the axial direction, and therefore all 
calculations are made per unit length of the cylinder. 


Differential Equation of the Temperature 4 


Distribution 


In steady state, the temperature at every point on the 
circumference of the cylinder, as determined by the angle 277 
in a coordinate system that is stationary in space, remains con- 
stant (Fig. 1). 

The problem can be visualized advantageously by unrolling 
the surface of the cylinder, as shown in Fig. 2, considering the 
surface stationary and allowing the heat source to travel in 
space, in the direction opposite to the cylinder rotation, in a 
manner that there is no change in the heat input. Fig. 2 il- 
lustrates the heat received (per unit area and unit time) ¢ at 
any time, as a function of the position parameter 7. 

Consider the heat balance during the time element dé for a 
small volume element represented in Fig. 3. The position of 
the element is determined by the coordinate x attached to the 
cylinder. The heat received from the radiating heat source is 
qdx dt. The heat outflow from the volume element due to 
heat conduction is 


re) 


(2) dex dt 


where o@e § 
= temperature 
The heat lost by radiation is 


where ¢ is the Stefan Boltzmann constant. 
The heat stored in the volume element during the time dt is 


ps dx c(dr/dt) dt 
The energy balance yields 


9. 
q dz dt — — dxdt — certdrdt = psc de dt 


It should be noted that the heat input due to radiation is 
moving along the circumference without ‘‘change in shape”’ 
and that at steady state, the temperature distribution is also 
moving along the circumference with no shape change. 
These observations can be expressed mathematically as 


It becomes apparent that a change of coordinates following a 
method described by Rosenthal (8) 
(€ = 2mrn) t=? [2] 
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should be useful. It follows that @ 


or or or ot’ or 
a ot’ dor 
az 
| 
ox? 
or Or OL’ or 
dt deat at’ at aé 
The last relation is due to the fact that 07/Ot’ = 0, because 


in the new coordinate system moving with the heat source, the 
temperature does not change with time. In the new coordi- 
nate system Equation [1] takes the form 


sk 


dn? a dyn sk 


where a = k/pCy. The value of qg is given by Lambert’s law 


q = ak, cos* 217 


where cos* 277 stands for a rectified cosine function as repre- 
sented in Fig. 2. 

The differential equation for the temperature distribution 
is then 


d*r dr 4r?r’aK, 
— — - “cost 2rn [4] 


dyn? a dn gk 


This equation is nonlinear. In order to facilitate the analysis 
it will be assumed that the fractional temperature fluctuation 
T, as defined in the following, is small compared to unity 


=T(1+ 7) 
(5) 


= (3/4 + 7) [6] 


This assumption permits the approximation 
= + 47) = 47 + 1/4) 
Introducing the notation 
T+1/4= 
Equation [4] is replaced by the approximating equation 


= +b aT A 7 
dn? dn 0 cos* [ ] 


where 
6 = 

(168*r2ae /sk)T 

Ay = 


“4 


c 


The 7 which solves [7] is thereby an approximation of 7' 
and yields an approximate value of the temperature 


= (3/4 + T) 


Average Temperature 


The temperature of the cylinder at any point is represented 
as 


where 


= reference temperature 
T = nondimensional temperature parameter, determining the 
temperature fluctuation 


May 1960 


= T(1+ T) 


lowing parameters: 


Let us choose 7) to be the average temperature of the 
cylinder 


It follows that the average value of 7 as given by S, Tdn 
must be equal to zero. 

Consider the equality of heat received by the cylinder and 
heat reradiated. The total heat received is 


Q, = 2raK, 
The heat reradiated is 
dQ. = cer*rd(24n) = 


Introducing Equation [5] into the last equation yields — { 


1 
Qo = 2mrraeT dyn+4 tn) = 2rroTo* 


Jquating Q, = Q2 results in 


T, = VaK,/weo [8] 
Introduction of this relation into Equation [7] results in 
aT dT 


+b — —cT = Ay cost 277 
dn? dn 


b = (2rr/a)v 


= 


—(m/4)c 
The solution of Equation [9a] with the boundary condition 
requiring periodicity 
T dT 
dn 


will yield the temperature distribution along the circumference 
of the cylinder. 

Equation [9a] is already in dimensionless form, but the in- 
terpretation of the solution and the meaning of the coefficients 
of the equation can be enhanced by the introduction of the fol- 


Thermal radius 


= V skT)/l67ak, 
Dimensionless radius 
po = 7/R = Vc 


Thermal velocity 


Dimensionless velocity 


Introduction of Equations [10 to 13] into Equation [9a] 
yields 


— — po? cost 277 [14] 
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Solution of the Linearized Temperature 
Distribution Equation 


Consider first the homogeneous part of Equation (14]. 
The two solutions of this equation are 


where 


= 1+ Vit 
= —1 — + 1/n?) 


In order to solve Equation [9a], it is advantageous to write 
it in the form 


— = A, cos* [17] 


where the subscript 7 indicates the differentiation with respect 
to 7. 

Green’s function for Equation [17] has to satisfy the equa- 
tion 


subject to the boundary conditions of periodicity in the in- 
terval 0 to 1 


The solution of Equation [9a] is then 


1 

T(n) = f, Ae” cos * 2rtg(n, t)dt 
or 
T Ao |= gate 
oLi-e™  1-e" 


em tt — — erm + — cost 2mtdt [22} 
Ao/(v; — v2) can be expressed in terms of vp and po 


Because of the form of cos* 277, the solution for 7 is different 
for the three regions: 


Region I, 0 
Region II, 4,< n< 
Region 


The temperature fluctuation 7 in the three regions is, respec- 


tively 


ev 1 ev 


o I= 16 Vout +1 (1/20)? + 1 (1/4) 


+ 1 e(1/4)v2 (1/4) 


6 = Dirac’s function 


The general solution of Equation [18], as described by 


Friedman (9) is ait, 
o(m, t) = AV(n) + BW(n) + 
= 9) + = 9 
J(W, V) 
where 
A, B = constants to be determined from the boundary 
conditions [19] 
V, W = two independent solutions of the homogeneous 
form of Equation [9a] ra 
H = Heaviside unit function : 


Substitution of computed values of A, B and J(W, V) into 
Equation [20] results in 


1 e717 4+ ert erm 
yt 
1 — e”! 1 — e” 
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where po, %, and are defined by Equations [10, 12, 


and 16]. 


Discussion of the Temperature Distribution 


aw 


The solution of the temperature distribution will be dis- 
cussed for three cases: Stationary cylinder, cylinder moving 
at very high speed, and rotating at intermediate speeds. 


Stationary Cylinder 
When the rotational velocity is zero, the dimensionless 


velocity also becomes zero 


v=0 


vo = 0 = = Po 
The equation for the temperature distribution is ~~ reduced 


to 
or 


Po cosh pon 
(po/2m)? + 1 — e~ (1/4) 0 


1 Po cosh po(y — 1/2) 
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Po cosh — 1) 
e(1/4)p0 _ — (1/4) 90 + 


cos 277 25c] 


The: maximum temperature occurs at n = 0 (the line on the 


cylinder facing the heat source) 


Po 1 __ Po 
Tmax = +3 (po/2m)? + 1 Ee — (1/4) p0 + 


The minimum temperature occurs at 7 = 3 (the line on the 


cylinder on the far side from the heat source) 


Po 1 
Tmin 8 (po/2m) 1 e (1/4) po e— (1/4) po 


The largest temperature difference is 


4 (po/2m)? + 
The tem perature distribution is symmetrical with respect to 
n and can be easily computed for a given value of po. 


Numerical example i 


Consider an aluminum alloy cylinder of: 


radius r = 1 ft 

wall thickness s = 1/200 ft 

thermal conductivity k = 100 Btu/ft hr R 

thermal diffusivity a 3 ft?/hr 

solar constant K, = 442 Btu/ft? hr 

radiation constant o0 = 0.1717.10~* Btu/ft? hr(R)4 


Assume the surface of the vehicle painted black and there- 
forea ~ e = 1. This assumption was chosen to maximize 
the radiation effects. 

The average temperature 7 can be computed — 


= VaK,/weo = 534.25 R 


A = 0.1096 ft 
Po = 9.12 
v* = = 8.71 ft/hr 
For the stationary cylinder 
vo = 0 
be Tmax = 705.48 R 
Tmin = 420.94 R 


Tmax — = 284.54 R 
The temperature distribution is shown in Table 1. 
The error in the energy balance for the case of the stationary 
cylinder was computed because the linearization of Equation 
[4] introduces the largest error in this particular case. 
Employing the computed temperature distribution 7 ob- 
tained from the linearized equation, the reradiated heat (on 
the basis of this approximation) is a oT‘ er2rdn. The heat 
received from the heat source is 2ra K,. Let the error in the 
energy balance be denoted by €g 


1 


The integration of this equation was performed numerically 
for the case of the stationary cylinder, resulting in €¢g = 0.185. 
To obtain an estimate of the error in the temperatur dis- 
tribution one can write a 


= + €7(7)] 
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Fig. 4 Temperature distribution along the wall of a rotating 
cylinder for several nondimensional surface velocities 


© 
8 


6+7/.46° 


SHIFT ANGLE, 


Fig. 5 Position of minimum and maximum temperature for a 
rotating cylinder as function of nondimensional velocity, as 


measured from the line 0-0’ 
where 
T = solution of the original Equation [4] 
€7(n) = temperature error 


The heat balance equation can be then rewritten as 


1 
f, 
fore 
0 


1 
(1 — €g) f, + 4€7()]T4e ~ 2aKk, = 
It follows that 


1 1 1_ 


where, by the mean value theorem, é7 lies betw een €r mini- 
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Table 1 Temperature distribution on the surface of a rotating cylinder for various nondimensional velocities 

vw = 0 705 .48 687 .05 636 . 83 567 .43 500.81 458.49 435.63 424.30 
=1 625.55 579.29 524.37 475.64 445.35 440.54 447 .27 458.28 
vo = V10 551.62 527 .52 505.38 490.51 487 .30 494.16 502.73 512.12 
V30 540.45 526.03 513.51 506.08 506.14 511.50 517.35 523 .52 
%=10 537.46 529.44 523.30 519.29 519.83 521.43 526.77 530.51 

= 0 420 .94 424.94 435.63 458.49 500.81 567.43 05 
= 1 473 .67 493 .27 517.37 548 .67 587 .30 629.83 2.91 
= V10 522.37 533 .57 545.79 559.14 573.71 585.15 72.20 
= 530.01 536.90 597 .46 551.61 559.75 565.02 3.68 
v = 10 534.25 537 .99 542.26 546.00 550.28 552.41 94 


mum and €, maximum. In this particular numerical ex- 
ample é; ~ 5.7 per cent. 


Cylinder Rotating at High Speed 


High speed is defined by the following conditions: v9? > 1, 
and vpo >> ?. Under these restrictive conditions, the equa- 
tion of the temperature distribution can be mgs by 


power series development and by neglecting terms 1/2? and 


as compared to unity. 
The following equations are obtained 


Ti = To 160, E ( 20 (» [26b] 


Po 


sin — ¢) 1286} 


w 


1 (po/2ir)? + 1 


The maximum temperature will occur in the fourth quad- 

rant where 07111/0n = 0. 
Using the same numerical example as in the previous section 
= 9.12, and selecting 


% = 10 v = 87.1 ft/hr = 0.0242 fps oe 
or 0.232 rpm, it follows k... 


tan @ = 0.107 Nmax = 0.814 


The maximum temperature occurs at a point on the cylinder 
shifted in the direction of rotation, 66.96 deg away from the 
point facing the heat source. 

Increasing the rotational speed to v = 871 ft/hr = 0.242 
fps or 2.32 rpm corresponding to v% = 100, shifts the point 
of maximum temperature within a fraction of a degree to cos~! 
(1/2), and the maximum temperature can be approximated by 


Tmax = + 0.084 Po/Vo) 


Table 2 Positions of maximum temperature 
Up 3 1 V10 30 10 © 


Tmax is 1.008 T> for the particular case considered. 

For v = 1000 corresponding to 23.2 rpm with po = 9.12 
and the data otherwise unchanged, the temperature fluctua- 
tion is approximately 0.1 per cent of 7». 


Temperature Distribution for Intermediate Values of tr 


In these cases, the temperature distribution can be calcu- 
lated from Equations [24]. These calculations were per- 


formed for v = 1, J/ 10, 4/30 and 10, and the values are 
presented in Table 1 and Fig. 4. 


General Observations 


The largest temperature fluctuations occur at zero velocity. 
Because these fluctuations are of considerable magnitude, the 
error in the values of the temperature due to the process of 
linearization is of the order of 5.7 per cent for the numerical 
example chosen. 

As the speed increases, the temperature fluctuations be- 
come smaller. In this case the linearization of the difterential 
equation is fully justified, and the result can be considered to 
be accurate. 

Also in this case, the position of highest temperature is 
shifted in the direction of rotation. This shift is denoted by 
6. At very high speeds @ approaches cos—(1/m) approxi- 
mately 71.46 deg, and the fluctuations vanish. The position 
of the maximum temperature was computed for po = 9.12 
and a series of values of vp, by solution of differentiated Equa- 
tion [24c] set equal to zero. The values of 7max and the cor- 
responding angle of shift 6 are presented in Table 2 and Fig. 5. 

The minimum temperature occurs at zero speed on the far 
side of the cylinder, and, as the speed increases, the position 
shifts in the direction of rotation until at infinite speed the 
value of ymin, at which 7 minimum is attained, moves to 
cos! 
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Technical Notes 
Solar-Lunar Perturbations of the Orbit orbital angular momentum vector of the disturbing body. 


The orbit of the satellite is described by the instantaneous 
osculating ellipse (3) which has a semimajor axis a, an eccen- 


of an Earth Satellite 


MILDRED M. MOE! tricity e, an inclination 7 (with respect to the orbital plane of 
the disturbing body), an argument of perigee w, and a longi- 
Space Technology Laboratories, Inc., Los Angeles, Calif. tude of the ascending node Q (see Fig. 1). If there were no 
perturbing forces, these orbital elements would remain con- 
T# GRAVITATIONAL fields of the sun and moon can stant (assuming that our coordinate system is fixed in inertial 
produce significant perturbations of the orbit of a highly space), and the orbit would be a Kepler orbit (4). If there is 
eccentric Earth satellite. The importance of their effect a perturbing acceleration having a radial component R, a 
on the lifetime of 1959 delta (Explorer VI) has been pointed transverse component S (measured in the plane of the in- 
out by Kozai (1),? who has made a study (2) of solar-lunar stantaneous osculating ellipse and taken positive when mak- 
perturbations of Earth satellites. The purpose of this note ing an angle less than 90 deg with the velocity vector), and a 
is to treat these perturbations by a different method, based normal component W (normal to the plane of the instan- 
on some simplifying assumptions. It is assumed that the taneous osculating ellipse and taken positive when making an 
angular velocity of the disturbing body (sun or moon) is small angle less than 90 deg with the z axis), the orbital elements 
enough compared to the angular velocity of the satellite that change at the rates (5) 


we may consider the disturbing body fixed during one revolu- 
tion of the satellite. This simplification makes possible the pp ees 
integration of the instantaneous rate of change of the orbital dt nV1 -e 
elements over one revolution of the satellite to obtain the — ee ae 
change in orbital elements per revolution. The magnitude de =- Vi-—ésinv R Vi = Ss 
of the error made by holding the disturbing body fixed is dt na nave 
estimated. This leads to an evaluation of the effect of the : 

di _ cos (w+v) 


sin » 
da _ _ 2esine 2aV/1 
S 
nr 


motion of the disturbing body on the semimajor axis (hence Reo ae 
energy) of the satellite orbit. The importance of the solar- narv/1 — e? 7 a 
lunar perturbations of orbits with high eccentricity is illus- sn ( ai va 
trated by computing the perturbations of the perigee distance QQ ran(@ts) ro! 
of 1959 delta. dt na?V/1 — sini 
Perturbations of the Orbital Elements dw V1 
The analysis of the perturbations is simplified if we choose dt dt — 
a geocentric coordinate system having its z axis along the = 3 
Received Dec. 11, 1959. 
1 Member of the Technical Staff. 
where 
ayes? = n = satellite’s mean angular motion 
; r = radial distance of the satellite from Earth center 
v = true anomaly (the angular distance of the satellite from 
perigee ) 


The rate of change of perigee distance g can be found from its 
relationship to a and e 


q = a(l — e) 
The components R, S and W are given by Moulton (6). If 
r is always small compared to the distance aa of the disturbing 


body from Earth’s center, these components may be ex- 
panded in powers of r/az. Keeping only the first-order term, 


one obtains ? 
R = Kar (1 + 3 cos 29) _ = 


S = —6Kar [cos y sin (w + v) — 
x sin Y cos (w + v) cos i] cos 
Fig. 1 Satellite orbit W = —6Kar cos @sinisin y 


Epitor’s Note: The Technical Notes and Technical Comments sections of ARS JouRNAL are open to short manuscripts describing 
new developments or offering comments on papers previously published. Such manuscripts are usually published without editorial 
review within a few months of the date of receipt. Requirements as to style are the same as for regular contributions (see masthead page). 
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where we have used the subscript zero to indicate that the 
disturbing body is held fixed. Now let us evaluate Ag when 
the disturbing body is moving, by taking y = Yo + nat where 


Ka = GMa/2a,' 


G = universal gravitational constant na is the angular velocity of the disturbing body. Equation 

M,q = the mass of the disturbing body [2] for Aq will contain integrals of the form 

= angle between rand ag 

Y = angle from the line of nodes to ag : [= . f(cos v, sin v, cos 2y, sin 2y) dv [10] 


We make an additional simplification by taking az constant. 
Letting e stand for any orbital element and Ae for the 

change in that element after one revolution of the satellite — 

(from perigee to perigee), we have 


where f stands for some function. We replace cos 27 and 
sin 27 by their Taylor expansions 


t=2x/n de 2x de dt sin 2y = sin 2 9 s 
de = [2] sin 2Y in + 2nat cos + 
and retain only the first two terms. Our integral now con- 
where ¢ is time measured from perigee passage of the satellite. tains both v and t as variables of integration. These may be 
Since Ae is supposed to be small compared to ¢, it is permissible related to the eccentric anomaly £ via the relations 


to approximate all variables in the expressions [1] for de/dt 


by the values they would have in the unperturbed orbit, and = (E — esin E)/n 7 
to approximate dt/dv by its relationship to the conservation 
sin v sie — e? sin E)/(1 — e cos E) a 


of angular momentum h ie 
dor <a cos v = (cos E — e)/(1 — e cos E) _ 


dv ih Equation [10] now takes the form 
where h = V1 — e? is assumed constant. Since the E, sin E, cos E) dE 
angular velocity of the satellite is usually large compared to 0 pinata ‘ 
the angular velocity of the disturbing body, we may assume enated : 
that is constant during the time the satellite takes to com- pe of integr al 
plete one revolution. Then integrals of the type [2] can be 
= sin 2 
Aa = 0 [3] n? 
+ 45¢e + 58e? + 


15K V1 2 cot + === 
n? x 15e(1 + e) V1 —e 
(4) [12] 


{sin 27 cos 2w cos i — [cos* y — sin? y cos? i] sin 2w} 
which shows that the additional term, contributed by the 


Ae = — (1/a)A 5 
; oom (5) motion of the disturbing body, is of the order na/n which 
. sini sin y is usually less than 0.03. 

Ai = 4/1 — 4 Ve. % It is also of interest to find how much the semimajor axis 
is perturbed, and this can be found only by considering the 
eee : j isturbing body. We return to the three- 
motion of the distur ving 

(6) dimensional problem (i 0) and find, by a procedure similar 


to that used to obtain Equation [12], that 


Aa = {(4 — 6e — 5e?)[sin 2w cos 27 cos i — (1 + cos?i) sin 2Y cos? w] + [2(1 — + (2 — be — 3e?) cos? sin 270} 


Example 
AQ = —6Ka 7 sin y To illustrate the size of the perturbations, let us compute 
n?V/1 — e? x ,% their effect on the perigee distance of 1959 delta (Explorer 
VI). Considering first the lunar perturbation and assuming 
{ 5 ¢? sin 2w cos y + [1 + 4e? — 5e? cos* w] cos i sin vt [7] that the orbital elements do not change significantly over 
Aw = —AQcosi — E — 3 sin? y cos? i — > sin 2y sin 2w cos i + (1 — 5 cos? w)(cos? — sin? cos? [8] 


If, now, we let y take values from zero to 27, we see how the 
orbital elements change as a function of the time of month 
or year. 


the first few weeks of the satellite’s life, we may average 
Equation [4] over one lunar month to obtain for the average 
value of Ag 

Error Due to Neglecting the Motion of the Disturbing Body | 


» 
i i —15Km 
To estimate the error we have made by taking y constant ameV/1 — (13) 


let us consider the special case of co-planar motion of the 2n? 
satellite and the disturbing body (i = 0) and set w = 0. 
Then Equation [4] for the change in perigee per revolution 
becomes simply 


where the subscript m indicates that the disturbing body is 
the moon. For 1959 delta the initial values of the orbital 


s fo, elements were approximately a = 15,000 nautical miles, e = 
ar Aq = [(15KyameV/1 — e?)/n?] sin 270 [9] 0.76, im = 41.5 deg and w, = 10.4 deg. With these values, 
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we obtain Ag,, = —0.06 nautical mile per revolution. Re- 
membering that the satellite made about two revolutions a 
day, we see that the moon’s perturbation can lead to a non- 
negligible lowering of perigee. 

The sun’s influence, on the other hand, tends to raise 
perigee at first. Relative to the sun’s plane (or ecliptic), 
the argument of perigee started at 2.7 deg, the orbital in- 
clination was 40.3 deg, and y, was about 48 deg. The initial 
change in perigee, determined from Equation [4], is therefore 
Aq, = 0.26 nautical mile per revolution. At the end of one 
month Ag, = 0.14 nautical mile per revolution. Eventually 
the sun moves to such an angular position that Ag, becomes 
negative. Then the combined solar-lunar forces produce an 
important lowering of perigee with the consequent reduction 
of lifetime. To obtain accurate results, one must evaluate 
Equation [4] for each revolution, taking into account the 


‘ 4 Ibid., pp. 82-83. 
Ibid., pp. 404-405. 


changing values of %,, im, @s, @m, ete., dependent on all 
perturbations including those caused by Earth’s oblateness 
and atmosphere. 
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Jet Flows With Shocks 


MILDRED M. MOE! and B. ANDREAS TROESCH! 


Space Technology Laboratories, Inc., Los Angeles, Calif. 


HE INVESTIGATION of supersonic jet flows from 

nozzles is of considerable interest in rocket technology. 
A knowledge of the jet flow field furnishes important informa- 
tion for problems of radiation, heat transfer and the attenu- 
ation of electromagnetic waves by exhaust gases. 


ear 


line and symmetry axis, and the entropy S. Then the equa- 


; tions for the directions of the characteristic curves are 


dr/dz = tan (6 — yp) (right running characteristic) —_[1] 


dr/dz = tan (@ + yu) (left running characteristic ) [2] 


where 
r = distance from nozzle axis 
z = distance along the nozzle axis 


= Mach angle, sin—(1/M) 


The corresponding characteristic relations turn out to be — 


dz VM? (3) 


on 
a0 


In this note we use a method for computing the flow field - 
(including shocks) of an axially symmetric supersonic jet 
expanding from a nozzle. The nozzle is assumed to be in still 
air (the case of ambient pressure) or to move through the 
surrounding air at supersonic speed. The latter case is illus- 
trated in Fig. 1. The gas properties at the exit plane of the 
nozzle are assumed known. For an underexpanded nozzle, 
the initial slope of the jet boundary is determined by match- 
ing the shock equations in the air and the Prandtl-Meyer 
relations in the jet. Downstream from the exit plane the 
boundary turns inward, causing a compression resulting in the 
formation of a shock within the jet. The problem is to com- 
pute the entire flow field by the method of characteristics 
supplemented by the Rankine-Hugoniot relations for the 
shocks inside the jet flow and in the surrounding air. The 
present method cannot be used in regions where the flow 
becomes subsonic, e.g., behind the Mach disk. 


Basic Equations 

The computation is based on the equations governing the 
steady-state flow of a compressible, ideal gas. Viscous effects 
are neglected, and it is assumed that the changes in gas proper- 
ties occur through adiabatic processes; i.e., no heat transfer 
takes place. Furthermore it is assumed that the air and gas 
obey the entropic equation of state (1).2 It can be shown 
(1) that these equations represent, for supersonic flow, a 
system of hyperbolic partial differential equations. Therefore 
an integration can be carried out along the characteristic 
curves (or Mach lines) of the flow. For numerical integration 
it has proved convenient (2) to introduce as dependent vari- 
ables the Mach number M, the angle 6 between the stream- 
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dS 


— 1)M? 


where 
-Y = ratio of specific heats 
Cp = specific heat at constant pressure - 


The upper sign refers to the right running characteristic 
and the lower sign to the left running characteristic. It 
should be noted that pressure and density do not appear in 
these formulas, but may be computed from the Mach number 
and entropy. 


Numerical Procedure 

The computation consists of several distinct parts (see 
Fig. 1): The determination of the flow at the nozzle edge, 
involving the Prandtl-Meyer expansion and the exterior 
shock; the determination of the points on this attached 
shock; the points on the interior shock; the boundary points; 
and all the interior points. 

The computation of interior points is performed by the 
usual finite difference procedure (3,4), and all formulas are 
iterated. On the axis, instead of Equations [1 to 3], slightly 
modified formulas must be used (5,6). 

At the boundary between jet and air, which represents a 
contact discontinuity, the pressures and streamline angles 6 
must be equal on both sides. The boundary point 7 (see Fig. 
2) can be determined by satisfying these conditions along 
with the characteristic relations [3] for the left running 
characteristic through the fixed point 6 and the right running 
characteristic through point 5, which is interpolated between 
its neighbors and changes its position during the iteration 
process. 

At the first boundary point (point 3 in Fig. 2), a slightly 
special treatment is necessary. We assume that the attached 
shock has constant strength up to point 2, where the right 
running characteristic extended backward from the boundary 
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Fig. 1 Flow pattern of a supersonic jet in a supersonic stream 
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Fig. 2 Characteristic net for determining a jet flow 
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Fig. 3 Computation of a typical shock point. (Points denoted 
by squares have fixed positions and properties in the course of 
the iterations) 
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point 3 intersects the shock. This simple treatment is per- 
missible because, in actual computations, the first left running 
characteristic is chosen especially close to the nozzle edge. 

The succeeding points on the attached shock are determined 
by using the shock equations and Equation [3] along a left 
running characteristic. For example, in Fig. 2, the left 
running characteristic from point 3 intersects the shock at 
point 4 where the Rankine-Hugoniot shock conditions (1,3). 
relating the Mach numbers, flow angles and entropies on the 
two sides of the shock, are applied. 

Turning next to the interior shock, we note that the com- 
putation of this shock is more involved, since the approach 
flow is not uniform and the starting point of the shock is more 
difficult to locate. The determination of the shock position 
and the properties on both sides is, therefore, the most com- 
plicated part of the method [(4), p. 675]. A sketch of the 
points involved in a typical computation is given in Fig. 3; 
the overall picture of the converging right running character- 
istics leading to the shock is shown in Fig. 2. Given the shock 
point 4 in Fig. 3, the next shock point 3 is determined by 
first computing an auxiliary or virtual point 5. This is the 
point determined as the interior point from 1 and 4 under the 
assumption that there is no shock present. It is used to 
obtain the upstream properties at 3 by interpolation between 
points 1 and 5. Although the position of point 3 and the 
shock angle there are not known, they may be determined by 
satisfying simultaneously Equation [3] along the right 
running characteristic from 2 to 3, along the left running 
characteristic from 1 to 3, and the shock equations. Addi- 
tional information is given by the slope of the shock at 
point 4. 

In a finite difference approach, the first point on the interior 
shock is taken as the point where two like characteristics first 
intersect. It can easily be shown that for shocks in jets it is 
sufficient to test the distance from the axis of two consecutive 
points on a left running characteristic and insert a shock when 
this distance decreases, i.e., when a foldback occurs. 

The elements of the computational procedure described in 
the preceding paragraphs are combined into a complete pro- 
gram in a way easily explained by referring to Fig. 2. The 
flow properties given at the nozzle exit are stored on an initial 
line, either at the first Prandtl-Meyer characteristic for a 
bell-shaped nozzle, or a circular are for a. conical nozzle. 
The integration proceeds along a left running characteristic, 
from the initial line toward the boundary and exterior shock, 
then back to the initial line or axis for the next characteristic. 
When two right running characteristics intersect, the interior 
shock computation is started. The integration continues 
until eventually a normal shock (7) must be inserted, and, 
therefore, a subsonic region appears. 

It is important that the program be able to insert new 
characteristic curves whenever the characteristics diverge 
rapidly. This situation is illustrated in Fig. 2 by the spread- 
ing of the left running characteristics originating on the axis. 
The insertions (designated by crosses in Fig. 2) are made 
whenever the test for maximum allowable mesh size is not met. 

The procedure outlined here applies to the determination 
of a supersonic jet in a supersonic stream, but, with the appro- 
priate simplifications, can be used for jet flows expanding to 
ambient pressure. The details of the numerical method are 
described in (8). 


Results 


The numerical procedure has been used to compute many 
jet shapes, three of which are shown in Figs. 4, 5 and 6. The 
first of these is an example of a jet expanding to a constant 
pressure boundary and the second is an example of a supersonic 
jet in a supersonic stream. They should be compared with 
Figs. 4 and 10 of (2) which show the same jets as computed 
by the method of characteristics without the shock conditions. 
Since the position and strength of the interior shocks have 
been determined in our computations, a greater portion of 
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the jet flow field has been found. At first the shocks are very 
weak, but they become stronger as they bend toward the 
axis. For example, in Fig. 4 the pressure ratio across the shock 
increases from 1.6 at z = 4.0 to 4.6 at z = 8.8. In Fig. 5, it 
increases from 1.3 at z = 5.2 to 3.5 at z = 19.4. Fig. 6 is 
another case of a jet expanding to ambient pressure. It shows 
good agreement with the photograph of a jet (7) having the 
same exit and boundary conditions. - 


The authors are indebted to C. J. Wang, J. B. Peterson, 
H. M. Lieberstein and R. Kramer for valuable discussions. 
Of particular help on the treatment of the shocks were the 
discussions with J. M. Boyer (Convair, Astronautics, San 
Diego), who explained his method to us in full detail. W. A. 
Sassaman carried out the programming, joined in the final 
stage of the work by M. N. Perry and L. Stoller. 
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Specific Impulse 


WILLIAM L. SENN Jr.! and RONALD V. REPETTP 
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necessitates the repetition of a rather laborious computation. 
A method of performing these calculations quickly and accu- 
rately is therefore highly desirable. A graphical method 
satisfying these conditions has been devised which facilitates 
the rapid solution of the following equation 


A graphical method is presented which makes possible I 1 k [1] 
the rapid determination of specific impulse once the vari- P, 
ous parameters of a rocket propellant system are known. 
where 
Hy : specific impulse, sec 
HE EVALUATION of propellant systems often requires 
the calculation of specific impulse from the various sys- = acce 
tem parameters. Frequently, it is of interest to examine the 
dependence of specific impulse upon each parameter, which = = average specific heat ratio of exhaust gases, ¢,/c,_ 
ff = chamber temperature, R 
Received Dec. 8 1959 M = average molecular weight of exhaust gases, lb/lb 
eceived Dec. 8, 1959. 
1 Chemist. Present address, Esso Research Laboratories, mole 
Baton Rouge, La. <> P,, P; = chamber pressure and exhaust pressure, respectively, 
2 Chemical Engineer. Member ARS. consistent units 
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Because of certain restrictions imposed by Equation [1], 
the method presented here involves the use of two nomo- 
graphs. One graph uses the known parameters k and P;/P2 
in order to solve the expression 


® = 1 — [2] 


The quantity ® will be recognized as being equal to the ideal 
cycle efficiency of a constant pressure engine cycle operating 
between the pressures P; and P:. Once the value of ® has 


been determined, the second graph is used to compute J 
using as input data: k, ® and T/M. 

The error involved in solving Equation [1] by these nomo- 
graphs is less than 1 per cent of the computed value. 

The procedure for calculating J,, by this method is as 
follows: 


8p? 


Connect & with P;/P2; read off value of ®. 


2 Fig. 2: Connect k with ® found above and mark inter- 
section with reference line; connect reference line point with 
T/M; read off value of /,, 


Initial Azimuths and Times for 
Ballistic Lunar Impact Trajectories 


W. C. RIDDELL! 

Convair (Astronautics Diy.), San Diego, Calif. 

Analytic expressions have been derived that show ap- 
proximately the launch azimuths and launch times for 
which a ballistic missile will hit the moon. Also shown 
are the conditions producing the smallest angle of inter- 
section between the missile plane and the lunar plane, and 
those giving an optimum launch azimuth. The formulas 
are applicable to missions other than lunar impact. 


HE PURPOSE of this paper is to indicate the time of 

launch and initial azimuth for which an Earth launched 
missile will intersect the moon. The contact is a function of 
other parameters. However, they are, in turn, functions of 
the two specified parameters. 

The formulas derived here are approximate and are 
necessary for the determination of refined values for opera- 
tional purposes. They may be used with slight modification 
for the formulation of missions other than the moon inter- 
section. For example, the missile might be asked to contact 
a space station, in which case the formulas will give the time 
of launch and azimuth for such an event to occur. 

The derived formulas are the consequence of a simplified 
version of the complex physical structure. This system con- 
sists of the Earth-moon space and a ballistic body under the 
influence only of a potential field due to a spherical Earth. 

Such a model may appear to be oversimplified. How- 
ever, digital computer simulation in the complex env ironment 
gives results for a moon contact compatible with analytic 
solutions in the simple model. Tf, then, first approximations — 
are desired, the simplified version of this paper is adequate. 

The first part of the paper gives the pertinent usable formu- 
las, with brief discussion. The Appendix gives the deriva- 
tion of these formulas. Where multivalued functions are 
encountered, the proper or desired values are to be used. For 
definiteness, the launch site is considered to be in the northern 


hemisphere. 
Pertinent Formulas a= 


Roman numeral I refers to case I: Roman numeral 


II refers tocaseII: 1 <7, 
cos 6 = sin a sin7, cos H’ + 


cos @ (sin 7, sin / sin H’ + cosl) [1] 
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for an assigned H’ or an assigned a 


sin A” 
{1 + cos? A” [cos?/ tan? (7 — 1) — 2 cos/ sin tan }1/2 
sin H’ = tan l/tan 7, cos @ = cos 7./cos | 6 = 0 deg 
II [2b] 
cos @ = cos O(sin A”/sin [3] 
cos @ = cos A” cosl cos 6 + sin / sin 6 [3a] 
saat sin A’ = cot 7, tan 6 a [3b] 
=> 
Minimum |a| 
leos < |— sind = sinleos@ =Odeg I [3d] 
sin 
leos @| S$ | cos (1 — 
6 = -i, 
| 
0 2 |a| S 90 deg 
sin cos a = 0 deg II [3f] 


Epo cot @ = sin] tan H’ + (cos 1 cot 7./cos H’) [la] 


tan H’ = cot asin [1b] 
Absolute 6 minimum 


H’ = 90 deg a = O deg — 1, I [le] 


tan COS 
sin H’ = - cos @ = 6=Odeg II [ld] 
tan 7, cos | 
sin 7, 
cos § = (cos H’ cos l sin 6 — cos A’ sin] cos 6) + 
sin 


—— sin A” cosleos6 [2] 


Minimum @ for a given target point 


H’' = 90 deg 


cos 6 = cos a cos (1 — te) 


= 


V cos? (l — — cos? d/eos (1 — i.) sind 


tan 6 = tan/ cos A” a = 0 deg 


H 
3 
“we 
3 
: 
tan 
3g] 


= tanéd = tanlcosA” a = Odeg [3i] 


The launch site ion er coincides with the lunar plane. The 
least angle of intersection is always produced by launching at 
a time when the launch site is closest to the lunar plane. The 
least possible angle is given by the conditions of I in Equation 
[lc]. These apply when entering a parking orbit. The 
direct shot conditions are given by I in Equation [2a]. The 
target point is defined by ¢ or A”. Either of these param- 
eters, especially ¢, also partially defines the missile orbit. 

The conditions for nearest due east launch are given by © 
I in Equation [3d], I in Equation [8e] or by I in Equation — 
[3g], I in Equation [3h]. The launch azimuth and launch — 
time are found by substituting into the general equations — 
[3 to 3c]. The intersection angle is found by substituting — 
into Equation [1] or [2]. 


Case II 


The launch site passes through the lunar plane. A co- 
planar launch may always be made. The conditions are 
given by the identical sets II in Equation [1d] and IT in Equa- 
tion [2b]. Under the ideal condition of 1 = 7., the launch is 
due east. 7 
A due east launch may always be made under the conditions 
of II in Equation [3f], or II [3i]. The launch time is ob- 
tained from Equations [3a, 3b and 3c], and the intersection 
angle from Equation [1] or [2]. 


Appendix: Derivation of Formulas 


The unit position vector of the launch site is, in equatorial 
coordinates 


cos H’ cos | 
q =| sin H’ cosl 
sin / i 
The orthogonal north and east vectors are =. | 7 


—cos H’ sin l 
lity = | —sin H’ sinl 


cos | 


From these the azimuth vector is formed 


H’ sin 1 sin a@ — sin H’ cos 


[A-1] 


—sin H’ sin] sin a + cos H’ cosa [A-4] 
cos | sin @ 


This with Equation [A-1] determines the missile orbital plane 


sin a sin H’ — cos asin cos H’ 
[A-5] 
The plane of the moon is defined by its normal — - ; 


—(sin a cos H’ + cos asin! sin H 
0 
ly =| —sin 
COS 


cos @ cos l 
The angle of intersection between the planes is given by - 


In = 12° la = 


[A-6] 


Im-ly = cos 8 = sin a sini, cos H’ + 


cos (sin 7, sin sin H’ + cosi,cosl) [A-7] 
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A stationary value of this function will give a minimum value 
of 

Setting the partial derivative with respect to a equal to 
zero results in 


cot a = sin/ tan H’ + (cot 7, cos l/cos H’) [A-8] 
This gives the optimum a for an assigned H’. 
_ In like manner we obtain © 
tan H’ = sin! cota [A-9] 


Solving Equations [A-8 and A-9] simultaneously gives, for 
the absolute minimum 


cos @ = cos 7,/cos sin H’ = tan//tan i, [A-10] 


If 1 > 2., these solutions are imaginary. But a geometric 
consideration verifies that then the optimum remains the 
same asfor] =7,. Thatis 
H'=90deg 
~The unit position vector of the target is mag) 

( os A’ cos 6 


lp = in ow 8) [A-12] 


sin 6 
This with Equation [A-1] gives the missile orbital plane in 


a=0 


terms of launch position and target position ~~ 
"sin 
— . . . . 
- 1 sin H’ cos 1 sin 6 — sin A’ sin 1 cos 6 
- cos A’ sin 1 cos 6 — cos H’ cos 1 sin 6 [A-13] 
sin ¢ sin A” cos 1 cos 6 


If this is to be identical with Equation [A-5], then it is 
identical component by component. Equating the z com- 
ponents gives a relationship between the launch azimuth 
and the target point 

cos a = (sin A”/sin ¢) cos 6 [A-14] 
A” and ¢ are related through 4, as is found by taking 


1z-lr = cos @ = cos A” cosl cos 6 + sin 1 sin 6 


[A-15] 
so that [A-14] may be rewritten > 
cos a = 1 cos? 6 — (cos @ — sin sin 6)?/sin @ cos 
[A-16] 
or 
cos @ = cos 6 sin A "IV/1 — (cos A” cosl cos 6 + sin 1 sin 6)? 
[A-17] 


The target position is defined by 6. Its connection with A’ is, 


from spherical trigonometry 


sin A’ = cot 7, tan 6 [A-18] 
It is logical to inquire as to what target position, for a given 
¢ or A”, givesa minimum |a|. This is found by setting the 
partial derivation of Equation [A-16] or [A-17], with respect 
to 6, equal to zero. The conditions are, respectively 


sin 6 = sin! cos [A-19] 


tan 6 = tanl cos A” [A-20] 


In case 1 > 7,, there are obvious restrictions on the values of 
cos @ and cos A” for which the equations can be fulfilled. 
When the equations cannot be applied exactly, the value of 6 
coming nearest to giving equality is to be used. This will be 
for A” close to 180 deg and is 6 = 
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| cos A | <1 6= -i, 
0 < |a| 90 deg 
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from Equations [A-6 and A-13] 


sin 72, 


lm-lyu = cos 6 = 


This equation is difficult to use in minimizing 6, since H’, 
A’, A", @ and 6 are all interdependent. A simple form that 
can be analyzed by inspection is obtained by using lunar- 
Earth system coordinates. These are defined in the Nomen- 
clature. In this system 


cos u’ cos X 


sin u’ cos A [A-22 
cosu) 
0 
ly ={0 [A-24] 
lm cos u sin [A-25] 
sin sin cos X 
9 
= cos = cos u” cos X [A-26] 
= cos 6 = sin wu” cos A/sin [A-27] 
ry 
on cos 86 = V cos? — cos? d/sin [A-28] 


By inspection it is evident that the smallest |@| is given by 

vit | is g 
the smallest |X|. Sin @ cannot equal zero unless cos? \ = 1, 
in which case cos? 6 = 1. 


The restriction is imposed of - 
wey 


cos? @ S cos? [A-29] 


If Equation [A-1] is rotated into the lunar-Earth system, 
we have 
cos H’ cos 1 
={ sin H’ cos? + sin 1 sin 
—sin H’ cos / sin 7. + sin cos 7, 


[A-30] 


The intersection angle in terms of target position is obtained 


> (cos H’ cos | sin 6 — cos A’ sin 1 cos 6) + yen = 
si 


sin A” cos | cos 6 


ig 


Comparison with Equation [A-22] shows that 


[A-31] 


sin \ = —sin H’ sin 7. + sin 1 cos i. 


Thus J varies from 1 — i, tol + 7. as sin H’ varies from 1 to 
—1. Itisthe H’ = 90 deg, X = 1 — i. condition that is of 
interest when / > 7,., and the A = 0, sin H’ = tan 1 cot 7. condi- 
tion otherwise. 


Nomenclature 
A = right ascension of the target a 


H = sidereal hour angle of the launch site 

A’ = A — Q, right ascension of the target, measured from the 
ascending lunar equatorial node 

H’ = H-Q, 

A” = A’—H’'=A — H, arc of the Equator intercepted by the 
meridian planes of the launch site and target point. 
It is also the negative of the hour angle of the target 
from the launch site 

Ze = inclination of date of the lunar plane to the equatorial 
plane 

= latitude of the launch site 4 

= launch site 

MM == moonor referring to moon 

m = missile 

T = target, the point on the moon’s orbit intersected by the 
missile plane 

u = angular position or the target measured in the lunar plane 
from Qe. 

uv’ = launch site coordinate measured as above 

“ue 

a = launch azimuth, relative toa nonrotating Earth, measured 
from due east, positive to north 

6 = declination of the target 

@ = angle between the launch site position vector and the 
target position vector 

\ = launch site angular elevation from the lunar plane 

6 = angle of intersection between the plane of the missile 
and the plane of the moon . 

Qe = ascending equatorial node of the moon — 


Evaluation of Preflight Risks by Means 
of Very High Speed Digital System 
Simulation 

JESSE R. BRINKERHOFF! 


Boeing Airplane Co., Seattle, Wash. 


This paper gives a brief description of the mathematical 
and programming techniques employed in using a compre- 
hensive digital computer simulation of a complex system 
to provide, in a short time, accurate data on the chance of 
successful operation of the system in a specific test. A 
mathematical model of the system is obtained by writing 
equations and logical relations specifying the behavior of 
each part of the system as a function of the states of the 
other parts and parameters determined by external factors. 
These equations are then programmed on a computer 
and after setting parameter values and starting condi- 
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tions, the equations are solved and the solution is checked 
against an applicable criterion to determine whether per- 
formance is satisfactory. In practice, this procedure may 
become extremely complicated, and it is often abandoned 
in favor of other methods which analyze the various parts 
of the system separately or which make broad simplifying 
assumptions and generalizations to enable simple analyti- 
cal treatment of the system. 


HE SUBJECT of this paper is work done at the Boeing 

Airplane Co. in designing a simulation tool and using it 
for predicting the results of some individual test flights of the 
BOMARC interceptor. 

Performance of the BOMARC system is strongly dependent 
on the setting of guidance and control parameters, and final 
optimization of many of these parameter values was not com- 
plete for much of the flight testing schedule. Performance is 
also dependent on the conditions of the test, many of which 
could not be controlled in a tactical situation, but are subject 
to control in a flight test. Further uncertainties arose in 
early flight tests when performance of some parts of the sys- 
tem was temporarily restricted. Under these conditions, a 


[A-21] 
= 


! 


method of accurately predicting the chance of success for a 
particular test was desired for use in making decisions on 
parameter settings and test conditions. This tool was es- 
pecially needed since tests were contemplated with conditions 
which were near or beyond the nominally expected perform- 
ance limits of the system. 

For this purpose, and for use in parameter optimization 
studies and analysis of overall system performance in the 
composite tactical situation, a digital computer simulation, 
the BUG (BOMARC Univae AN/GPA-35), was designed. 
It was programmed in machine language for the Remington 
Rand Univae Scientific 1103-A computer, and AN/GPA-35 
is the type designation of the weapon control equipment 
utilized in the flight test program. The simulation had to 
satisfy two basic requirements: The system representation 
had to be complete and detailed enough so that results ob- 
tained would be accurate and dependable; it had to be fast 
enough to allow generation of statistically significant quanti- 
ties of data on short notice. 

The following paragraphs enumerate some of the major 
difficulties which were met in attempting to satisfy both re- 
quirements, the techniques which were utilized, the manner 
in which this simulation was used, and the type of results 
which were obtained. 


Development of the Model 


The BOMARC system can be roughly separated into three 
parts: 


1 Search radar and tracking. 
2 Guidance calculation and ground control. 
3 The missile, including homing guidance. 


The characteristics of each of these parts can be described, to 
a close approximation, by elementary mathematical equa- 
tions. However, the equations for one part depend on those 
of the other parts, and in many cases this interdependence is 
not simple. For example, there are some roughly periodic 
occurrences in each of the parts, e.g., rotation of the search 
radar, transmission of successive command signals and search- 
ing operation of the missile’s homing guidance device. The 


phase relationships between these quasi-periodic functions 
may be quite important, and the only information available 
about them is that they are probably random at some starting 
point of time and are determined thereafter by the relations 
between the quasi-periods of these occurrences. In the de- 
velopment of the model, these phase relationships are so rep- 
resented. Two more random items appear in the search radar 
and tracking simulation. The fact that each time the radar 
beam passes the target a return may or may not be obtained 
is taken into account in the model by requiring the specifica- 
tion of the nominal blip-to-scan ratio as an external parameter 
and making a random test each time to determine whether a 
return is obtained. In the same way, parameters are required 
for the standard deviation of radar noise, and a random noise 
error is ascribed to each return accordingly. 

Representation of the random processes in this manner 
poses the requirement that a statistically significant sample of 
runs is needed for a single set of input conditions in order to 
determine the probability of successful performance. How- 
ever, the estimate attained in this manner should accurately 
represent the actual probability. 

With the exception of the random sampling problems noted, 
the development of the model for the search radar and track- 
ing, and guidance calculation and ground control, sections of 
the system presented no major problems. With the missile 
model, however, two major problems appeared. 

The equations of this section were mainly the differential 
equations needed to specify the missile path as a function of 
time and of the operation of the airborne guidance and con- 
trol system. These had to be integrated numerically. For 
the greatest accuracy, a complex iterative integration scheme 
would be suggested, with a small mesh size. However, be- 
cause of the computer speed needed, the iterative integration 
scheme could be used only with a large mesh size, and it was 
decided to use a very simple integration scheme (Euler’s 
method) with a fairly small mesh size. This provided an in- 
tegration error smaller than the random errors in the system, 
and the greater accuracy which would be provided by a more 
complex integration method was not necessary. The small 
mesh size also proved advantageous in successfully represent- 
ing the relatively fine grained effects near the end of the in- 
tercept while the missile is under homing guidance. 

Another problem in developing the missile model was the 
fact that some quantities, for example, the thrust of the ram- 
jet engines, were given in graphical form and could be repre- 
sented mathematically only as arbitrary functions of several 
variables. One possible method of representing these func- 
tions in the computer would be to enter tables, and program 
a table read routine rather than a function evaluation routine. 
However, the tables needed would be too large for the high 
speed computer storage, and in order to achieve the computing 
speed required, the program had to operate almost exclusively 
from the high speed magnetic core storage. 

To alleviate this difficulty, a curve fitting process was de- 
veloped to provide simple polynomial approximations to these 
arbitrary data. Intuitively, the process appears somewhat 
unsound, and examples can be made up where it fails com- 
pletely, but in experience with actual data it seems that the 
process is remarkably versatile. It should be stated that the 
author’s experience with this process has been mainly with 
functions of only two variables. Although the extension of 
the method to functions of more than two variables is simple, 
it is not clear that the accuracy obtained would be satisfactory. 
An explanation of the procedure will be given here for a func- 
tion Z of two variables X and Y. 

Suppose the function Z is given on a graph of the form 
shown in Fig. 1. The first step in obtaining a polynomial 
approximation to the function is to find a first-degree (as 
shown) or second-degree (if more accuracy is needed) curve 
fit for each curve where Y equals constant. These curve fits 
are shown as dashed lines in Fig. 1, and have the equations 
shown there in terms of A; and B;,, the curve fit coefficients. 
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Fig. 2 Curve fits of zero-degree coefficients against second 
variable, Y 


This curve fitting can be done by least squares or any other 
method. 

The next step is to plot the zero-degree curve fit coefficients 
A; against the variable Y, as shown in Fig. 2. As before, a 
first- or second-degree curve fit is obtained for this curve, 
shown as the dashed line in Fig. 2, and its equation is given in 
terms of a;, the curve fit coefficients. 

In the same manner, the first-degree coefficients B; are 
plotted against the variable Y, as shown in Fig. 3, and the 
curve fit coefficients b; are obtained. The polynomial ap- 
proximation for Z is now completely specified, as 


A; B; 
Z= (do + (bo + bY )x 


In this explanatory example, first-degree curve fits were 
used throughout, but in the simulation problems both first- 
and second-degree fits were used. They can be intermixed at 
each stage of the process, and higher degree curve fits can be 
used if desired. If the process does not work too well, it may 
be improved by inverting the order of curve fitting, fitting 
first against the variable Y and then against X. In particu- 
larly bad cases the data may be broken into two or more 
classes, according to the value of one or both variables. 

Asan example of the applicability of this method, the actual — 
data shown in Fig. 4 were successfully represented in three — 
sections with only first- and second-degree fits. : 

Even with these simplifications, the missile model takes a — 
large proportion of the machine time, but the simulation is 
fast enough to enable a complete interception to be simulated 
in about 30 seconds. In contrast, the only BOMARC simula- 
tion available previously, programmed by Boeing on the — 
IBM 701, required about 30 min for the same job. 

The missile model developed for the BUG simulation has 
been used in essentially the same form for all the compre- | 
hensive digital simulations of the BOMARC IM-99A system — 
which have been written to date. Useful data were obtained © 
from the BUG simulation as early as June 1957. 


Use of Simulation for Preflight Risk Analysis 

To find the probability of success for specific test, the target 
paths, guidance and control parameters, and noise and blip-to- 
scan ratio parameters are inserted into the simulation. Then 
a number of runs are made, with only the random quantities 


ficult to determine, and it varies with conditions. For ex- 
ample, if the target is performing a periodic maneuver, the 
phase of the maneuver is an additional random variable, and 
more runs are needed. It has been determined from experi- 
ence that 20 or 30 runs will give a good estimate for a non- 
maneuvering target. 

After all the runs are obtained, the probability of success — 
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varying from run to run. The number of runs needed is dif- ; 
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Fig. 3 Curve fits of first-degree coefficients against second 
variable, Yo ¥ 
is just the ratio of the number of successful runs to the total 
number. 

Much additional data are provided by the simulation. Data 
are given from which correlation can be estimated between 
system performance and many variables of the system which 


© 


can be controlled by adjusting parameters. a 


Conclusions 


By using the technique described, it was possible to obtain 
a very high speed digital simulation as a tool for accurately 
evaluating the chance of success or failure of a flight test 
under particular conditions on short notice. This procedure 
was very useful in providing information for preflight de- 
cisions on test conditions and parameter settings. : 
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Prediction of Peak Temperature for 


Satellite Entries With Lift 


R. H. EDWARDS! and G. S. CAMPBELL? 


Hughes Aircraft Co., Culver City, Calif. 


LARGE share of the present understanding of lifting 

re-entry technique is based on numerical calculation. 
In order to avoid the necessity of repeating such computations 
with every change of parameters, it seems desirable to have 
a simple but accurate method of trajectory prediction. The 
approximations obtained might also be useful for developing, 
at a later time, lift programs to minimize equilibrium surface 
temperatures. 

The method used in the present note is essentially an 
extension of the skip trajectory analysis of (1). The 
formulas obtained give surface temperature, altitude and 
velocity at the bottom of the first skip for satellite entry at 
constant lift. Results of the approximate method are 
compared with IBM calculations from a design study re- 
ported in (2), 


The problem is to obtain, for given initi: al conditions at 
the edge of the atmosphere, an accurate estimate of the 
peak temperature during the first skip of a lifting satellite 
entry. Lift and drag coefficients are constant, and the 
flight path angle is assumed small. 


First approximation 
As a first approximation, velocity is constant and weight 
is balanced by centrifugal force, asin (1) 


(CiS/m)q — {1] 


The usual relations between density, altitude and flight path 
angle are assumed 


For constant velocity, the solution of these equations is 


« 


= —AXo tanh T 


where 
BVo Yo 
T = — — tanh — 
in i, 
2 PoC LS 
Bm 


Quantities with subseript zero denote initial values at t=0. 
After linearizing the relation for Xo, these equations are 


simplified 
Y = —Yo tanh T [4] 
Y= 
B cosh To ae [5] 
sech? T (6 
qo 
= + log (€/2) (7] 
= poCrS/2Bmy.2 [8] 
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The next approximation to the first skip trajectory requires 
an estimate of the velocity variation, which is obtained from 


the drag equation 
V = gy — (CpS/m)q [9] 


This equation may be integrated by substituting 
[4 and 6] to give 


4, = (V = Vo)/Vo [10] 
(y — yo) 
= 1 + tanh 7) [11] 


Vo? 

The first term on the right is the drag contribution to slow- 
down and is the only term in the skip analysis of (1). The 
remaining term represents the component of gravity in the 
drag direction. 


Second approximation 


The necessity for a higher approximation may be shown 
from the relation for equilibrium surface temperature given 


in (3) 


The peak temperature is achieved at the bottom of the skip 
at altitude (y:,, + hm), where y:,, is the first approximation 
to minimum altitude and h,, is a small correction. Peak 
density is therefore 


= — Bhm) 


from Equation [6] with 7 = 0. This relation, along with 
Equation [10], gives the approximate maximum equilibrium 
temperature 


Bm 3/4 1 3 
Tw=K Yo ''Vo Bhm + Wm | [13] 


The constant K can be readily evaluated from (3), for 
example. 

For representative satellite entry conditions, the secondary 
terms in Equation [13] provide temperature reductions up 
to 200 F. The necessary corrections are the first-order 
velocity term (Eq. [11] with 7 = 0) and the second approxi- 
mation to the minimum altitude. Development of the 
relation for Bh, follows. 

The complete lift equation, for small flight path inclina- 
tions, is 


The integral of this equation is, for the exponential atmos- 
phere 


CLS po 


7? = Yo? Bm — 


e —8y0) 


In the first approximation, denoted by subscript 1, the last 
term was neglected. The relation for the altitude correction 
hm is obtained by putting y = 0 in Equation [15], subtracting 
the first approximation, and assuming Bh» small i 


Entry altitudes are a small fraction of Earth’s radius, and 
so the gravity and radius variations may be linearized. For 
entry at velocities close to surface circular satellite velocity 


V let 


Yo" 


u = (Vo — Vgere)/Vo [17] 
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The integrand of Equation [16] finally becomes 
1 
r V2 


Comparison with exact numerical results has shown that the 
integration of Equation [16] may be carried out by replacing 
y and y,, by their first approximations. The resulting expres- 
sion for the altitude correction is 


1 C 
(“ -2 Yo + 2u) log 2e| [18] 
= @. e 4 


The difference between yo/r, and (ji, + Yo)/2r. was neglected 

in obtaining Equation [18]. The altitude at the bottom of 
. . 

the skip is the sum of y;,, and h,,, so that : 


Bym = Byo + log 2€ + Bhn [19] 
The velocity decrease to the bottom of the skip is 


Maximum stagnation temperature is calculated from 
Equation [13] with the proper values of hp» and Wy. 


Results 


The formulas from the preceding section have been used 
to compute peak temperatures during lifting satellite entries 
for a range of initial flight path angles. The results of the 
calculations are compared in Fig. 1 with IBM solutions 
obtained during the design study of (2). For entry angles 
above 1 deg, Equation [13] predicts peak temperature within 
10 F of the IBM solutions. An earlier method for peak skip 
temperature prediction (1) is subject to 180 F inaccuracy, 
and so a marked improvement is shown. 

The minimum altitude for the first skip is shown in Fig. 2, 
and velocity decrease is given in Fig. 3. Initial flight velocity 
Vo varied from 26,155 to 25,717 fps over the 0 to 4 deg range 
of yo because the initial conditions were established by single 
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Fig. 2. Minimum altitude reached during first skip of a lifting 
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Fig. 3 Velocity decrease during first skip of a lifting entry from 
60 nautical miles. 
C, = 0.54, Cp = 0.52, W/S = 15 lb/ft? 


retrothrust from a 300-nautical mile orbit. Equation [19] 
predicted minimum altitude within } mile for yo greater 
than $ deg. The predicted velocity decrease is quite in- 
accurate at the lower values of yo, but the effect of this 
inaccuracy on temperature is relatively small. 

One of the primary uses for the approximate skip trajectory 
equations is to estimate the effects of small parameter 
changes on temperature, dynamic pressure, ete. For 
example, differentiation of Equation [13] gives he effect 
of initial over-velocity on peak temperature 


oT w 2e 4 :) 
Br. Yo? 4 7 
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This formula confirms the numerical results of (2): That 
velocities greater than local satellite velocity can cause 
decreased stagnation temperatures. This effect is caused 
by increased centrifugal force with higher velocity. How- 
ever, it is seen from Equation [21] that for entry angles 
greater than 3 or 4 deg, depending on initial weight to lift 
ratio, stagnation temperature is increased by initial over- 
velocities. 
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Polymerization of Solid Propellants’ 
GROETZINGER 
i 
Rohm and Haas Co., Huntsville, Ala. 


EGLECTING internal heat transfer by convection, the 
problem of coupled heat transfer and polymerization 
during curing in simple, one-dimensional, solid propellant 

grains with constant thermal properties can be stated as 
a t) 0°7(s, t) Q 
Os Os? c 


at 
 aP(s, 0) 


P) 


with appropriate initial and boundary conditions, where 


G = geometric factor (0, 1 and 2 for slabs, cylinders and 
spheres, respectively ) 

7 = temperature 

P = polymer fraction (fraction of monomer converted to 
polymer) 

t = time 

s = distance from reference (centerline for cylinders, center 
for spheres and arbitrary for slabs) 

a = distance from reference to inside propellant surface (in- 


side radius for hollow cylinders and hollow spheres, 
zero for solid cylinders and solid spheres, arbitrary for 
slabs) 

b = distance from reference to outside propellant surface 
(outside radius for cylinders and spheres and arbi- 
trary for slabs) 

= propellant thermal diffusivity 

propellant specific heat 

propellant rate of polymerization function 

= total heat of polymerization per unit mass of propellant 


il 


This problem has been studied previously by investigators 
at Aerojet-General Corp. (1)* and Jet Propulsion Laboratory 
(3)* for solid cylindrical geometries. The Aerojet investiga- 
tors solved the problem numerically on a card calculator using 
explicit finite difference equations, whereas the JPL investi- 
gators solved the problem numerically on a differential ana- 
lyzer using difference-differential equations. The interest of 
these investigators was apparently temperature oriented. 
Results were reported principally in terms of centerline tem- 
perature. 

This problem was studied recently at Rohm and Haas Co. 
(7); the solution was achieved numerically on an electronic 
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digital computer using explicit finite difference methods. At- 
- tention was given to polymer fraction as well as to tempera- 
ture. The nature of the polymer fraction distribution is dis- 
cussed in this note. 

As an illustration, consider the typical problem of a hollow 
circular cylinder with an inside radius of 7 in. and an outside 
radius of 18 in. (1l-in. web, 0.85 loading fraction) initially 
filled with an unpolymerized propellant (P = 0) at 150 F, hav- 
ing a thermal diffusivity of 6(10~*) ft?/hr. The ratio of heat 
of polymerization to specific heat Q/c is taken to be 50 F, 

and the polymerization is assumed to progress at the rate 
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Fig. 3 Radial polymer fraction distribution. Hotter surfaces 


f(T, P) = BQ — P) exp (—A/RT), where the frequency fac- 
tor B is 10"/hr and the activation energy A is 3(10*) Btu/Ib 
mole. The boundaries are assumed to be held constant at 
150 F. 

The solution to this problem is shown in Figs. 1 and 2 which 
give the resulting temperature and polymer fraction distribu- 
tions, respectively. 

Taking this problem as a standard, the effect on the poly- 
mer fraction distribution of changes in wall temperature, 
ratio of heat of polymerization to specific heat, and tempera- 
ture sensitivity of the polymerization rate can be examined. 
Fig. 3 shows the effect of increasing the wall temperature 
25 F, Fig. 4 the effect of a similar decrease. Fig. 5 shows 
the effect of increasing the ratio of the heat of polymerization 
to specific heat by 25 F. Fig. 6 shows the effect of increasing 
the temperature sensitivity of the polymerization rate 
[achieved by increasing the activation energy from 3(10*) to 
4(104) Btu/lb mole while simultaneously increasing the fre- 
quency factor from 10" to 3.85 (10'*)/hr to keep the poly- 
merization time approximately the same]. 

The effect of temperature and polymer fraction distribution 
during curing on the quality of the finished grain is of im- 
portance to the propellant producer. The effect of tempera- 
ture gradients as well as the effects of thermal contraction 
and polymerization shrinkage on the stress-strain distribution 
in the grain-case combination under the assumptions of uni- 
form polymer fraction distribution and of elastic response 
upon the initiation of polymerization has been studied pre- 
viously (8,9). However the effect of nonuniform polymer 
fraction distributions, such as those illustrated here, on the 
stress-strain distribution both during curing and after the 
completion of curing presents an aspect not previously inves- 
tigated. Also, the nonuniform polymer fraction distributions 
presumably affect the strength properties of the polymerized 
propellant in critical areas, such as at the case bond and at the 
highly stressed inside propellant surface. 
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Real-Time Analysis, New Approach in 
Flight Testing 


GUENTHER HINTZE* 


A ‘White Sands Missile Range, N. Mex. 


The possibility of increasing the information contents of 
flight records through controlled exercise of missiles during 
firing tests is discussed. The discussion is based on a 
brief description of an experiment which is presently being 
planned at White Sands Missile Range for the accurate de- 
termination of aerodynamic coefficients under controlled 
and trimmed flight conditions. The best solutions for the 
analysis equations will be obtained when the exciting fin 
deflection and the corresponding angle of attack are ap- 
proximately 90 deg out of phase. This leads to the re- 
quirement to exercise the missile with a frequency close 
to its natural frequency. The basic relations of natural 
frequency and damping ratio as functions of missile ve- 
locity and altitude are derived, and typical trends of natural 
frequency with altitude for different Mach numbers, as 
obtained from flight simulation studies, are presented. 
Firing tests can be made into experiments yielding more 
information concerning the in-flight component and mis- 
sile behavior when through the real-time analysis of flight 
data exercises of the missile can be introduced which are 
adapted to the momentary flight condition. 


S THE first practical problem in connection with the 

development of real-time flight analysis and control at 
White Sands Missile Range, the measurement and determina- 
tion of certain aerodynamic coefficients have been chosen. 
The aerodynamic problem requires a minimum of change to 
the existing flight test program. The quality of present tele- 
metered data appears to be adequate for the proposed 
problem. The only special instrumentation required is an 
angle of attack meter and a triple linear potentiometer for 
sensing the fin position. Both of these have been used in En- 
gineer-User (E-U) flights. 

Presently, so-called aerodynamic rounds are flown to de- 
termine aerodynamic coefficients by flight tests whereby the 
missile is made to respond to pre-programmed commands with 
the guidance system kept inactive. However, this method 
does not exercise a closer control of the missile, keeping certain 
variables such as bank angle and undesired motion in several 
planes to a minimum. Real-time connection with the flight 
will control the missile into aerodynamic conditions which en- 
hance the measurement of specific aerodynamic coefficients 
by determining the type of maneuvers the missile must per- 
form to obtain the most useful data and by momentarily 
minimizing the effect of other interfering parameters. For 
the planned experiment a two-dimensional problem was 
chosen to simplify initial setup procedures and problem 
solution. 


Analysis Equations 


The experiment shall find a solution with respect te Mach 
number of the normal force coefficient 


Cy = Cyga + 


where 
a = angle of attack 
6 = fin deflection 


For simplicity the analysis is being restricted to the com- 
putation of the aerodynamic coefficients from the linearized 
normal force equation 


ma 


+ Cc = [2] 


Presented at the ARS Flight Testing Conference, March 23- 
25, 1959, Daytona Beach, Fla 

1 Chief, Flight Simulation Laboratory. 


where 


= mass of the missile 
= normal acceleration determined by an accelerometer 
= distance from the missile c.g. to the accelerometer 
= pitch angular acceleration 
= dynamic pressure 
S = reference area 
In Equation [2] all quantities except Cyvq and Cys can be ob- 
served directly or deduced from observed data. Thus, two 
sets of data observed at judicious times, t; and f, close enough 
to assume that Mach number has remained unchanged, will 
yield corresponding values of Cyq and Cys by solving the 
equations 
+ = 
+ 
The solution for the unknowns Cyg and Cyzg is found by 
Cramer’s rule 


> 


aN, 
Nob. lan Ne a6; | 
Cna A Cys A A= 252) [4] 


From Equations [4] it can be seen that good solutions for 
the two unknowns Cy, and Cys are obtained when the values 
for a and 6 have markedly changed between the observations 
at times 4; and tz, and when, in addition, a and 6 are offset by 
a distinct phase shift. Otherwise the determinant A may tend 
toward zero, resulting in singularity or near singularity of the 
matrix. On the other hand, it has been stated that the ob- 
servations at t; and #2 must be close enough to justify the as- 
sumption of a constant Mach number. 

The element of possible control in the experiment offers the 
opportunity for achieving data sets which will satisfy the above 
conditions. Real-time analysis of flight data provides the 
basis for a controlled exercise of the missile as compared with 
missiles which go into larger exercises only when some mal- 
functions occur. 


Derivation of Excitation Commands 


The excitation could be done either by step inputs or by a 
sinusoidal input for the commands on the fin deflection 6. 
Sinusoidal excitation has been chosen for several reasons. A 
step command introduces transients in the missile response 
which probably do not settle out in 3 sec for the missile 
under investigation, making very questionable the filtering 
process by which the values of the variables are averaged over 
the period for which the step command is being held. Sinus- 
oidal signals can be more easily and accurately smoothed to 
remove their noise content. In addition, the computation for 
transient response is more involved than the corresponding 
computation of the response of a linear system to steady-state 
sinusoidal excitation. With a sinusoidal input for 6 and as- 
suming a linear system, the quantities a, a and Q will also 
vary sinusoidally with the information of interest contained 
in their respective amplitudes and phases. 

In order to keep the controlled excitation as simple as pos- 
sible in this experiment, the missile will be exercised in climb 
only with the command 


G. = Go + G sin Wal [5] 
where 
Go = +0.5g is a fixed command to keep the missile on a long 


G, = amplitude 

@, = frequency of the sinusoidal excitation 

Additional safeguards will be provided in the real-time 
analysis experiment with the live missile, in case the planned 
exercise has to be changed back to normal missile control 


(1).? 


_ 2 Numbers in parentheses indicate References at end of paper. 


ARS JouRNAL 


+ high trajectory 


: 

| 
] 
| ( 

\ 

i 

1 

| 
| 

* 

» | 
\ 
: 


tory motion of the missile around its center of gravity. With 
Pe PERMISSIBLE the missile motion restricted to a vertical plane and consider- 
ACCELERATION ing the oscillation around the pitch axis only, the respective 

; moment equation is as follows 


SELECTED \ 
Max ACCELERATION 


Ip@ = qSc (cou [6] 


where, in addition to the previously defined quantities 


! = missile body attitude in pitch 
= Ip moment of inertia around pitch axis 


Fig. 1 Acceleration G,,,x vs. altitude = reference length 
aerodynamic moment coefficient with respect to @ 


ll 


Cus = 

For best resolution of the measured data the amplitude of Cm 6 2v aerodynamic moment coefficient with respect to 
the excitation should be maintained during the flight at the c6/2V 
maximum permissible level. For the test vehicle with which V = total missile velocity relative to air 
the experiment is being planned, the lateral acceleration is not Taking the Laplace transform and assuming for the 
allowed to go beyond a fixed number of G ~ determined by the momentarily considered oscillation a constant velocity V and 
altitude according to a curve acceleration Gnax vs. altitude R, initial conditions to be zero, Equation [6] can be written 
as shown in Fig. 1. 

With the exception of the rightmost segment, the computer gSc cs8(s) sa 
is programmed to select the total climb command G, accord- = E La) + [7] 
ing to the dotted line. This allows a small amount of re- 
serve in the command value so that the missile is not strained In order to obtain an oscillatory equation in @ only, the 
to its limits. The difference between solid and dotted line term L(a) has to be replaced by an expression in @. This 
varies from zero to about 3g. From the total command G,, expression can be derived from the normal force equation 


determined in this way, the constant upward command G) 5 . 6 (8] 
which compensates for gravity has to be subtracted in order to Binge ae — a) 
obtain the size of the sinusoidal exercising acceleration com- where no fin deflection is assumed, and 


mand G; sin wnt. 


The next step is the determination of the excitation fre- w = lateral acceleration in vertical plane 
quency w,. It has been stated that good solutions for u = longitudinal velocity component ~ V 
Equations [4] are obtained when the fin deflection 6 and the For small angles the angle of attack a can be expressed by 
angle of attack are clearly out of phase, the most desirable with collecting the with 
situation being a phase shift of 90 deg betweendanda. Fora Equation [8] a rearranged , 
second- order linear system, which can be assumed with suf- aca _ 
ficient accuracy for the tested missile, oscillations with the gS 6 
natural frequency w, will cause 6 and a to be 90 deg out of si — ' 7 [9] 
phase, with the sharpness of resonance determined by the 7 
damping ratio ¢ of the oscillating body. A universal chart Taking the Laplace transform 
which is applicable to any system with persistent oscillation, Vs6(s) (10) 

= V: 


illustrating these conditions (2), is shown in Fig. 2. 


Natural Frequency and Damping Ratio of Missile 
sO(s) 


w 


The natural frequency and of a are 


1.0 10.0 


Fig. 2 Phase plot with u = w/wn 
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Substituting this in Equation [7] gives 


= 


Cma*80(s) 


8 + (gS/mV ] 


Cm V 


Now letting 


c 
Sc gS gSc* ‘ 3 
= nm = = [13 
2 
the oscillatory equation for @ is obtained as follows 4 > “ee u 
fe] 
+ (6 + + (4 + 62)6 = 0 [14] 
with the natural frequency 
a 
w, = Va+ be = Cm, + eva) ~ 
[15] 1 1 j 
: : 40K 60K 80K 100K 120K 
and the damping ratio Ta rye ALTITUDE R_IN FEET ———~ 4 
z 
t= [16] Fig. 3 Natural frequency w,, vs. altitude 
_ From Equations [15 and 16] it can be seen that w, and ¢ a maximum of useful information without upsetting the flight. 


are mainly functions of the dynamic pressure g and missile 
velocity V; that means, since g = 3pV?, the natural frequency 
and the damping ratio of the missile vary with velocity and 
air density as a function of altitude. Typical trends of the 
variation of the natural frequency w, with altitude as obtained 
from a series of flight simulation studies are shown in Fig. 3. 

The curves can be expressed as powers of e. For easier 
computer solution the e-functions can be approximated by : 
power series for which only the first two terms have to be con- 
sidered when the curve is divided in small enough segments 
along the altitude scale. For each interval correction factors 
are predetermined and stored together with the known center 
values of altitude, Mach number and frequency for the inter- 
val. With these data, prepared in advance, it is easy to 
compute in real-time the natural frequency corresponding 
to the actual altitude and Mach number (3). 


Utilization of Real-Time Data 

The preceeding discussion brings out several points for the 
need of flight data analyzed in real-time. In order to obtain 
flight records with a high information content, missiles have 
to be exercised in flight. This excitation has to be done ac- 
cording to the momentary flight condition when it shall render 


Since missiles have a relative low damping ratio the phase 
shift around the natural frequency will be very steep, as shown 
in Fig. 2. Therefore, in order to obtain close to 90 deg phase 
shift between excitation and response it is necessary to excite 
the missile with a frequency very near to the natural fre- 
quency. According to Fig. 3 the natural frequency changes 
considerably with altitude, and a satisfactory phase shift can- 
not be obtained for all altitude and Mach number conditions 
using a single excitation frequency. 

Flight data, analyzed and displayed in real-time, will not 
only provide an immediate knowledge of the events and the 
behavior of the missile in flight, but they will also furnish auto- 
matically the proper inputs to adaptive excitation program- 
mers when such programs are being executed. 
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Ballistic Effect of Pyrolyzed Liner in 
Solid Propellant Motor Firings 


L. J. GORDON’ 


mek 
Aerojet-General Corp., Sacramento, Calif. 


Any combustion chamber liner material lost during 
burning of a solid propellant motor contributes to the 
momentum of the exhaust even though it adds no energy 
to the stream. A relationship yielding the true propellant 
specific impulse from the measured thrust-time integral, 
the propellant weight and the weight of liner lost in a static 
test firing is derived. A relationship between vehicle 
burnout velocity and these parameters is also given. 


Received Dec. 18, 1959. 
1 Senior Engineer, Thermodynamics Section, Solid Propellant 
R&D Laboratory. Member ARS. 


Test Firings 


N MANY cases the weight loss upon test firing of solid 
_& propellant motors is found to be greater than the weight 
_ of the propellant alone; the additional loss is due to pyrolysis 
or erosion of liner. It is necessary, however, to determine 
the specific impulse of the propellant from the experimental 
data even though some liner is lost. This discussion will 
show that the specific impulse of the propellant itself may be 
nearly correctly calculated from the experimentally obtained 
thrust-time integral and the weight loss by 


I, = FaT/(W,(W, + [1] 


where 
W, = weight of liner lost upon firing 
W, = propellant weight 


This result is first derived for a mathematically simplified 
case, and then a discussion of possibly more accurate methods 
of data reduction are discussed. 
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Assume for the simple case that: 

1 The liner is pyrolyzed at a uniform rate throughout the 
firing ; 

2 the specific heat and molecular weight of the propellant 
combustion products are equal to those of the liner decompo- 
sition products and are constant; 

3 there is no pressure drop during the pyrolysis and com- 
plete mixing of gases. 

A heat and material balance then yields 


W,C,(T; — T2) = Wig + W1C,(T2 — T'1) (2] 
where 
Cy, = specific heat of combustion products 
q = heat required for pyrolysis of unit weight of liner 
7; = combustion temperature of propellant 
7_ = initial temperature of liner 
7. = resultant temperature upon mixing of propellant and 


liner decomposition products 


Solving for yields 


* 


_= 3 
BI 


And since there has been no change in pressure, molecular 
weight, or in C, 


I, (7) 2 E + W1(T1/71) 


where 


I,, = specific impulse of propellant 
J, = specific impulse of mixed gases as expanded through 


nozzle 


This expression can be simplified by two approximations: 

1 Since < W, and T, < 7, then < W, 
so that W,(7./T,) can be neglected. 

2 If the heat of pyrolysis and interaction can be ignored 


I, - ( Ww, 5] 
I, \We+ Wi 


Experimentally one finds 


JS = (W, + [6] 


may be written ; 


1, = J 


Equation [1] may also be derived from simple energy- 
momentum considerations by realizing that addition of mass 


neglecting any heat of pyrolysis or interaction, a kinetic 
; energy balance yields 

=(W,+ Wt) — [8] 
ag 


where 
l V, = velocity that the propellant combustion products would 
have at the nozzle exhaust plane if no liner were pyro- 
lyzed 


V. = velocity the mixed products would have if the liner de- 
composition products obtained their velocity at the 
energy expense of the propellant gases 


One may then solve for J,, as follows (using the same 
nomenclature as in Equations [1 to 7]) 


‘ (Wet 
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and solving Equations [5 and 6] for the desired J,. yields — 
Equation [1]. For small values of W,/W,, Equation [1] | 


to the exhaust stream increases the momentum of the stream, © 
although the total kinetic energy remains constant. By — 


I, = + Wi) 


(W, + Wz) 


[10] 


(WW, + 


and Equation [11] is recognized to be identical to Equation 
[1]. 

It may be shown that this result is independent of the time 
period during which the liner is pyrolyzed. Assume, as a 
simple example, that the liner is pyrolyzed uniformly during 
the second half of the burning only, and that the propellant 
burns uniformly throughout. Applying Equation [1] to the 
second period and the normal thrust relationship to the first 
period 


T T/2 
FdT = FdT + FdT = 
0 0 T/2 
W, W, Wi 


and this result is identical to that in Equation [8]. It is then 
obvious that whatever the pyrolysis-time relationship, it can 
be broken into such periods. 

The only reasonable way to include the heat of pyrolysis 
when estimating the effect of pyrolyzed liner on a measured 
thrust is to assume that the decomposition products of the 
liner, which is of the same material as the propellant binder, 
come to equilibrium with the propellant combustion products. 
The gases that flow out the nozzle may then be considered to 
be identical to the combustion products of a composite pro- 
pellant having a greater binder proportion than the base 
propellant. The heats of pyrolysis and interaction are there- 
fore included in the specific impulse calculated for this com- 
posite propellant. In order to obtain the specific impulse of 
the base propellant from the thrust time integral obtained 
for the composite propellant, one must know the relative 
theoretical specific impulse and the combustion efficiencies of 
the base and composite propellants. For similar propellants, 
one generally expects combustion efficiencies for a given motor 
to be independent of propellant composition so that 


[13] 


where the primed numbers are theoretical values, and the un- 
primed numbers are experimental values. This yields the 
relationship 


JS Far 


[14] 


It is seen from the previous discussion that the usual ex- 


loaded vehicle, 


where A is the weight fraction expelled as liner and propellant. 
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Modulated Entry Into Earth’s 
Atmosphere 


JOHN V. BECKER! 
_ Langley Research Center, NASA, Langley Field, Va. 


N INTERESTING recent paper by Lees, Hartwig and 
Cohen (1)? discusses the use of modulated aerodynamic 
lift to reduce the peak resultant force experienced by a space 
vehicle entering the atmosphere. By reducing only the lift 
force in the low drag region of the force polar (Fig. 1) they ob- 
tained an attenuation factor (G/Gtixea tm) of about 4 for a 
vehicle capable of a lift-drag ratio of 3. It is eletenate 
that this mode of operation leads to difficult, if not impracti- 
cal, heating problems for the relatively low drag vehicles re- 
quired to achieve large attenuation factors. Since both the 
peak heat rate and the total convective heat load in laminar 
flow tend to increase as 1/WCp (where Cp is the average ef- 
fective drag coefficient), the problem encountered in this low- 
drag type of modulated entry is not readily dealt with by 
either the radiation cooled or heat absorbing structural ap- 
proaches. A recent study by Chapman (2) leads to this same 
conclusion. 

Grant (3,4) has investigated the modulation of both 
aerodynamic lift and drag according to realistic schedules corre- 
sponding to operation in various -regions of the force polar. 
He finds that maximum attenuation is achieved by an entry 
starting at Cz... and extending down the force polar into the 
low lift, low drag region (Fig. 1). The maximum attenuation 
ratio achieved with this schedule is considerably greater than 
for the “lift-only” mode, being about $ for a body capable of 
a lift-drag ratio of 3. If the region of the polar below maxi- 
mum lift-drag ratio is ruled out because of the heating diffi- 
culties, it is still possible to achieve an attenuation factor of 
about 3 by a schedule extending from Cimx to (L/D) max. 
The heating penalties associated with such a schedule are the 
subject of the remainder of this note. 

We consider first the case where force modulation is used to 
permit steeper entry path angles for a specific resultant force 
limitation of 10 g. It is assumed that all of the entries have 
the same initial velocity and the same velocity decrement in 
the pull-up maneuver. Using an unmodulated entry at L/D 
= 1 for reference, we obtain the comparisons for the total 
convective (laminar) heat load increment (AQ) given in 
Table 1. 

It is seen that the heat load penalty associated with modu- 
lation in the entry in which the low drag region of the polar 
is avoided is only about one quarter of the penalty found for 
the low drag entry. 

The increased heat flux rate in these modulated entries 
will, in general, cause no difficulty if the ablation technique is 
employed to absorb the heat pulse. However, if the heat 
pulse is absorbed by an internal coolant the increased heat 
rate may be an important factor in design. 

The guidance corridor widths corresponding to these various 
types of entry are shown in Fig. 2. (The “corridor” width 
Ah is defined as the difference in miles between the highest 
possible perigee for a single-pass entry at Cimx = —0.7 and 
the perigee of the approach trajectory corresponding to the 
10g entry.) It is evident that the modulation according to 
Grant (3) yields a corridor some three times wider than the 
fixed L/D entry, at the expense of some 60 per cent in the 
pull-up heat load. As indicated previously the “lift-only” 
modulation is less effective in widening the corridor and much 
more costly from the heating standpoint. 
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Fig. 1 Regions of force polar considered in modulated entry 


studies of (1 and 3) 
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Fig. 2. Guidance corridor widths for various modes of entry at 


parabolic velocity. 


Table 1 

Type of modulation i” 
fixed L/D = 1 

variable lift (1) 


(L/D) max CL = 0 
variable lift and drag (3) 
(L/D) max 


AQ/AQ(L/D=1) 
1.0 


° Table 2 
Type of entry Qn/AQn 1) 
fixed L/D = 1 5.69 1.0 
variable lift and 
drag (3), 
Cimax — (L/D) max 2.09 1.5 
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A second comparison has been made holding the entry angle 
as well as the other initial entry conditions constant, i.e., 
utilizing modulation to reduce the resultant force rather than 
to permit a steeper entry. In this study the vehicle surface 
was assumed to have a peak allowable temperature of 2500 F, 
and an emissivity of 0.9. The net laminar heat load (AQ,) 
which would have to be absorbed by the vehicle to hold the 
surface temperature at 2500 F was evaluated. The assumed 
entry conditions are: Velocity, 25,500 fps, entry angle 
—6 deg and altitude 400,000 ft. The vehicle had a flat- 
bottomed 70-deg swept delta wing with a blunt leading edge. 
Appropriate allowances were made for changes in distribution 
of the heating rates with angle of attack. The results ob- 
tained are shown in Table 2. 

The penalty in net heat load increments associated with 
modulation is seen to be about 50 per cent, or, relatively, 


the same moderate penalty as found in the previous example. 
The net heat load for the modulated entry in this case, in ab- 
solute terms, is only about 50 Btu per lb of vehicle weight for 


a vehicle in the 5000-lb class. It thus appears that modu- 
lated entries, properly executed, will not involve prohibitive 


a) = 


1 Lees, L., Hartwig, F. W. and Cohen, C. B., ‘‘Use of Aerodynamic Lift 
During Entry Into the Earth’s Atmosphere,’’ ARS JourNat, vol. 29, no. 9, 
1959, pp. 633-641. 

2 Chapman, D. R., ‘‘An Analysis of the Corridor and Guidance Require- 
ments for Supercircular Entry Into Planetary Atmospheres,’’ NASA TN 
D-136, 1959. 

3 Grant, F. C., ‘‘Importance of the Variation of Drag with Lift in Mini- 
mization of Satellite Entry Acceleration,” NASA TN D-120, Oct. 1959. 

4 Grant, F. C., ‘Analysis of Low-Acceleration Lifting Entry from Es- 
cape Speed,’” NASA TN D-249, 1960. 


References 


Electrically Charged Missile in Vertical 
Descent 


"pe 


MARVIN H. HEWITT! 


The Martin Co., Denver, Colo. 


The Lorentz deflection of an electrically charged missile 
in vertical descent from low altitude is determined and 
compared in nature with the Coriolis displacement in the 
northern hemisphere. It is found that a 10? kg missile 
having an initial altitude of 10 km at latitude 60 deg N and — 
an initial speed of Mach 10 is negligibly displaced. 


ECENT work? has indicated that small electrically 
charged missiles moving in the Earth’s magnetic field 
are perturbed and afford a small measure of control over 
the trajectory. Most work of this nature has necessitated 
machine computation or numerical methods. It will be 
shown here that perturbed vertical descent at low altitude 
can be determined without resort to a computer. The 
high Mach numbers now attainable and the fact that the 
Lorentz force is velocity-dependent elicit an interest in- 
ascertaining how effective this force is in producing de 
flection of a fast, relatively low altitude missile. 
The equation of motion governing the motion of a nega- 
tively charged missile in vertical descent over a rotating 
Earth with magnetic field is, neglecting drag 


ma = mg — 2mwxv—QxB 
where 
m = mass of missile t= 
@ = charge on missile in electromagnetic units 


= vertical velocity 


B = magnetic field strength due to terrestrial dipole 

9 = gravitational acceleration 
= = apparent acceleration of missile relative to Earth 

w 


= Earth’s rotational velocity, directed from south to north 
(dipole situated at Earth’s center) 
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The terrestrial field as treated here has components in a 
spherical coordinate system (r, 6, 


B, = B, sin 8 
Be = —(1/2)B, cos 0 [3] 
By = 0 [4] 
B, = 2M/r3 : [5] 


[2] 


where 


M = dipole moment = 8.1 X 10” gauss-cm* 
6 = latitude 
distance from Earth’s center to center of charged missile 


ll 


and 


B = + Bo? + By? = [6] 


where 
= V1 + 3 sin? 
A magnetic line of force has the equation 
r = ro cos? 8 \ = constant [8] 


where ro is the distance from the dipole center to where the 
line of force intersects the equatorial plane (9 = 0). The 
angle between the line of force and the radius vector is given 
by 


The direction of the velocity is at all times nearly down- 
ward. Consequently, the angle a can be treated as a con- 
stant. Because of its special interest, consider impact in 
the northern hemisphere where the Coriolis force gives an 
eastward deflection. Equation [1] becomes 


ma= ™g + 2mav cos OX — 


[9] 


Mve 
10) 


where X is a unit coordinate vector in the eastward direction. 
We see that the magnetic field causes a westward deflection 
of the missile in the northern hemisphere, for the velocity 
we have chosen. To a sufficient degree of approximation, 
the velocity is equal to the product of the gravitational 
acceleration by the time ¢ during which the missile has been 
falling, increased by the initial velocity ». We can now 
write the differential equation for the total deflection 


QM (gt + vo) cos 


mrs 


2w(gt + v9) cos — [11] 
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The height through which the body has fallen is very closely 
h = (1/2)gt?; consequently 


— (1/2)gt? [12] 


where r, is the mean Earth radius taken equal to 6371.221 
km. We consider those cases where r can be well approxi- 
mated by r,. Then, 


Mvp cos 


mr, 


Integrating Equation [13], we obtain the deflection 


3 6mr,3 2mr,3 


We recall that t = 2h g and therefore ae 
» (8h3\'/2 QM (8h3\'? 2h Mv 
r= (= cos 6 — (= cos 6 + avo (2 Q 
3\g9 6mr.3 \ g g 2mr. 


The westward deflection is 


om (ay\" Mos (23 It differs from the Coriolis displacement in being mass- and 
(*) (*)| on [16] charge-dependent, and, in the magnetic model assumed 


6mré \ g 2mré \ g here, is independent of Earth’s rotation. In addition, the 


Lorentz deflection is in an opposite direction, for negatively 
To illustrate the magnitude of A, choose a missile with a charged bodies, to the Coriolis displacement. cig : 
mass of 10? kg, an initial height of 10 km, at latitude 60 deg — 
N, and an initial velocity of Mach 10 with an eastward _ — 
Coriolis deflection. We find that A is less than a centimeter 
and therefore of no importance. 
It is interesting to note that the Lorentz deflection is The author expresses his gratitude to D. E. Muir for 
latitude-dependent in the same way as the Coriolis deflection. information concerning the problem. | 
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ROCKETRY LETTER SYMBOLS NOW AVAILABLE 


_ The American Society of Mechanical Engineers has just published the new Ameri- 
can Standard Letter Symbols for Rocket Propulsion, ASA Y10.14-1959, approved 
last September by the American Standards Assn. Sponsored by ASME, in collab- 
oration with the AmericAN Rocket Society and the Institute of the Aero- 
nautical Sciences, the standard letter symbols are available from the American 
Society of Mechanical Engineers, 29 West 39th St., New York 18, N.Y., at a price 
of $2 per copy, with discounts available. 
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New Patents 


. 
Sampling device (2,906,125). F.  B. 
Jewett Jr., New Canaan, Conn., assignor 
to the U.S. Navy. 

Deflated bag carried aloft by a balloon. 
The top of the bag is held open so the 
sample will enter and inflate the bag as it 
ascends. 

Nozzle (2,907,527). R. L. Cummings, 
Berea, Ohio, assignor to Thompson 
Ramo Wooldridge, Inc. 

Swirl-type fuel nozzle for directing a 
jet of pressurized fuel vapor through a 
chamber and flow orifice, to atomize and 
disperse the liquid fuel into a finely 
divided spray. 

Magnetic torque motor (2,908,171). L. 
Biberman, China Lake, Calif., assignor to 
the U. 8. Navy. 

Means for processing a gyroscope in a 
given direction at a predetermined rate. 
A hollow dome of magnetic material is 
mounted on the inner gimbal. Pairs of 
electromagnets produce a torque on the 
dome. 

Surface cooling means for aircraft (2,- 
908,455). H. H. Hoadley, Manchester, 
Conn., assignor to United Aircraft Corp. 

Porous surface exposed to an airstream 
of high supersonic Mach numbers subject- 
ing it to high temperatures. Cooling 


material exuded through the surface has 
high heat absorption characteristics. 
Guidance system (2,922,600). J.B. Craft, 
Los Angeles, Calif. 


Sensing elements responsive to the 
electrical conductivity of the surrounding 
air to control electrically operated controls. 
During movement of the vehicle, the ele- 
ments are disposed radially outward. 


Rocket motor (2,913,983). A. L. Lytle, 
Pacoima, Calif. 

Projectile with a releasable cover at one 
end, and a rocket motor in the other end. 
A partition plug carries pins to engage 
rifling in a gun barrel when launching. 
Safety and arming device (2,913,984). 

M. Meek, H. B. Smith Jr. and E. T 
Winston, Falls Church, Va., assignors to 
the U. S. Army. 

Remotely located rotor actuating means 
for use in a missile. A coil spring between 
a movable unit and the rotor holds apart 
mating slots when there is no acceleration. 


Eprtor’s Nore: Patents listed above 
were selected from the Official Gazette of 
the U. S. Patent Office. Printed copies 
of patents may be obtained from the 
Commissioner of Patents, Washington 25, 
D. C., at a cost of 25 cents each; design 
patents, 10 cents. = 
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Low noise electron gun (2,914,699). 
M. R. Currie, Beverly Hills, Calif., as- 
signor to Hughes Aircraft Co. 

Profile-shaping electrode at a potential 
positive with respect to an accelerating 
electrode. Emitted electrons are con- 
strained by a magnetic field to flow in a 
selected direction. 

Propyl nitrate monofuel (2,914,910). C. 
H. King Jr. (Member ARS), Jupiter, Fla., 
assignor to United Aircraft Corp. 

Method for decomposing monopropel- 

lant in a jet engine by adding a nickel salt 
soluble in propyl nitrate outside of the 
combustion chamber, converting the salt 
into finely divided particles. 
Controlling the effective thrust produced 
by ramjet units or athodyds (2,914,911). 
G. Richter, Nievre, France, assignor to 
SNECMA Co. 

Fuel valve responsive to the position of 

the inlet shock wave formed when the unit 
travels at supersonic speeds. Two static 
pressure taps in the intake are spaced at 
both sides of a selected location of the 
wave, and a third tap at a location where 
it is not affected by the wave. 
Apparatus for jet propulsion through water 
(2,914,913). Dr. F. Zwicky (Member 
ARS), Pasadena, Calif., assignor to Aero- 
jet-General Corp. 

Internally streamlined channel housing 

a propellant injector into which pres- 
surized water reactive propellant is re- 
leased. 
Three-dimensional variable exhaust 
nozzle (2,914,914). C. KE. Vandenberg, 
Fullerton, Calif., assignor to North Ameri- 
ean Aviation, Inc. 

Forward and aft elongated members 

pivotally mounted at one of their ends to 
form a jet engine exhaust tube. The 
other ends may be moved radially inward 
and outward when rotated. A longitu- 
dinal gap is caused between forward and 
aft members. 
Radiantly cooled inlet (2,914,915). E. A. 
Sziklas and C. M. Banas, Manchester, 
Conn., assignors to United Aircraft 
Corp. 

Opening and throat for conducting 
fluid at high supersonic velocities. A 
plurality of surfaces of the inner wall are 
parallel to the flow axis; alternating por- 
tions are at an angle to the axis. One set 
of surfaces of high emissivity faces toward 
the opening for radiating heat in that direc- 
tion. 

Control of fluid flow (2,914,916). P. Gelin 
and M. Kadosch, Paris, France, assignors 
to SNECMA Co. 

Two opposite auxiliary nozzles opening 
toward the interior of a jet propulsive 
nozzle and supplying pressure fluid trans- 
versely into the stream to form a partial 
obstruction to the flow and deflection of 
the thrust. 

Rotary jet engine (2,914,920). EK. J. 
Paradiso, Valley Stream, N. Y. 

Fuel and ignition carried by a fixed 
shaft to the combustion chamber within 
an elongated casing. Spaced nozzles 
project from the casing and open into the 
combustion chamber for receiving exhaust 
gases to rotate the casing. 

Rotating flame spreader (2,914,921). W. 
D. Pouchot (Member ARS), Media, Pa., 
assignor to Westinghouse Electric Corp. 


George F. McLaughlin, Contributor 


Blades moving in a direction transverse 
to the flow direction of high velocity air 
under pressure in a combustion chamber. 
Portions of blades are movable into and 
out of the wake of a stationary flameholder. 
High speed wind tunnel (2,914,941). 
O. Frenzl, Dammarie-les-Lys, France, 
assignor to SNECMA Co. 

Ejection-type tunnel. An accumulator 
contains water at a temperature greater 
than its boiling point at atmospheric 
pressure, and under a pressure sufficient to 
prevent boiling. The water passing 
through a nozzle partially vaporizes, ex- 
pands, and forms a fine spray. 

Variable sweepback airplane (2,922,601). 
B.N. Wallis, Effingham, England, assignor 
to Vickers-Armstrongs (Aircraft) Ltd. 


Delta fore-wing with tapered cantilever 
main wings pivotally connected at the 
rear. Means for driving main wings along 
an are to permit their movement in sweep 
during flight. 

Aerodyne with wings having variable 
sweepback (2,915,261). B. N. Wallis, 
England. 

Fore-wing of delta-shaped plan. Two 
pivoted main wings extending from the 
base ends are adjustable in sweep during 
supersonic and low subsonic flight. 
Rocket head design (186,798). L. Y. 
Booharin, San Francisco, Calif. 

Disk-like shroud below a conical head, 
and provided with four exhaust nozzles 
near its edge for imparting rotary motion 
to the vehicle. 

Supersonic airfoil (2,916,230). W. R. 
Nial, Troy, N. Y., assignor to General 
Electric Co. 

Airfoil having an upper and lower por- 
tion with opposing surfaces defining a 
diffusion channel. In the downstream 
direction, the channel has an inlet, a 
supersonic section of decreasing flow area, 
a throat section of minimum flow area, 
and a subsonic flow section of increasing 
flow area. 

Propellant powder (2,916,775). J. J. 
O’Neill (Member ARS), Alton, IIl., as- 
signor to Olin Mathieson Chemical Co. 

Pourable paste made by mixing particles 

of gelatinized nitrocellulose propellant 
and a liquid plasticizer. The composition 
is maintained at atmospheric pressure 
until the nitrocellulose is dissolved and 
the mass has solidified. 
Rocket powder grain (2,916,776). J. J. 
O’Neill (Member ARS) and C. E. Silk, 
Alton, Ill., assignor to Olin Mathieson 
Chemical Co. 

Plastic molding material subjected to 
changing temperature while confined in 
one dimension by a solid plastic sheath, 
and in the transverse dimension by a solid 
plastic member. All have the same 
thermal expansion coefficient. 
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High Temperature Materials, edited by 
R. F. Hehemann and G. M. Ault, John 
Wiley and Sons, New York, 1959, 
xvi + 544 pp. $17.50. 
Reviewed by LEon GREEN JR. 
Aeronutronic Division, Ford Motor Co. 


This book is a collection of papers pre- 
sented at an AIME Symposium in 1957. 
Part I, Cobalt and Nickel-Base Alloys, 
contains the chapters “Status and Future 
of Nickel-Base Alloys’ by F. L. Ver 
Snyder; ‘Investigation of High-Tempera- 
ture Properties of Nickel-Base Alloys 
Using Balanced-Experiment Designs” by 
T. L. Robertshaw; “Evaluation of a Cast 
Nickel-Base Alloy for Use as a Bucket 
Material at 1650° F’”’ by R. A. Signorelli, 
J. R. Johnston and J. W. Weeton; ‘W545 
—A New Higher-Temperature Turbine 
Disk Alloy” by J. T. Brown and J. Bulina; 
“History, Status, and Future of Cobalt 
Alloys” by G. A. Fritzlen; ‘“ML—1700 
Cast Turbine Bucket Alloy” by R. J. 
Morris and “Temperature Dependence of 
the Hardness of Secondary Phases Com- 
mon in Turbine Bucket Alloys” by J. H. 
Westbrook. 

Part II, Cermets and Intermetallics, 
contains the chapters, “Cermets and 
Bonded Hard Metals” by J. T. Norton; 
“Recent Developments in Sintered Ti- 
tanium Carbide Compositions’ by J. 
Wambold and J. C. Redmond; ‘Recent 
Advances in Infiltrated Titanium Car- 
bides’”” by H. W. Lavendel and C. G. 
Goetzel; ‘‘Wettability and Microstructure 
Studies in Liquid-Phase Sintering’ by 
N. M. Parikh and M. Humenik Jr.; 
“Modes of Fracture and Slip in Cemented 
Carbides”’ by N. M. Parikh; ‘Materials 
Based on Intermetallic Compounds” by 
R. Steinitz and ‘“‘The Use of Cermets as 
Gas-Turbine Blading,”’ by G. C. Deutsch. 

Part III, Refractory Metals, contains 
the chapters ““Chromium-Base Alloys for 
High-Temperature Applications’ by R. W. 
Fountain and M. Korchynsky; “Transi- 
tion Temperatures of Chromium and 
Chromium-Base Alloys’ by E. P. Abra- 
hamson II and N. J. Grant; ‘“Molyb- 
denum, Its Alloys and Its Protection” by 
J.J. Harwood and M.Semchyshen; ‘‘Proc- 
essing and Properties of Selected Experi- 
mental Molybdenum-Base Alloys’ by 
T. K. Redden and “Refractory Metals: 
Tungsten, Tantalum, Columbium, and 
Rhenium”’ by J. W. Pugh. 

Part IV, Strengthening by Dispersion 
of Insoluble Particles, contains the 
chapters ‘Observation of SAP and Present 
Theories for Its Remarkable High-Tem- 
perature Strength’ by F. V. Lenel and 
“The Development of SAP-Type Struc- 
tures in Titanium for Elevated Tempera- 
ture Service by Powder Metallurgical 
Techniques” by R. W. Jech, A. D. 
Schwope and E. P. Weber. 

Part V, Vacuum Melting and Its Effects 
on Properties, contains the chapters 


“Effects of Vacuum Melting on Gas Con- 
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Base Alloys’ by C. M. Hammond and 
R. A. Flinn; “Evaluation of Vacuum- 
Melted Vacuum-Investment-Cast Nickel- 
Base Alloys” by P. W. Beamer and J. J. 
Eisenhauer; ‘Consumable - Electrode 
Vacuum Remelting of High-Temperature 
Alloys” by R. K. Pitler, E. E. Reynolds 
and W. W. Dyrkacz; ‘Relations of High- 
Temperature Properties of a Nickel-Base 
Precipitation-Hardening Alloy to Con- 
tamination by Crucibles’”’ by R. F. Decker, 
J. P. Rowe and J. W. Freeman and ‘‘The 
Effect of Vacuum Melting on the High- 
Temperature Properties of Nickel-Base 
Superalloys” by F. M. Richmond. 

Part VI, Effect of Testing Environment 
on Properties, contains the chapters ‘‘The 
Effect of Environment on High-Tempera- 
ture Creep Properties of Metals and 
Alloys’ by D. A. Douglas; ““A Comparison 
of the Creep-Rupture Properties of Nickel 
in Air and in Vacuum”’ by P. Shaninian 
and M. R. Achter; “Effect of Environ- 
ments of Sodium Hydroxide, Air, and 
Argon on the Stress-Rupture Properties 
of Nickel at 1500° F’”’ by H. T. McHenry 
and H. B. Probst and ‘‘The Effect of Test- 
ing Atmospheres on the Creep-Rupture 
Properties of Molybdenum-Base Alloys at 
1800° F” by D. D. Lawthers and M. J. 
Manjoine. 

Part VII, Oxidation Resistance, con- 
tains the chapters “‘Oxidation Resistance 
of Binary Chromium-Base Alloys” by E. 
P. Abrahamson II and N. J. Grant; ‘‘The 
Oxidation of Tungsten and Molybdenum 
from 1800 to 2500° F’’ by J. W. Semmel 
Jr. and “High-Temperature Oxidation of 
Gas-Turbine Alloys” by J. F. Radvich. 

As is evident from the contents, the 
symposium was concerned with materials 
for a specific application, namely, air- 
craft gas turbines, and a large part of the 
book consists of a collection of metallurgi- 
cal and engineering data of use to turbine 
designers. This information on specific 
alloys or cermets usable up to about 2000 F 
will be of interest to readers in the rocket 
field who are concerned with the design 
of turbine drives for propellant pumps, 
auxiliary power supplies and the like. To 
readers concerned with thermal protection 
against re-entry heating or with the design 
of nozzles and controls for large solid pro- 
pellant rockets, on the other hand, the 
title of this book will appear much too 
general. In this sense the title is doubly 
misleading, since in general (with the ex- 
ception of the chapter by Pugh and pos- 
sibly that by Harwood and Semchyshen) 
the temperature levels involved might 
more properly be classed as ‘‘intermedi- 
ate,’”’ and since the refractory nonmetals 
were not considered except as constituents 
in cermets or dispersion-hardened metals. 

Despite these shortcomings, the book 
will be of general interest because of the 
presence of chapters which are general or 
basic in character, treating principles 
which may guide the development of new 


tent and Metal Structure of Cast Nickel- 


high temperature materials. The discus- 
sions by Westbrook, Norton, Parikh and 
Humenik, Parikh and Lenel are of this 
nature. Also, the interpretations by 
Douglas of specific data on Inconel may 
have general validity. It is believed that 
this book will be desired by institutional 
libraries as a reference work, but because 
of its high price, its appearance in personal 
libraries will probably be limited. 


Soviet Space Science, by Ari Shternfeld 
(translated from the Russian by the 
Technical Documents Liaison Office, 
WPAFB), Basic Books, Inc., New York, 
1959, xviii + 349 pp. $6. 

Reviewed by C. C. MigssE 
Armour Research Foundation 


This book, recently translated from the G 
Russian by the United States Air Force, 
provides a broad general summary of 
international progress in space technology 
for the novice, and a compendium of 
astronautical facts and figures for the 
lecturer. The author displays his math- 
ematical genius for the calculation of 
orbits, orbital periods and _ trajectories, 
and the reader is presented with a magnifi- 
cent display of mathematical space-gems 
which illustrate the unexpected phenom- 
ena associated with optimal trajectories. 
Although repeated reference is made to the 
characteristics and problems associated 
with the first two Sputniki, and a re- 
vealing review of recent Russian achieve- 
ments in space technology is presented to 
the reader, no undue emphasis is placed on 
Soviet plans or progress. Hence, the title 
ascribed to the translated edition is some- 
what misleading, and should probably be 
replaced by a direct translation of the 
original Russian title: ‘Artificial Satel- 
lites.” 

The opening chapters present the laws of 
celestial mechanics as applied to artificial 
satellites. There follows next a descrip- 
tion of the rocket propulsion system and 
launching procedures, liberally augmented 
by the mathematical details of rocket 
staging and optimum launching charac- 
teristics. Brief discussions of satellite 
construction and the survival problems to 
be encountered by the first cosmonauts are 
followed by careful analyses of the ob- 
servations possible from a satellite plat- 
form, and of the mechanics of re-entry 
and nongeocentric _ satellites. After 
several chapters on the problems of com- 
munication, and the utilization of artificial 
satellites, the author concludes with some 
consideration of international cooperation 
and space law. 

Perhaps the best critique of the treatise 
is provided by the author himself, in his 
Preface to the Second Soviet Edition: 
‘Within the scope of this book it has not 
been possible to discuss the entire complex 
of questions connected with artificial 
satellites. We have given major space to 
questions of celestial mechanics and 
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Maintenance 


Multi-Use 
Automated 


The recent demonstration of multi-purpose 
test equipment (MPTE), developed by 
RCA under a series of Army Ordnance con- 
tracts, highlights a new dimension in auto- 
mated multi-use systems support and culmi- 
nates a long-term RCA effort in this field. 
This General Evaluation Equipment is an 
automated, transistorized, dynamic check- 
out system. It contains a completely modu- 
~ larized array of electronic and mechanical 


Tmk(s) ® 


evaluation equipment, capable of checking 
a variety of electromechanical devices, 
ranging from radar subassemblies to missile 
guidance computers. MPTE provides the 
stimuli, programming, control, measure- 
ment and test functions for the NIKE AJAX, 
NIKE HERCULES, LACROSSE, HAWK 
and CORPORAL missile systems and has 
been extended to other weapons systems 
related to our defense efforts. _ 
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rocketry, as the basic elements of the road 
toward realization of artificial satellites, 
and have relegated to secondary place 
questions such as remote control, radio 
communications, and physiology. This 
book presents not only a discussion of 
propositions already known but also re- 
sults of the author’s own research.” 


Explosions, Detonations, Flammability 
and Ignition, AGARDograph 31, by 
B. P. Mullins and 8S. S. Penner, Per- 
gamon Press, London, 1959, 287 pp. 
$10. 

Reviewed by Ropert A. Gross 
University of California 

This single volume really comprises two 
books. The first, by Penner, is a 107-page 
summary of the current analytical under- 
standing of aerothermochemistry. The 
second, by Mullins, is a 162-page summary 
of experimental findings in the area of 
flammability, ignition and _ explosion. 
There is, unfortunately, little interaction 
between the two parts other than that 
they are examining, from different view- 
points, the same subject. This is not 
surprising since it is indicative of the state 
of affairs in this complex subjec and, 
indeed, is one of the values of the book. 
The subject matter treated concerns ad- 
vances in the subject made recently, but 
prior to 1958. The reader is expected to 
be familiar with the classic treatments and 
to have access to or previous knowledge of 
AGARDographs on related subjects. 

In Part I are presented the funda- 
mentals of aerothermochemistry, explo- 
sion limit theory for premixed gases with 
and without wall effects, detonation 
theory, flammability, quenching, ignition 
and single droplet burning theory. As 
befits analysis, the problems are simplified 
and models are frequently employed. The 
list of references is very good and includes 
some of the important Russian progress in 
this field. 

In Part II data on flammability and igni- 
tion are critically reviewed and correlated. 
This is done using a simple thermal ap- 
proach with no chain reaction concepts 
employed. Thus, limiting combustion 
phenomena are thought to occur when the 
rate of heat developed by chemical trans- 
formation equals the rate of heat lost to 
the surroundings. This simple concept 
frequently (but not always) gives good 
correlation for flammability limits of 
homogeneous and _ heterogeneous 
oxidizer mixtures, spark ignition, and 
spontaneous ignition phenomena.  Fi- 
nally, two very good and welcome chapters 
are presented on explosion prevention with 
particular emphasis on techniques used in 
large engine development test facilities. 
These final chapters will be helpful to 
those persons who are responsible for the 
running and safety of large test facilities. 

Unfortunately the authors do not give 
a separate listing of symbols and for this 
subject, which contains so many disci- 
plines, the notation is often complex and 
confusing. Since this book will serve 
primarily as a reference rather than a text 
book (to be studied from cover to cover), 
this lack of a listing of notation is par- 
ticularly annoying. It struck this re- 
viewer as surprising that, in this subject 
which is so concerned with rate phenom- 
ena, there is no mention made, either 
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theoretically or experimentally, of collision 
cross sections. Neither is there mention 
made of the microscopic point of view to 
aid in understanding the macroscopic 
events. In general this book will appeal 
to, and be used by, persons active in this 
type of research or those who wish to 
apply the results of such research. 


Handbook of Supersonic Aerodynamics, 
Vol. 3, Section 6, Two-Dimensional 
Airfoils, prepared by the Applied 
Physics Laboratory, Johns Hopkins 
University and the Defense Research 
Laboratory, University of Texas. 
NavOrd Report 1488 (Vol. 3), 1957, iv 
+96 pp. $1.50. 

Reviewed by Arruur A. Kovrrz 
Northwestern University 


This section is another in the impressive 
series, published under cognizance of the 
Department of the Navy, Bureau of 
Ordnance, titled Handbook of Supersonic 
Aerodynamics. 

The object of Section 6 is to provide the 
designer with data pertaining to certain 
two-dimensional airfoils. Included are 
data in graphical form for double sym- 
metrical sharp edge airfoil sections of 


double wedge, flattened double wedge and 
biconvex types. Additional data are 
provided for modification of the types 
mentioned, including the effects of camber 
and position of maximum thickness. Fi- 
nally, data for estimating the effect of vis- 
cosity on drag coefficients are also pre- 
sented. 

As a background for the data the 
theories applicable to these airfoils are 
given in short résumé. Comparison of the 
theories with available experimental data 
is also presented in graphical form. 

This section is designed to serve as an 
introduction to Section 7, Three-Dimen- 
sional Airfoils, which treats more practical 
design considerations. As such the re- 
viewer feels it has done its job very well. 
The theoretical resume has sufficient de- 
tail to adequately review and renew the 
reader’s knowledge of two-dimensional 
supersonic flow. At the same time it is 
not so lengthy as to cloud the main pur- 
pose of the section, i.e., presentation of 
data. 

The data is presented on clearly labeled 
graphs that can stand apart from the text. 
This adds greatly to the usefulness of the 
section. 

This reviewer finds it a source of com- 


fort to know that at his fingertips is a 
readily accessible, easily readable hand- 
book of two-dimensional supersonic air- 
foil theory. 


Book Notices 


Magnetic and Electrical Fundamentals 


(Franklinian Approach), by Alexander 
Efron, John F. Rider, Publisher, Inc., 
New York, 1959, 132 pp. $2.50. 


This booklet is written for the high 
school and the junior college level. It is 
profusely illustrated and descriptive in na- 
ture. 

The Philosophy of Space and Time, by 
Hans Reichenbach, translated by M. 
Reichenbach and J. Freund, Dover 
Publications, New York, 1958, 295 pp. $2. 

This is the English translation from the 
German of ‘Philosophie der Raum-Zeit 
Lehre.” 

Fuels and Lubricants, by M. Popovich 
and C. Hering, John Wiley & Sons, New 
York, 1959, 312 pp. $8.50. 

This volume constitutes a compilation 
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